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CHAPTER  ONE 
INTRODUCTION 

Importance  of  III~V  Semiconducting  Materials 
The  development  of  the  solid  state  electronics  industry  has  princi¬ 
pally  centered  around  the  semiconductor  material  Si  due,  primarily,  to 
its  good  electrical  properties,  available  high  purity  and  the  relatively 
simple  chemistry  surrounding  the  fabrication  of  silicon  semiconductor 
devices.  The  rapid  growth  experienced  by  the  semiconductor  industry  has 
imposed  further  demands  on  device  characteristics,  such  is  extremely 
high  speed  and  optical  emission  properties,  which  cannot  be  met  by  silicon 
as  a  result  of  its  relatively  low  electron  mobility  and  fixed  indirect 
bandgap.  The  desire  to  satisfy  these  new  demands  has  sparred  the 
development  of  compound  semiconductors  composed  of  the  elements  Al ,  Ga 
and  In  from  column  Illb  and  P,  As  and  Sb  from  column  Vb  of  the  periodic 
table.  Semiconductors  formed  from  combinations  of  these  elements 
provide  a  large  selection  of  electrical  and  physical  properties  such 
as  electron  mobility,  bandgap  energy,  possible  direct  bandgap  and  crystal 
lattice  parameter.  Table  1-1  compares  these  parameters  for  all  of  the 


III-V  binary  semiconductors  and  silicon.  Electron  mobilities  ranging 
from  80  to  100000  cm^/V-s  and  bandgap  energies  ranging  from  0.18  to 
2.45  eV  are  available  using  the  binary  III-V  semiconductors.  This 
tremendous  range  of  electrical  properties  contrasts  sharply  with  the 


single  available  set  of  properties  provided  by  silicon. 
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Semiconductors  consisting  of  eianents  from  the  group  III  and  V 
columns  of  the  periodic  table  need  not  be  limited  to  simple  binary 
compounds.  Indeed,  the  use  of  ternary  and  quaternary  solid  solutions 
provides  a  means  of  varying  the  physical  and  electrical  properties  of 
III-V  materials  continuously  between  the  limits  in  Table  1-1.  The 
vastly  increased  electron  mobility,  which  is  characteristic  of  many 
of  the  III-V  semiconductors  over  that  of  silicon,  makes  these  materials 
ideally  suited  to  advancing  the  current  capabilities  of  high  speed  de¬ 
vices  (e.g.  oscillators  and  mixers  in  the  100  GHZ  range,  central  pro¬ 
cessors  and  computer  memories). 

Currently,  the  most  important  application  for  semiconductors  made 
from  III-V  materials  is  in  the  production  of  optoelectronic  devices. 

The  energy  band  structure  for  many  III-V  materials  is  of  a  direct 
nature;  that  is,  a  third  particle  (e.g.  phonon)  is  unnecessary  for  the 
generation  and  recombination  of  free  carriers.  This  allows  for  the 
efficient  conversion  of  electromagnetic  radiation  into  electrical  pov/er 
(photovoltaic  cells)  or,  in  the  opposite  situation,  emission  of  radia¬ 
tion  from  the  semiconductor  device  (light  emitting  diodes,  lasers). 

The  available  bandgap  energies  for  binary  III-V  devices  having  a 
direct  band  structure  result  in  device  emission  characteristics  which 
range  from  the  infrared  into  the  visible  region  of  the  electromagnetic 
radiation  spectrum. 

Figure  1-1  plots  the  lattice  parameters  of  many  possible  III-V 
compounds  against  the  bandgap  energies  of  these  compounds  at  300K. 

The  solid  dots  represent  binary  compounds  and  the  lines  connecting  each 
dot  represent  ternary  solid  solutions  of  intermediate  composition  be¬ 
tween  the  binaries.  Solid  lines  signify  direct  bandgap  materials 
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while  broken  Tines  indicate  indirect  materials.  Essentially,  the  entire 
area  enclosed  in  Figure  1-1  is  accessible  to  the  designer  when  employing 
ternary  and  quaternary  III-V  solid  solutions.  This  flexibility  is  ex¬ 
tremely  useful  for  the  optimal  design  of  new  solid-state  electronic 
devices. 

Currently,  the  quaternary  systan  is  receiving  much 

attention.  The  available  ranges  for  the  lattice  parameter  and  bandgap 
energy  in  this  system  are  given  by  the  cross  hatched  surface  shown  in 
Figure  1-1.  Thus,  the  lattice  parameter  and  bandgap  energy  may  be 
specified  independently  with  the  composition  of  the  quaternary  solution 
chosen  to  meet  these  specifications.  The  availability  of  a  degree  of 
freedom  in  the  lattice  parameter  is  extremely  important  since  currently 
only  GaAs,  GaSb,  GaP,  InP  and  InSb  are  available  in  bulk  crystal  form 
for  use  as  substrate  materials.  A  disparity  of  greater  than  O.IS  be¬ 
tween  the  lattice  parameter  of  the  substrate  and  epitaxial  layer  induces 
the  formation  of  interface  defects  in  the  crystal  structure  which  can 
degrade  the  device  performance.  One  important  application  of  the  qua¬ 
ternary  ’s  in  the  development  of  heterojunction  lasers 

for  use  as  transmitters  in  optical  fiber  communication  systems  [2]. 
Currently,  available  optical  fibers  exhibit  minima  in  attenuation  and 
dispersion  characteristics  for  radiation  of  approximately  1  eV  [2]. 
Choosing  the  values  of  x=0.8  and  y=0.35  yields  an  emission  energy  of 
-1.11  eV  [3]  with  a  lattice  parameter  which  closely  matches  that  of 
the  InP  substrate. 


Epitaxy  of  III-V  Semiconducting  Materials 
There  currently  exists  three  primary  methods  for  growing  epitaxii 
III-V  films:  liquid  phase  epitaxy  (LPE),  molecular  beam  epitaxy  (MBE 


and  chemical  vapor  deposition  (CVO).  LPE  is  the  growth  of  thin  single 
crystal  layers  from  a  liquid  solution.  The  driving  force  for  film 
deposition  is  most  commonly  provided  by  decreasing  the  temperature  of 
the  substrate  relative  to  that  of  the  liquid.  The  advantages  of  LPE 
include: 

1)  The  method  is  capable  of  growing  multicomponent  layers  with 
a  high  reactivity  disparity  among  the  elements. 

2)  The  equipment  is  relatively  simple  and  inexpensive. 

3)  A  large  selection  of  dopants  is  available. 

4)  The  process  is  near  equilibrium  at  the  surface  thus  allowing 
reproducibility- 

5)  The  growth  occurs  below  the  film  melting  temperature. 

6)  The  growth  rate  can  be  high. 

7)  The  impurity  distribution  coefficients  are  generally  favorable. 
There  are,  however,  several  drawbacks  with  LPE.  Frequently,  the  pres¬ 
ence  of  surface  defects  such  as  incomplete  melt  removal,  terraces, 
pinholes  and  miniscus  lines  degrade  the  material.  The  thickness 
uniformity  can  be  poor  and,  for  solid  solution  films,  inherent  composi¬ 
tion  gradients  are  present.  Furthermore,  LPE  is  a  small  scale  batch 
operation  and  heteroepitaxy  can  be  difficult. 

Molecular  beam  epitaxy  is  a  method  for  growing  epitaxial  thin 
films  of  semiconductors  by  impinging  one  or  more  thermal  energy  beams 
of  atoms  or  molecules  onto  a  heated  substrate  under  ultra-high  vacuum 
conditions.  The  distinguishing  characteristic  of  MBE  is  the  slow 
growth  rate  (0.1-2  pm/hr)  that  permits  precise  control  of  layer 
thickness,  composition  and  doping  profiles.  Furthermore,  it  is  possi¬ 
ble  to  achieve  spatial  resolution  not  offered  by  other  techniques.  As 
with  LPE,  the  growth  temperatures  are  lower  than  those  encountered  in 
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CVO  systems-  Additionally,  in  situ  analysis  of  the  surface  structure 
and  reaction  conditions  during  growth  are  possible.  However,  the  equip¬ 
ment  is  very  expensive  and  the  throughput  is  low. 

Conmercially,  the  most  successful  technique  for  the  growth  of 
epitaxial  films  has  been  chemical  vapor  deposition.  The  major  advan¬ 
tages  given  by  CVO  over  MBE  and  LPE  are  relatively  fast  growth  rates 
(0.1-1  um/min.),  the  ability  to  scale  up  research  equipment  to  accommo¬ 
date  many  large  substrates,  the  availability  of  a  variety  of  source 
gases  and  the  ease  in  which  dopants  may  be  changed  during  the  deposition 
process. 

Three  source  chemistries  dominate  the  CVO  process  for  11 I -V 
materials;  group  III  metalorganic  (MOCVO),  group  V  hydride  and  group 
V  halide  sources.  The  MOCVO  technique  involves  an  irreversible  pyroly¬ 
sis  reaction  in  which  a  group  III  metalorganic  gaseous  species  is  fed 
to  a  cold-wall  reactor  along  with  a  group  V  species  (usually  a  hydride). 
These  species  then  contact  a  heated  substrate,  decompose  and  deposit 
an  epitaxial  layer  onto  the  substrate.  Due  to  the  nature  of  the  source 
gas,  carbon,  which  is  liberated  during  the  decomposition  of  the  group 
III  metalorganic  species,  is  incorporated  as  an  unintentional  acceptor 
in  MOCVO  grown  epitaxial  layers.  Other  unintentional  dopants  include 
Mg,  Zn  and  Si  which  are  usually  present  as  impurities  in  the  metalor¬ 
ganic  source  material.  The  cold-wall  design  employed  in  the  MOCVO 
system  results  in  large  amounts  of  deposition  on  the  reactor  wall 
and  dopants  which  have  been  deposited  on  the  wall  often  desorb  during 
later  stages  of  the  deposition  process.  This  makes  it  difficult  to 
grow  epitaxial  layers  with  abrupt  doping  profiles  using  MOCVO.  A  major 
advantage  associated  with  the  MOCVO  process  is  its  ability  to  success¬ 
fully  grow  epitaxial  layers  which  contain  aluminum. 
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The  most  successful  CVD  techniques  for  the  growth  of  III-V 
epitaxial  layers,  which  do  not  contain  aluminum,  have  been  the  hydride 
and  chloride  systems.  A  schematic  representation  of  the  chloride  CVD 
system  is  shown  in  Figure  1-2.  The  reactor  consists  of  source,  mixing 
and  deposition  zones  which  are  usually  operated  at  100  IcPa  pressure. 

Due  to  the  exothermic  nature  of  the  overall  deposition  reaction,  the 
reactor  is  hot-wall  design,  the  temperature  of  the  mixing  zone  is 
greater  than  or  equal  to  that  of  the  source  zone  and  the  deposition 
zone  temperature  is  normally  less  than  that  of  the  source  zone.  Typi¬ 
cally,  hydrogen  is  used  as  the  carrier  gas  with  the  concentration  of 
the  group  V  chloride  in  the  inlet  vapor  being  on  the  order  of  IS  by 
volume.  The  group  III  source  is  either  the  III-V  stoichiometric  com¬ 
pound  (in  order  to  avoid  an  initial  source  transient)  or  the  group  III 
liquid  metal  saturated  with  the  group  V  element  (generally  available  in 
higher  purity).  Upon  entering  the  reactor,  the  group  V  chloride  decom¬ 
poses  to  form  primarily  V2»  V^  and  HCl  vapor  species.  The  HCl  then 
reacts  with  the  group  III  source  to  form  III-Cl  and  other  high  chlorides. 
The  mixing  zone  allows  the  species  in  the  vapor  to  eouilibrate  while 
being  transported  to  the  lower  temperature  deposition  zone  where  the 
group  III  and  V  vapor  species  react  at  the  substrate  surface  to  deposit 
an  epitaxial  layer.  When  a  liquid  group  III  source  is  used,  the  ratio 
of  group  III  to  group  V  atoms  in  the  vapor  is  fixed  at  approximately  3 
since  essentially  all  of  the  chlorine  atoms  on  the  group  V  chloride 
react  to  predominantly  form  III-Cl.  The  use  of  a  III-V  stoichiometric 
compound  as  the  group  III  source  limits  the  III/V  ratio  to  a  maximum 
value  of  1  since  one  group  V  atom  is  released  from  the  solid  for  each 
group  III  atom  that  reacts  to  form  III-Cl. 
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The  hydride  CVD  process  is  shown  schematically  in  Figure  1-3. 

The  source  zone  of  the  hydride  system  consists  of  two  mass  transfer 
isolated  inlets,  one  for  the  group  V  species  and  one  for  the  group  III 
species.  The  group  V  hydride,  at  a  typical  inlet  composition  of  15, 
is  introduced  into  the  source  zone  where  it  decomposes  to  form  primarily 
V2,  and  H2.  As  in  the  chloride  process,  the  group  III  element  is 
transported  principally  as  the  mono-chloride  by  the  reaction  of  HCl 
(typically  the  HCl  inlet  concentration  is  15)  with  the  liquid  group  III 
metal.  One  major  advantage  the  hydride  system  provides  over  the  chlo¬ 
ride  is  the  ability  to  vary  the  vapor  III/V  ratio  by  adjusting  the  in¬ 
let  compositions  or  flowrates  of  VH^  and  HCl.  Typically,  the  source 
and  mixing  zones  in  the  hydride  system  are  operated  at  higher  temper¬ 
atures  than  those  of  the  chloride  system  in  order  to  increase  the  rate 
of  VHj  decomposition  (the  decomposition  kinetics  of  are  much  slower 
than  VClj).  Again,  hydrogen  is  usually  used  as  the  carrier  gas  and  the 
mixing  and  deposition  zones  provide  functions  equivalent  to  those  in 
the  chloride  system.  Indeed,  the  equilibrium  chemistry  of  the  two 
systems  are  identical  after  the  source  zone. 

Both  the  chloride  and  hydride  systems  are  hot-wall  designs  (heated 
from  the  outside  of  the  reactor  tube  by  conduction).  Therefore,  depo¬ 
sition  does  not  occur  on  the  reactor  wall  and  sharp  dopant  profiles  and 
heterostructures  may  be  grown.  Unfortunately,  the  hot-wall  design 
allows  interactions  between  the  vapor  and  the  reactor  wall  (usually 
quartz)  which  results  in  the  introduction  of  Si  and  0  into  the  vapor 
phase.  The  unintentional  incorporation  of  Si  into  hydride  and  chloride 
grown  epitaxial  layers  is  a  major  problem  in  these  processes. 
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The  source  materials  used  in  the  chloride  and  hydride  processes 
are  the  purest  of  all  the  III-V  CVO  systems.  All  of  the  systems  use 
H2  as  a  carrier  gas  which  is  usually  diffused  through  a  palladium 
barrier  in  order  to  obtain  extremely  pure  H2.  The  VCl^  used  in  the 
chloride  system  is  a  liquid  which  is  introduced  to  the  reactor  by 
bubbling  H2  through  it  at  a  controlled  temperature.  The  impurities 
present  in  VCl^  liquid  are  typically  very  small  and  the  bubbling  pro¬ 
cess  tends  to  further  purify  the  inlet  gas  over  that  of  the  liquid. 

The  HCl  used  in  the  hydride  system  is  a  source  of  many  impurities  if 
it  is  taken  from  a  high  pressure  gas  cylinder.  This  problem  can  be 
circumvented,  however,  by  using  AsCl^  as  a  source  of  chlorine  atoms, 
cracking  the  AsCl^  at  high  temperature  in  the  presence  of  H2,  and 
depositing  the  As  as  a  solid  in  a  trap  before  introduction  to  the 
hydride  CVO  reactor. 

An  understanding  of  the  chemistry  involved  in  the  chloride  and 
hydride  CVO  processes  is  essential  in  order  to  advance  these  technolo¬ 
gies.  In  this  study,  a  complex  chemical  equilibrium  analysis  is  pre¬ 
sented  which  identifies  the  principal  vapor  phase  species  which  must 
be  accounted  for  in  order  to  understand  these  CVO  processes.  The 
effects  of  reactor  temperature,  pressure  and  inlet  concentration  on 
the  equilibrium  chemistry  of  each  process  are  evaluated.  Further, 
species  which  may  be  added  to  or  removed  from  the  system  in  order  to 
suppress  the  amount  of  silicon  incorporated  into  the  epitaxial  layer 
are  identified.  The  influence  exerted  by  the  vapor  phase  on  the 
point  defect  structure  of  the  epitaxial  layer  is  discussed  relative 
to  the  unintentional  incorporation  of  silicon. 
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Ban  [4]  has  suggested  that  the  sluggish  decomposition  kinetics 
of  PHj  and  AsH^  may  prohibit  the  use  of  a  thermodynamic  analysis  of 
the  hydride  process.  Therefore,  the  decomposition  rates  of  NH^,  PH^ 
and  AsHj  were  studied  In  order  to  ascertain  the  degree  to  which  an 
analysis  based  on  the  assumption  of  chemical  equilibrium  can  be  applied 
to  the  hydride  process - 

Both  the  chloride  and  hydride  processes  were  Investigated  for  the 
deposition  of  homoepitaxial  6aAs  and  InP.  Due  to  the  application  of 
a  consistent  basis  set  of  operating  parameters,  direct  comparisons 
between  the  systems  are  made. 
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CHAPTER  TWO 

REVIEW  OF  THE  LITERATURE 

Impurities  in  GaAs  Epitaxial  Filins  Grown  by  the  Chloride  Process 
The  feasibility  of  applying  the  chloride  system  CVD  technique  to 
the  epitaxial  growth  of  high  purity  GaAs  was  first  demonstrated  by 
Knight  et  al .  [5]  and  Effer  [6].  Initially,  the  commercially  avail¬ 
able  AsClj  contained  sufficient  impurities  to  cause  significant  con¬ 
tamination  of  the  epitaxial  layers  and  therefore  the  purity  of  the 
feed  materials  was  believed  to  be  the  controlling  factor  for  this 
system  [7].  As  better  quality  AsCl^  became  available,  Cairns  and 
Fairman  [8,9]  and  OiLorenzo  et  al,  [10]  found  that  an  increase  in  the 
AsClj  mole  fraction  in  the  inlet  gas  stream  resulted  in  a  decrease  in 
unintentional  impurity  incorporation  in  the  epitaxial  layer. 

For  materials  grown  in  their  laboratory,  OiLorenzo  and  Moore  [11] 
identified  the  primary  unintentional  dopant  as  being  silicon,  through 
the  use  of  photoluminescence  spectra.  Further,  they  proposed  a  thermo¬ 
dynamic  model  for  the  generation  of  vapor  phase  chlorosilanes  as  a 
result  of  the  interaction  of  HCl  with  the  quartz  (Si02)  reactor  wall 
and  presented  an  expression  for  the  activity  of  solid  silicon  (i.e.  as 
an  impurity)  as  a  function  of  the  partial  pressures  of  the  chloro¬ 
silanes.  Their  model  showed  that  increasing  the  vapor  HCl  concentra¬ 
tion  (e.g.  as  a  result  of  AsCl^  decomposition)  decreased  the  condensed 
phase  silicon  activity  by  further  stabilizing  the  silicon  species  in 
the  vapor  phase  in  the  form  of  chlorosilanes.  Additionally,  their 
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model  predicted  that  the  generation  of  vapor  phase  silicon  species 
could  be  suppressed  by  the  introduction  of  H^O  vapor  into  the  system. 

Rai-Choudhury  [12]  performed  a  thermodynamic  analysis  on  the  in¬ 
corporation  of  silicon  into  GaAs  epitaxial  layers.  His  results  reflected 
those  of  OiLorenzo  and  Moore  [11]  when  considering  the  effects  of  H2O 
and  HCl,  but  he  also  showed  that  higher  temperatures  increased  the 
amount  of  vapor  phase  sil icon  species. 

The  work  of  Ashen  et  al,  [13]  further  supported  the  conclusion  that 
silicon  was  an  impurity  in  GaAs.  A  8N  lined  reactor  was  used  to  grow 
epitaxial  layers  from  liquid  Ga  sources  which  were  doped  with  Si.  Com¬ 
paring  the  electrical  characteristics  of  these  epitaxial  layers  to  layers 
grown  from  pure  Ga  sources  provided  qualitative  evidence  for  the  presence 
of  Si  in  GaAs.  The  effect  of  AsCl^  concentration  on  the  amount  of  Si  in¬ 
corporated  into  the  epitaxial  layer  was  also  verified  by  their  experiments. 
Additionally,  these  studies  provided  evidence  which  indicated  that  Si  was 
much  more  likely  to  reside  on  a  Ga  site  than  an  As  site  and  therefore, 
behaves  as  a  donor.  This  conclusion  was  also  supported  by  Beiden  et  al . 
[14]. 

Wolfe,  Stillman  and  Korn  [15]  have  identified,  through  intentional 
doping  and  determination  of  ionization  energies,  three  unintentional 
impurities.  Si,  C  and  one  unknown  (possibly  Te),  in  GaAs  grown  by  the 
chloride  CVD  system.  Also,  due  to  the  results  of  Solomon  [16]  which 
showed  that  oxygen  may  be  a  shallow  donor  in  GaAs,  they  attempted  to  dope 
the  epitaxial  layer  with  oxygen  by  adding  6020^  to  the  liquid  gallium 
source.  The  oxygen,  however,  was  not  incorporated  into  the  epitaxial 
layer,  but  did  reduce  the  amount  of  silicon  which  was  incorporated  into 
the  layer.  This  reduction  in  background  doping  due  to  the  presence 


of  oxygen  was  also  investigated  by  Palm  et  al.  [17]  by  injecting  a 
hydrogen-oxygen  mixture  into  the  mixing  zone  of  a  chloride  system 
CVD  reactor.  Using  silane  as  an  intentional  dopant,  the  presence  of 
oxygen  was  found  to  reduce  the  incorporation  of  silicon  in  the  epitaxial 
layers  by  as  much  as  four  orders  of  magnitude. 

Seki  et  al.  [18]  performed  a  thermodynamic  analysis  of  the  GaAs 
chloride  process  in  order  to  identify  the  effects  of  inerts,  HCl  and 
substrate  temperature  on  the  activity  of  silicon  in  the  epitaxial 
layers.  The  analysis  predicted  that  increasing  the  HCl  content  or  de¬ 
creasing  the  substrate  temperature  lowered  the  silicon  activity.  In 
addition,  replacing  the  hydrogen  carrier  gas  with  an  inert  gas  caused 
a  very  large  reduction  in  the  silicon  activity. 

The  effect  of  replacing  the  hydrogen  carrier  gas  with  an  inert  was 
investigated  by  Ozeki  et  al .  [19].  Through  far  infrared  photocon¬ 
ductivity  measurements,  it  was  determined  that  sulfur  was  the  dominant 
residual  donor  present  in  epitaxial  GaAs  when  Ng  was  used  in  place  of 
H2  as  the  carrier  gas.  It  was  also  found  that  the  dominant  residual 
donor  when  H2  was  used  as  the  carrier  gas  was  sometimes  Si  and  some¬ 
times  S.  A  possible  source  of  S  in  the  system  was  not  discussed 
(although  it  was  presumably  in  the  feed  gases}  and  elaboration  on  the 
growth  conditions  which  caused  Si  or  S  to  be  dominant  was  not  provided. 

A  thermodynamic  analysis  of  the  chloride  CVO  system  performed  by 
Boucher  and  Hoi  Ian  [20]  assumed  solid  GaAs  as  the  group  III  source 
material.  By  comparison  with  experiment,  it  was  found  that  the 
dominant  group  III  and  group  V  species  present  in  the  vapor  were  GaCl 
and  As^.  Under  the  experimental  conditions  investigated,  the  deposition 
process  was  kinetically  controlled  with  an  activation  energy  of  -40 
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kcal/mole,  and  reproducible  growth  conditions  could  be  attained  only 
if  105  or  less  of  the  thermodynamically  available  GaAs  was  deposited 
from  the  vapor  phase . 

Centner  et  al .  [21]  also  studied  the  chloride  process  experimentally 
and  presented  a  thermodynamic  analysis  over  a  greater  range  of  tempera¬ 
ture,  pressure  and  inlet  ASCI3  composition  than  did  previous  investiga¬ 
tors.  They  found  that  AS2  became  the  dominant  group  V  species  below 
10  kPa  pressure  and  that  GaCl  was  always  the  dominant  group  III  species. 
At  large  AsCl^  inlet  compositions,  the  higher  gallium  chlorides  became 
more  pronounced  but  never  competed  with  the  monochloride  as  the  dominant 
species.  They  concluded,  based  on  a  kinetic  model  [22],  a  mass  transfer 
model  and  experimental  results,  that  the  deposition  of  GaAs  was  kinet- 
ically  rather  than  mass  transfer  controlled. 

Shaw  [23]  studied  the  transport  kinetics  of  the  GaAs  chloride 
system  in  the  source  and  deposition  zones.  He  found  an  activation 
energy  of  49.1  kcal/mole  in  reasonable  agreement  with  that  of  Boucher 
and  Hollan  [20],  for  a  surface  reaction  associated  with  the  deposition 
process. 

Impurities  in  InP  Epitaxial  Films  Grown  by  the  Chloride  Process 

The  epitaxial  growth  of  InP  using  a  chloride  CVD  system  was  first 
demonstrated  by  Clark  et  al .  [24]  and  later  by  Hales  et  al.  [25].  Both 
groups  of  investigators  reported  limitations  on  the  purity  of  their 
epitaxial  layers  due  to  unintentional  dopants.  Joyce  and  Williams  [26] 
tentatively  identified  the  impurities  as  being  Si  and  Zn  acceptors. 

They  also  found  evidence  of  a  donor  which  was  thought  to  be  amphoteric 
Si. 


The  dependence  of  unintentional  doping  on  PCl^  mole  fraction  in 
the  InP  chloride  system  was  first  reported  by  Clark  [27].  The  simi¬ 
larity  between  the  6aAs  and  InP  chloride  system  reactors  combined  with 
the  analogous  dependencies  on  the  group  V  hydride  mole  fraction  supported 
the  belief  that  Si  was  an  impurity  in  InP  epitaxial  layers.  Clarke  [28] 
later  studied  the  effect  of  III/V  ratio  in  the  vapor  phase  on  the  un¬ 
intentional  doping  of  InP  epi -layers  and  found  p-type  conductivity  for 
III/V  <  3  and  n-type  for  III/V  >  3,  with  a  minimum  in  the  free  carrier 
concentration  at  III/V  «  3.  No  explanation  was  offered  for  these 
observations. 

Easton  [29]  investigated  the  unintentional  doping  of  InP  epitaxial 
layers  grown  by  the  chloride  system  and  concluded  that  S  (acting  as  a 
donor)  was  the  major  impurity  and  that  the  origin  of  the  S  was  the  PCI2 
liquid  source.  Using  mass  spectrometric  analysis,  Easton  found  Si,  S 
and  Zn  present  in  the  PCl^  source  at  levels  between  1  ppm  and  10  ppm 
(atomic)  and  Fe,  Cu,  Cd  and  Sn  at  '0.7  ppm.  These  same  elements  were 
found  in  the  unused  bulk  In  liquid  at  levels  below  0.1  ppm.  Analysis 
of  the  used  In  source  liquid  revealed  impurity  levels  approximately  10 
times  larger  than  those  in  the  unused  liquid. 

These  results  support  the  work  of  Weiner  [30]  who  proposed  models 
for  the  contamination  of  a  Ga  liquid  source  by  Si  in  the  GaAs  and  GaP 
systems.  Weiner's  results  showed  that  the  liquid  group  III  metal  im¬ 
purity  level  increased  as  the  exposure  to  the  CVD  environment  increased. 
He  also  found  the  level  of  Si  contamination  to  be  inversely  related  to 
the  partial  pressure  of  H2O  in  the  system. 

Fairhurst  et  al .  [31]  studied  the  InP  halide  system  using  both 
PCI 2  and  P8r2.  They  found  that  oxyhalide  salts  were  present  in  both 
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phosphorous  liquids  at  approximately  the  100  ppm  level.  The  presence 
of  oxygen  was  expected  to  decrease  the  level  of  impurity  incorporation 
in  the  epitaxial  layers.  This  effect  was  not  observed  however,  pre¬ 
sumably  due  to  this  level  of  oxygen  contamination  being  too  low  to  be 
significant.  Equilibrium  calculations  were  performed  which  showed 
InCl  and  to  be  the  dominant  group  III  and  V  species  in  the  vapor 
over  a  temperature  range  of  850  K  to  1150  K  and  an  inlet  PCl^  mole 
fraction  range  of  O.IS  to  6*.  These  results  agreed  with  those  of 
Boucher  and  Hoi  Ian  [20]  for  the  analogous  GaAs  system. 

Hales  and  Knight  [32]  investigated  the  effect  of  introducing 
oxygen  into  the  system  in  order  to  reduce  the  level  of  impurities  in 
InP.  They  observed  a  monotonic  decrease  in  free  electron  density  for 
additions  of  O2  up  to  24  ppm.  The  electron  mobility  (measured  at 
77  K)  however  reached  a  very  broad  maximum  at  approximately  16  ppm 
of  added  O2,  which  suggests  that  oxygen  was  becoming  incorporated 
into  the  epitaxial  layer  and  that  there  is  a  limit  to  the  degree  of 
benefit  which  may  be  obtained  through  oxygen  addition.  They  also 
observed  POCl^  to  be  an  impurity  in  the  liquid  PCl^  used  in  the 
chloride  system. 

Investigations  of  the  dependence  of  impurity  incorporation  on 
PClj  inlet  composition,  total  flowrate  and  deposition  zone  temperature 
were  carried  out  by  Chevrier  et  al .  [33].  They  observed  a  decrease 
in  free  carrier  concentration  with  increasing  PCl-j  concentration, 
as  did  other  investigators,  but  also  found  that  the  impurity  concen¬ 
tration  increased  with  increasing  total  flowrate.  This  velocity  effect 
had  not  been  reported  before  and  suggests  the  presence  of  a  mass 
transfer  resistance  at  the  group  III  source  (if  impurities  are  picked 
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up  from  the  liquid  metal)  or  at  the  substrate  in  the  deposition  zone. 
They  also  studied  the  intentional  doping  of  InP  as  a  function  of 
deposition  zone  temperature  using  SnCl^.  Low/er  free  electron  con¬ 
centrations  and  higher  electron  mobilities  were  observed  when  the 
deposition  zone  temperature  was  decreased  from  950  K  to  910  K.  Thus, 
the  uptake  of  group  IV  impurities  (Sn,  Si,  etc.)  was  apparently  re¬ 
duced  by  lowering  the  deposition  zone  temperature. 

Cardwell  et  al .  [34]  found  kinetic  limitations  in  both  the  source 
and  deposition  zones.  The  previously  reported  effect  of  PCl^  mole 
fraction  on  impurity  levels  in  the  epitaxial  layers  was  observed. 
Intentional  doping  of  InP  using  Sn  followed  the  same  behavior  as  that 
of  unintentional  dopants  and  therefore  supports  the  use  of  Sn  for 
studies  regarding  the  reduction  of  unintentional  impurities.  In  con¬ 
trast  to  Chevrier  et  al .  [33],  no  dependence  of  impurity  uptake  on 
total  flowrate  was  found. 

A  thermodynamic  analysis  comparing  the  GaAs  and  InP  chloride 
systems  using  the  stoichiometric  III-V  solid  as  the  group  III  source 
material  was  reported  by  Shaw  [35].  His  results  also  confirmed  GaCl, 
As^,  InCl  and  P^  to  be  the  dominant  group  III  and  V  vapor  species  in 
these  systems.  Further,  the  degree  of  supersaturation  in  the  deposi¬ 
tion  zone  was  calculated  to  be  less  for  InP  than  for  GaAs  under 
analogous  conditions.  Since  solid  III-V  source  materials  were  em¬ 
ployed  etching  conditions  were  predicted  whenever  the  deposition  zone 
temperature  was  greater  than  that  of  the  source  zone. 

Impurities  in  GaAs  Epitaxial  Films  Grown  by  the  Hydride  Process 

The  feasibility  of  applying  a  hydride  CVD  system  for  the  growth  of 
high  purity  GaAs  was  demonstrated  by  Enstrom  and  Peterson  [36].  Since 
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the  hydride  system  consists  of  a  hot-wall  quartz  reactor  with  the 
elements  H,  Cl,  Ga  and  As  present  in  the  vapor,  one  would  expect  it 
to  show  an  unintentional  impurity  incorporation  problem  similar  to  that 
of  the  chloride  system.  Pogge  and  Kemlage  C37]  investigated  the 
effects  of  HCl ,  AsH^  and  PH^  on  the  unintentional  doping  of  GaAs 
and  GaP  grown  by  the  hydride  system.  They  found  that  the  free  carrier 
concentration  decreased  with  increasing  HCl,  AsH^  or  PH^  composition. 

The  effect  of  HCl  was  less  than  that  of  the  group  V  hydrides  and 
changes  in  PH^  showed  larger  effects  than  did  AsH^.  They  concluded 
that  the  HCl  effect  on  the  vapor  phase  composition  was  similar  to  that 
of  the  chloride  system.  Further,  they  concluded  that  As^  and  P^  caused 
blockage  of  the  available  surface  sites  on  the  substrate  due  to  the 
large  size  of  these  molecules.  The  unintentional  dopant  was  assumed 
to  be  Si  generated  from  reactions  with  the  quartz  wall. 

Kennedy  et  al .  [38]  investigated  the  effect  of  HCl  inlet  com¬ 
position  and  additions  of  HCl  downstream  of  the  source  zone  on  the  un¬ 
intentional  doping  of  GaAs  grown  in  a  hydride  CVD  reactor.  Increasing 
the  HCl  inlet  composition  greatly  reduced  the  free  carrier  density  in 
the  epitaxial  layer.  In  contrast  to  this  result,  however,  when  HCl 
was  added  downstream  of  the  source  zone  the  free  carrier  density  was 
found  to  increase.  This  led  to  the  conclusion  that  the  equilibrium 
model  as  proposed  by  OiLorenzo  and  Moore  [11]  for  the  chloride  system 
was  not  applicable  to  the  hydride  system.  However,  the  HCl  which  was 
injected  may  not  have  been  as  pure  as  that  which  was  generated  from 
the  decomposition  of  AsCl^  in  the  chloride  system  and  therefore  may 
have  introduced  additional  impurities  into  the  epi-layer.  These  effects 
were  also  observed  by  Enstrom  and  Appert  [39]. 
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The  work  of  Skromme  et  al.  [40]  identified  some  of  the  unintentional 
donors  and  acceptors  present  in  6aAs  and  InP  prepared  by  the  hydride 
CVD  system.  They  found  C,  Zn,  Cu  and  Mn  as  acceptors  and  Si,  S  and  Ge 
as  donors  in  GaAs.  Epitaxial  InP  was  found  to  contain  Zn,  C  or  Mg  and 
an  unidentifiable  acceptor  along  with  Si  and  S  as  donors.  Additionally, 
in  one  of  the  laboratories  (Honeywell)  where  the  epi-layers  were  grown, 
an  increase  in  the  impurity  concentration  in  epitaxial  GaAs  was  noted  as 
the  HCl  gas  cylinders  "aged".  This  effect,  however,  was  not  observed 
at  the  other  laboratory  (Hanscom  AFB). 

The  effect  of  pressure  was  studied  experimentally  by  Putz  et  al . 

[41]  from  1  kPa  to  100  kPa.  They  found  that  the  unintentional  doping 
of  GaAs  was  reduced  at  pressures  below  IQO  kPa. 

Impurities  in  InP  Epitaxial  Films  Grown  by  the  Hydride  Process 

Growth  of  InP  epitaxial  layers  using  the  hydride  system  has  been 
demonstrated  by  Olsen  [42]  and  Hyder  [43]  among  others.  Both  of  these 
investigators  observed  unintentional  impurity  incorporation  similar  to 
that  occurring  in  the  GaAs  system.  Hyder  also  found  that  for  the  ter¬ 
nary  In^Gai_^As  (x=0.53),  a  maximum  in  electron  mobility  occurred  when 
the  III/V  ratio  in  the  vapor  phase  was  approximately  2,  but  the  effect 
of  III/V  ratio  on  free  carrier  concentration  was  not  discussed. 

Zinkiewicz  et  al.  [44]  also  studied  the  growth  of  InP  and  the  ternary 
In^Gai_^As  in  the  hydride  system.  They  found  Zn,  Cu  and  Hg  to  be  pre¬ 
sent  as  unintentional  donors. 

Anderson  [45]  studied  the  hydride  system  for  InP  growth  in  order 
to  determine  the  effect  of  HCl  mole  fraction,  Hg  flowrate  and  mixing 
zone  temperature  on  unintentional  impurity  incorporation.  He  found  that 
these  parameters  caused  only  minor  changes  in  the  electrical  behavior 
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of  the  InP  epitaxial  layers.  This  suggests  that  the  InP  hydride  sys¬ 
tem  may  perform  somewhat  differently  than  the  GaAs  hydride  system. 

Jones  [46]  performed  a  thermodynamic  analysis  of  the  InP  hydride 
system  in  order  to  understand  the  effect  of  process  parameters  on  unin¬ 
tentional  Si  incorporation.  The  calculations  predicted  that  decreasing 
temperatures  lowered  the  silicon  activity  in  the  epitaxial  layer.  Addi¬ 
tionally,  the  silicon  activity  was  decreased  by  increasing  the  PH^  and/or 
HCl  inlet  composition.  Very  little  effect  was  noted  when  HCl  was  added 
downstream  of  the  source  zone.  His  analysis  used  liquid  In  as  the  group 
III  source  material  and  showed  InCl  and  P^  to  be  the  dominant  group  III 
and  V  vapor  species. 

Ban  and  Ettenberg  [47]  coupled  a  mass  spectrometer  to  a  hydride 
system  reactor  used  for  the  growth  of  Inj^Ga^_j^P.  They  compared  measured 
vapor  species  to  those  predicted  by  a  thermodynamic  model  and  concluded 
that  the  model  yielded  an  acceptable  representation  of  the  system.  The 
major  shortcomings  of  the  model  were  an  over predict ion  of  the  amount  of 
InCl  generated  from  the  heterogeneous  reaction  of  HCl  and  In  liquid,  and 
a  predicted  greater  degree  of  dissociation  for  PH^  than  was  measured. 

Due  to  the  slow  decomposition  kinetics  of  PH^  and  the  potential  mass 
transfer  and  kinetic  limitations  associated  with  heterogeneous  reactions, 
these  discrepancies  were  not  surprising.  Their  mass  spectrometric 
measurements  identified  the  major  vapor  phase  species  as  being  InCl, 

GaCl ,  HCl ,  PH3,  P2.  P4  and  H2. 

Usui  and  Watanabe  [48]  investigated  the  effects  of  temperature  and 
oxygen  additions  on  the  unintentional  doping  of  hydride  grown  InP.  Ad¬ 
ditions  of  O2  in  the  1  ppm  to  10  ppm  range  decreased  the  free  carrier 
concentration  about  one  order  of  magnitude,  but  further  additions  caused 


the  free  earner  concentration  to  increase,  presumably  due  to  uptake 
of  oxygen  by  the  epitaxial  layer.  The  liquid  In  source  that  was  used 
in  these  experiments  was  found  to  have  a  gettering  effect  on  impuri¬ 
ties  in  the  inlet  gases.  Unused  In  showed  less  than  1  ppm  levels  of 
Si,  S,  Sn,  Te,  Zn,  Fe  and  Cu,  while  In  exposed  to  the  CVD  environment 
contained  increased  levels  (-2  ppm)  of  Fe,  Cu  and  Sn.  Increasing  the 
source  zone  temperature  caused  a  decrease  in  the  free  carrier  concen¬ 
trations  in  InP  epi-layers  due  to  an  increased  ability  of  the  In 
liquid  to  getter  impurities  at  high  temperature.  Thus,  the  purity  of 
source  gases  still  appears  to  be  a  major  problem  in  the  hydride  system. 

The  Thermal  Decomposition  of  NH^ 

The  thermal  decomposition  of  the  trihydrides  of  N,  P  and  As  have 
been  studied  by  many  investigators  and,  for  temperatures  below  1500  K, 
a  general  consensus  exists  that  these  reactions  are  almost  entirely 
heterogeneous  in  nature.  Bamford  and  Tipper  [49]  have  reviewed  the 
literature  relevant  to  the  homogeneous  pyrolysis  of  ammonia  at  temper¬ 
ature  above  2000  K  and  found  the  reaction  to  be  characterized  by  an 
activation  energy  of  approximately  100  kcal/mole.  Based  on  the  ob¬ 
served  activation  energy  and  the  results  of  experiments  with  deuterated 
ammonia,  the  initiating  step  in  the  pyrolysis  reaction  sequence  was 
proposed  to  be: 

NH3  +  M  — >  NH  +  H2  +  M  2-1 

where  M  represents  any  gas  molecule.  They  also  found  evidence  that  a 
reaction  which  forms  NH^  is  likely  to  be  present  in  the  decomposition 
chain  reaction  sequence,  but  were  unable  to  identify  the  nature  of 
this  reaction. 
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The  decomposition  of  NH^  in  a  quartz  vessel  was  first  studied  by 
Bodenstein  and  Kranendieck  [50]  using  a  manometric  method.  The  amount 
of  surface  area  present  in  the  reactor  was  varied  by  the  addition  of 
quartz  fibers.  They  concluded  that  within  the  temperature  range  of 
their  study  (1063  K  to  1153  K) ,  the  reaction  appeared  to  be  first  order 
and  was  entirely  heterogeneous  in  nature.  Further,  they  found  no 
change  in  the  reaction  rate  when  or  N2  additions  were  made  to  the 
system. 

Hinshelwood  and  Burke  [51]  investigated  NH^  decomposition  in  a 
quartz  vessel  for  temperatures  as  high  as  1323  K.  They  also  concluded 
that  the  reaction  was  dominated  by  the  heterogeneous  component  and 
demonstrated  a  first  order  dependence  on  NH^.  Additions  of  H2  to  the 
reactor  decreased  the  reaction  rate  while  N2  additions  were  ineffective. 

Christiansen  and  Knuth  [52]  suggested  the  following  mechanism 
for  the  heterogeneous  pyrolysis  of  NH^  in  a  quartz  vessel: 


+  S  NH  +  H2  +  S 

2-2 

NH  NH* 

2-3 

< — 

+  NH3  N2  +  2H2 

2-4 

where  S  represents  a  surface  site.  Their  experiments  were  carried  out 

in  a  reactor  vessel  having  a  surface  to  volume  ratio  (S/V)  of  1  and  a 

2 

surface  area  of  0.02  H  .  They  concluded  that  the  forward  component  of 
reaction  2-2  was  the  rate  limiting  step  and  over  the  temperature  range 
of  1062  K  to  1132  K,  this  reaction  was  characterized  by  an  activation 
energy  of  43  +  5  kcal/mole  with  an  Arrhenius  type  frequency  factor  of 
4.5  X  10^  s‘^ . 

The  investigation  of  Russow  and  Pewsner  [53]  into  the  decomposi¬ 
tion  of  in  a  quartz  reactor  demonstrated  that  the  reaction  followed 


a  first  order  dependence  with  respect  to  NH^  partial  pressure.  They 
reported  an  activation  energy  of  38.2  kcal/mole  for  the  pyrolysis 
reaction. 

The  decomposition  of  NH^  on  quartz  was  reported  by  Szabo  and 
Ordogh  [54]  to  follow  a  1/2  order  dependence  with  NH^  pressure  at  913  K 
and  a  first  order  dependence  at  1013  K.  The  reduced  order  of  reaction 

at  913  K  was  reported  to  be  a  result  of  H2  competing  with  NH^  for  ad¬ 

sorption  sites  on  the  quartz  reactor  wall.  The  presence  of  or  O2 
in  the  system  was  found  to  decrease  the  reaction  rate  while  no  reaction 

rate  changes  were  observed  when  N2  was  added.  Initial  partial  pressures 

2  4 

of  NH^  ranging  from  6.6  x  10  Pa  to  2.6  x  10  Pa  were  tried,  but  no 
changes  were  observed  in  the  order  of  reaction  over  this  range. 

An  activation  energy  of  34  +  2  kcal/mole  for  the  decomposition  of 
NHj  on  quartz  was  reported  by  Voelter  and  Schoen  [55].  The  frequency 
factor  associated  with  an  Arrhenius  type  temperature. dependence  was 
560  s"\  but  they  did  not  report  the  surface  area  of  the  reactor  used. 

Over  the  temperature  range  of  1023  K  to  1173  K,  the  reaction  was  found 

to  be  first  order  with  respect  to  pressure. 

Kelvin  [56]  utilized  an  infrared  spectrometer  to  measure  the  out¬ 
let  NHj  composition  from  a  plug  flow  quartz  reactor  over  the  tempera¬ 
ture  range  833  K  to  1373  K.  He  found  that  the  reaction  rate  varied 
with  the  reactor  S/V  to  the  0.75  power,  confirming  the  heterogeneous 
nature  of  the  reaction.  Hydrogen  was  found  to  exhibit  a  strong  in¬ 
hibitory  effect  on  the  decomposition  reaction  nature  due,  primarily, 
to  a  reaction  between  NH  radicals  and  H2,  which  forms  NH^.  Additions 
of  N2,  Ar  or  He  to  the  reactor  resulted  only  in  a  dilutent  effect. 
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The  Thermal  Decomposition  of  PH-j 

The  decomposition  of  PH^  in  a  quartz  reactor  was  first  studied  by 
van't  Hoff  and  Kooj  [57].  They  found  the  reaction  to  be  first  order 
in  PHj  pressure  over  the  temperature  range  310  K  to  512  K.  The  reac¬ 
tion  was  believed  to  be  heterogeneous  as  a  result  of  the  increase  in 
reaction  rate,  which  occurred  upon  the  addition  of  quartz  fibers  to  the 
system.  An  activation  energy  of  46.4  kcal/mole  was  reported. 

Trautz  and  Bhandarkar  [58]  reported  a  transition  from  hetero¬ 
geneous  to  homogeneous  reaction  kinetics  at  940  K  for  the  decomposition 
of  PH^  in  a  porcelain  bulb.  They  reported  activation  energies  of  59 
kcal/mole  for  temperatures  below  940  K  and  116  kcal/mole  for  tempera¬ 
tures  above  940  K.  Based  on  the  large  degree  of  scatter  in  their  re¬ 
sults  and  the  fact  that  other  investigators  have  not  seen  this  transi¬ 
tion,  it  is  doubtful  that  a  homogeneous  decomposition  reaction  was 
actually  observed. 

Hinshelwood  and  Topley  [59]  studied  the  decomposition  of  PH^  in  a 
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155  cm  quartz  bulb  with  surface  areas  of  210  cm  to  1600  cm  .  They 

concluded  that  in  the  temperature  range  848  K  to  1042  K,  the  reaction 
was  first  order  and  behaved  in  a  heterogeneous  manner  with  an  acti¬ 
vation  energy  of  46  +  4  kcal/mole.  The  reaction  rate  was  found  to  in¬ 
crease  with  increasing  S/V  to  the  0.8  power. 

Oevyatykh  et  al .  [60]  studied  the  decomposition  of  PH^  on  glass 

and  Si  over  the  temperature  range  740  K  <  T  <  822  K.  They  found  PH^ 
decomposition  to  be  first  order  with  activation  energies  of  44.2  kcal/ 
mole  and  55.3  kcal/mole  on  the  glass  and  Si  surfaces,  respectively. 

The  decomposition  of  SbH^  was  studied  on  an  antimony  surface  and  was 
found  to  have  an  activation  energy  of  7.7  kcal/mole.  The  SbH^  decom¬ 
position  reaction  was  investigated  over  the  temperature  range  293  K  < 
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T  <  364  K  and  the  reaction  order  was  observed  to  change  from  half  or¬ 
der  at  298  K  to  first  order  at  364  K. 

The  Thermal  Decomposition  of  AsH^ 

The  decomposition  of  AsH^  on  glass.  As  and  Sb  has  been  studied  by 
Tamaru  [61].  The  reaction  was  found  to  proceed  most  rapidly  on  the  Sb 
surface  and  slowest  on  glass.  Adding  H2  to  the  system  had  no  effect 
on  reaction  rate  and  no  isotopic  exchange  was  observed  when  D2  was 
added.  Tamaru  proposed  that  the  reaction  mechanism  consisted  of  AsH^ 
adsorbing  on  the  surface  followed  by  sequential  stripping  of  the  hydro¬ 
gen  atoms  off  of  the  As  atom.  He  believed  the  rate  determining  step 
to  be  the  removal  of  the  first  H  atom  and  assigned  an  activation  energy 
of  23.2  kcal/mole  to  this  reaction.  Tamaru  [62]  later  attempted  the 
calculation  of  the  reaction  rate  constants  for  AsH^  and  SbH^  decomposi¬ 
tion  on  As  and  Sb  surfaces  using  a  model  based  on  activated  complex 
theory.  His  predicted  rate  constant  for  AsH^  decomposition  was  six 
orders  of  magnitude  below  the  observed  value  while  the  predicted  rate 
constant  for  SbH^  was  two  orders  of  magnitude  low.  The  restrictive 
assumptions,  which  required  all  of  the  hydrogen  bond  energies  to  be  the 
same  and  all  partition  functions  to  have  the  value  of  one,  were  proba¬ 
bly  the  reasons  for  the  poor  results. 

Devyatykh  et  al .  [63]  found  that  the  decomposition  of  AsH^  on  a 
Si  surface  obeyed  first  order  kinetics  and  was  characterized  by  an 
activation  energy  of  50.9  kcal/mole.  Their  experiments  were  conducted 
from  a  temperature  of  659  K  to  707  K. 

Kedyarkin  and  Zorin  [64]  reported  an  activation  energy  of  25.5 
kcal/mole,  in  good  agreement  with  that  of  Tamaru  [61],  for  the  decom¬ 
position  of  AsH^  on  As.  They  found  the  reaction  to  be  first  order  with 


respect  to  AsH^  pressure  and  heterogeneous  in  nature.  The  tempera¬ 
ture  range  of  their  study,  however,  was  somewhat  restricted  (543  <  T  < 
583  K). 

The  only  reported  study  of  AsH^  decomposition  on  a  quartz  sur¬ 
face  was  undertaken  by  Frolov  et  al.  [65].  Their  investigations  were 
carried  out  on  quartz,  Ge,  Ga  and  GaAs  surfaces  over  a  temperature 
range  of  699  K  to  909  K.  The  activation  energy  for  AsK^  decomposition 
on  quartz  was  found  to  be  32.6  kcal/mole  with  an  Arrhenius  frequency 
factor  of  513  s”^ .  The  activation  energies  for  AsH^  decomposition  on 
Ge,  Ga,  Te  doped  GaAs  and  Cr  doped  GaAs  were  reported  as  54,  30,  45 
and  27  kcal/mole,  respectively,  in  the  presence  of  H2.  When  He  was 
used  as  the  carrier  gas  in  place  of  H2,  these  activation  energies  in¬ 
creased  slightly.  No  explanation  was  provided  for  this  observation. 
Their  experimental  apparatus  consisted  of  an  open  tube  coupled  to  an 
infrared  spectrometer.  The  surface  area  or  surface  to  volume  ratio  of 
the  reactor  was  not  reported. 


CHAPTER  THREE 

FORMULATION  OF  AN  EQUILIBRIUM  MODEL 
FOR  THE  HYDRIDE  AND  CHLORIDE  PROCESSES 

Introduction 

Performing  a  detailed  analysis  on  CVO  processes  as  complex  as  the 
hydride  and  chloride  processes  presents  a  formidable  task  if  results 
which  allow  direct  comparison  of  the  systems  are  desired.  A  rigorous 
treatment  would  require  the  solution  of  the  mass,  energy  and  momentum 
equations  with  variable  properties  and  kinetic  expressions  for  many  of 
the  homogeneous  and  heterogeneous^ reactions  which  are  present.  The 
current  lack  of  knowledge  surrounding  the  fundamental  chemistry  which 
underlies  these  processes  precludes  this  type  of  analysis.  Many  of 
the  essential  differences  between  these  processes  may,  however,  be 
elucidated  by  the  application  of  a  model  based  on  the  assumption  of 
chemical  equilibrium.  Due  to  the  high  temperatures  employed  in  these 
processes,  it  is  expected  that  many  of  the  chemical  reactions  proceed 
at  a  very  high  rate.  Therefore,  homogeneous  reactions  are  expected  to 
be  near  equilibrium.  Heterogeneous  reactions  may  be  impeded  by  the 
existence  of  mass  transfer  barriers.  These  barriers  may  be  approxi¬ 
mated  by  either  allowing  the  reaction  to  reach  equilibrium  or  by  ne¬ 
glecting  the  reaction.  The  results  related  to  heterogeneous  reactions 
are  therefore  expected  to  be  semiquantitative  in  nature  in  that  they 
provide  limits  to  the  system  composition.  These  results  may  then  be 
compared  to  experiments  in  order  to  determine  the  degree  of  influence 
exerted  by  the  heterogeneous  reactions  on  the  systems. 
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The  resources  employed  in  order  to  affect  the  analysis  per¬ 
formed  in  this  study  include  a  computer  code  for  the  calculation 
of  multiphase  equilibria  in  systems  having  many  chemical  species 
(typically  more  than  20),  models  for  each  zone  of  these  CVO  processes 
which  are  based  on  chemical  equilibrium  constrained  to  account  for 
actual  mass  transfer  or  kinetic  limitations,  a  model  for  the  point 
defect  chemistry  in  the  solid  epitaxial  layer  and,  finally,  a  set 
of  consistent  thermodynamic  properties  for  each  chemical  species  in 
the  system.  The  computer  code,  reactor  zone  models  and  solid  state 
defect  chemistry  are  discussed  in  this  Chapter.  The  choice  and 
analysis  of  a  set  of  thermodynamic  properties  are  described  in  Chap¬ 
ter  Four. 

Method  of  Calculation  for  Complex  Chemical  Equilibrium 

The  calculation  of  complex  chemical  equilibrium  in  nulticomponent, 
multiphase  systems  has  been  reviewed  most  completely  by  Smith  [66]. 
Essentially,  there  are  two  statements  of  the  solution  to  this  problem. 
Nonstoichiometric  methods,  such  as  the  popular  Rand  algorithm  [67], 
directly  minimize  the  Gibbs  energy  of  the  total  system  in  order  to  ob¬ 
tain  a  solution  without  recourse  to  a  specific  set  of  formation  re¬ 
action  equations.  Stoichiometric  methods  [68]  require  that  an  inde¬ 
pendent  set  of  chemical  reactions  be  in  equilibrium.  Generally,  a 
formation  reaction  is  written  for  each  species  present  in  the  system 
and  the  corresponding  equilibrium  constant  for  each  reaction  is  gen¬ 
erated  from  the  Gibbs  energy  change  of  the  reaction. 

An  extension  of  the  Rand  algorithm  to  include  not  only  a  gas 
phase  with  an  inert  speces,  but  also  a  multicomponent  solution  and 
pure  condensed  phases,  was  developed  by  Anderson  [69].  This  algorithm 


was  initially  applied  to  the  hydride  and  chloride  CVD  systems,  but 
was  susceptible  to  finding  local  minima.  In  particular,  some  com¬ 
ponent  mole  fractions  sought  were  as  low  as  0.1  ppb.  The  contribu¬ 
tion  to  the  system  Gibbs  energy  for  these  components  was  quite  small 
and  the  resulting  component  mole  changes  were  not  capable  of  releasing 
the  Gibbs  energy  of  the  system  from  the  local  minima. 

A  stoichiometric  algorithm,  presented  in  Appendix  A,  was  there¬ 
fore  developed  which  was  superior  to  the  extended  Rand  algorithm  since 
a  linearized  Gibbs  energy  function  was  not  required.  The  stoichio¬ 
metric  algorithm  performed  well  for  all  of  the  systems  studied  (the 
results  were  compared  to  other  investigator's  calculations  and  were 
found  to  be  independent  with  respect  to  initial  guesses)  and  yielded 
results  which  were  in  agreement  with  those  of  the  extended  Rand  al¬ 
gorithm,  when  it  could  be  successfully  applied.  The  amount  of  com¬ 
puter  memory  required  for  the  stoichiometric  algorithm  was  found  to 
be  much  less  than  that  required  by  the  extended  Rand  algorithm  in 
order  to  solve  identical  systems. 

The  input  data  which  was  required  in  order  to  perform  the  cal¬ 
culations  consisted  of  the  standard  enthalpy  and  entropy  of  formation 
and  heat  capacity  for  each  species  along  with  the  system  temperature, 
pressure  and  inlet  composition. 

Aside  from  numerical  difficulties,  the  two  major  factors,  which 
must  be  considered  in  determining  the  overall  accuracy  of  the  calcu¬ 
lated  equilibrium  compositions,  are  the  choice  of  species  postulated 
to  be  present  in  the  system  and  the  accuracy  of  the  thermodynamic 
data  chosen  to  represent  each  species.  Choosing  an  appropriate  set 
of  species,  which  accurately  represent  the  system  at  equilibrium,  is 
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an  inherent  difficulty  in  the  calculation  of  multicomponent  equilib¬ 
rium.  A  true  calculation  of  equilibrium  in  a  given  system  must  in¬ 
clude  any  chemical  species  which  is  formed  from  any  combination  of 
the  elements  present  in  the  system.  The  compilation  of  such  a  com¬ 
plete  thermodynamic  data  base  can  represent  a  nearly  impossible  task, 
even  for  systems  consisting  of  only  a  few  elements. 

It  is  important  to  reali2e  that  anytime  a  possible  species  is 
not  included  in  the  data  set,  a  constrained  equilibrium  calculation 
will  result.  This  is  most  easily  understood  if  the  calculation  of 
chemical  equilibrium  is  considered  from  the  viewpoint  of  the  Rand 
algorithm.  In  the  Rand  algorithm,  multicomponent  equilibrium  repre¬ 
sents  the  optimal  distribution  of  a  given  quantity  of  elements  among 
a  set  of  chemical  species.  The  optimizing  function  for  a  constant 
pressure  system  is  the  minimization  of  the  total  Gibbs  energy.  There¬ 
fore,  as  the  number  of  available  chemical  species  is  decreased,  the 
elements  are  constrained  to  reside  in  a  smaller  selection  of  mole¬ 
cules.  This  causes  a  shift  in  the  calculated  compositions  in  order 
to  satisfy  the  atom  balance  while  keeping  the  Gibbs  energy  of  the 
system  as  low  as  possible.  In  general,  the  exclusion  of  a  species 
will  impact  the  equilibrium  composition  of  the  remaining  species  con¬ 
taining  similar  atoms  that  are  in  the  vicinity  or  below  the  equilib¬ 
rium  composition  of  the  excluded  species. 

The  procedure  for  developing  a  species  list  first  excludes  those 
species  not  expected  to  be  present  because  of  severe  kinetic  limita¬ 
tions.  In  practice,  this  species  set  is  developed  by  including  only 
those  species  which  have  been  experimentally  observed  in  the  system  or 
in  appropriate  subsystems.  As  an  example,  mass  spectroscopic  work  in 


the  CVD  of  GaAs  has  indicated  approximately  10  species,  but  observa¬ 
tions  in  the  subsystem  Si-H-Cl  indicate  approximately  15  additional 
species.  The  next  step  consists  of  generating  an  initial  thermodynamic 
data  base  for  all  species.  This  includes  thermodynamic  compilations 
(such  as  the  JANAF  tables),  data  bases  of  other  investigators  for  sim¬ 
ilar  systems  and  crude  estimates  for  the  remaining  species.  Initial 
equilibrium  calculations  are  then  performed  at  the  extreme  limits  of 
the  study  and  those  species  whose  compositions  are  more  than  three  or¬ 
ders  of  magnitude  in  mole  fraction  below  the  range  of  interest  are  ex¬ 
cluded.  Finally,  the  initial  thermodynamic  data  base  is  completed  by 
referring  to  the  literature  and  the  results  of  internal  consistency 
tests. 

The  sensitivity  of  the  results  to  errors  in  the  thermodynamic  data 
was  investigated  by  Smith  [70]  in  terms  of  a  Jacobian,  which  relates 
the  changes  in  the  calculated  results  to  changes  in  the  input  parameters. 
The  first  order  approximation  to  the  result  was 
N  5n. 

Mn.|  <11^1  I6u»l  3-1 

Here,  n.  is  the  number  of  moles  of  specie  i  present,  is  the  standard 
chemical  potential  of  specie  j  and  N  is  the  total  of  components.  This 
expression,  while  simple  in  form,  is  extremely  difficult  to  evaluate 
due  to  the  complicated  and  implicit  nature  of  the  function  n.(pp  for 
all  values  of  i .  If  problems  seem  apparent  for  some  species,  this 
function  can  be  numerically  evaluated.  The  work  of  Sirtl  and  Hunt  [71] 
and  similar  calculations  performed  in  Chapter  Four  showed  by  means  of 
a  calculated  example  the  effects  of  changes  in  the  enthalpy  of  formation 
of  SiHClj  on  the  predicted  equilibrium  ratio  of  SiCl ^/SiHCl 2-  This 


ratio  was  found  to  change  by  approximately  two  orders  of  magnitude 
for  a  ICZ  change  in  the  standard  enthalpy  of  formation.  The  shape 
of  the  curve  relating  this  ratio  to  tanperature  was  also  found  to 
change  markedly.  Therefore,  it  is  extremely  important. to  critically 
review  the  thermodynamic  data  set  in  order  to  perform  meaningful 
eouilibrium  calculations.  The  absolute  composition  of  the  calculated 
solution  can  be  no  better  than  the  data  set  employed.  Extreme  care 
must  also  be  used  when  comparing  calculated  equilibrium  compositions 
to  experimental  process  compositions  as  the  latter  include  possible 
kinetic  limitations.  However,  these  calculations  are  of  great  value 
in  predicting  the  directions  of  compositional  changes,  particularly 
at  the  high  temperatures  and  low  pressures  encountered  in  this  study. 

Thermodynamic  Models  of  CVD 
Models  for  the  CVD  Source  and  Pre-Source  Zones 

Each  of  the  CVD  systems  under  investigation  may  be  separated, 
based  on  temperature  or  composition,  into  the  pre-source,  source,  mixing 
and  deposition  zones.  The  pre-source  zone  was  investigated  as  a  source 
of  Si  by  considering  the  equilibrium  gas  phase  Si-content  in  the  sys¬ 
tem:  Si02(c)  in  excess,  carrier  gas  (H2  or  inert)  and  vapor  reactant 
(VH3,  VCI3  or  HCl). 

Historically,  the  chloride  process  has  employed  two  different 
aroup  III  source  materials  (III-V(c)  and  III(l))  and  thus  required 
two  separate  model  formulations.  The  III-V(c)  chloride  source  zone 
model  considered  the  system:  Si02(c)  in  excess,  III-V(c)  in  excess, 
carrier  gas  (H2  or  inert)  and  VCl3(g). 

The  chloride  process  which  used  a  group  III(l)  as  the  group  III 
source  material  was  a  more  complicated  situation.  Shaw  [23]  has 
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studied  this  source  zone  and  found  that,  following  an  initial  transient, 
a  constant  rate  of  mass  loss  of  material  occurred.  An  overall  mass 
balance  on  the  source  boat  yields  the  following  expression: 

+  My)  +  n  (xMjjj  +  (l-x)My)] 

n'^(yMiii  +  (l-y)My)  *  constant  3-2 

while  a  group  III  component  mass  balance  on  the  source  boat  produces  the 
constraint: 

^  [i  3-3 

In  these  expressions:  n*^  and  n^  are  the  moles  of  atoms  in  the  solid 
GaAs  crust  and  the  lijquid  III-V  mixture,  Mjjj  and  My  are  the  molecular 
weights  of  the  group  III  and  V  elements,  x  and  y  are  the  mole  fractions 
of  the  group  III  element  in  the  liquid  and  vapor  phase,  and  is  the 
molar  rate  at  which  vapor  species  are  formed.  If  it  is  assumed  that  the 
solid  and  liquid  phases  are  in  equilibrium,  the  liquid  phase  mole  frac¬ 
tion  is  a  function  of  temperature  only  and  therefore  constant  for  a 
given  process  condition.  Furthermore,  if  the  actual  kinetic  processes 
produce  a  steady  state  value  of  the  crust  thickness,  the  first  term  on 

the  left  side  of  both  equations  is  zero.  With  these  assumptions,  equa- 

•  y  dn^ 

tions  3-2  and  3-3  can  be  solved  to  show  that  n  =  and  x=y.  That 
is,  the  rate  at  which  group  III  and  V  atoms  are  introduced  into  the 
vapor  phase  is  equal  to  the  rate  of  loss  for  the  melt  and  the  III/V 
vapor  atom  ratio  is  the  same  as  that  in  the  melt.  Another  way  of  view¬ 
ing  the  situation  is  to  consider  the  three  phase  equilibrium  problem. 

The  activity  of  the  group  III  and  V  elements  in  the  solid  compound  can 
vary  greatly  even  though  the  stoichiometry  (~1:1)  can  be  very  small 
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and  therefore  the  sum  of  these  two  activities  is  nearly  constant.  The 
large  amount  of  melt  in  equilibrium  with  the  solid  will,  however,  fix 
the  activity  of  each  element  in  the  solid,  with  the  group  III  activity 
being  considerably  higher  than  that  of  the  group  V,  which  in  turn  fixes 
the  vapor  phase  fugacity.  In  the  event  that  the  assumption  of  constant 
crust  thickness  is  not  valid,  the  dn^/dt  terms  in  equations  3-2  and  3-3 
must  be  retained  and  the  result  is 


-  (2x-l)  ^  „  (2x-l)  dn^ 
n  -  constant  -  2(xlyj  dt  2yH  dt 


The  limit  that  =  0  implies  that  y  *  1/2  and  the  source  can  be  con¬ 
sidered  to  be  pure  solid  compound.  This  limit  is  simply  the  first  case 


examined  (solid  compound  source).  Thus,  an  investigation  of  the  two 


source  zones  described  here  should  establish  the  limits  of  operation  for 


the  liquid  source  in  the  chloride  process.  In  practice,  the  conditions 
of  operation  may  lie  somewhere  in  between  with  the  observed  III/V  ratio 
providing  an  indicator  of  the  relative  rates.  However,  if  x  is  a  con¬ 
stant  as  determined  by  the  condition  of  solid-liquid  equilibrium  and  y 
is  also  a  constant  as  witnessed  by  a  constant  growth  rate,  it  follows 
that  both  45-  and  45-  are  constant.  If  dnVdt  is  dependent  upon  the 
crust  thickness,  n^,  (i.e.  a  diffusion  limited  process),  then  it  is 
impossible  for  dnVdt  to  be  constant  for  a  finite  value  of  dn^/dt. 


which  implies  operation  at  one  of  the  limits. 

The  above  considerations  motivated  a  model  for  the  liquid  group 
III  source  zone  to  consist  of  an  ideal  vapor  phase  in  equilibrium  with 
excess  Si02(c)  and  IIIjjV^_^(l ) .  The  gas  inlet  stream  contained  VClj 
and  carrier  gas  (H2  or  inert).  The  development  of  a  thermodynamic 
data  set  for  the  hypothetical  specie  IUx'^l-x^^^  presented  in  Chapter 
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Four.  Thus,  with  this  relation,  the  compound  crust  does  not  contri¬ 
bute  elements  to  the  system. 

Two  source  zones,  one  for  the  thermal  decomposition  of  the  group 
V  hydride  and  one  for  volatilization  of  the  group  III  liquid,  are  found 
in  the  hydride  CVD  process.  The  group  V  source  zone  was  modeled  as  an 
ideal  vapor  phase  in  equilibrium  with  excess  Si02(c).  The  inlet  gas 
reactants  were  the  hydride  (VH^)  and  carrier  gas  (H2  or  inert)  at  con¬ 
stant  temperature  and  pressure.  The  group  III  source  zone  included 
excess  pure  group  III  liquid  in  equilibrium  with  excess  Si02{c)  and  an 
ideal  vapor  phase  (HCl  plus  carrier  gas). 

Models  for  the  CVD  Mixing  and  Deposition  Zones 

Since  the  only  differences  between  the  chloride  and  hydride  sys¬ 
tems  exist  in  the  source  regions,  the  mixing  and  deposition  zone  models 
were  identical  in  both  systems.  An  ideal  vapor  phase  in  equilibrium 
with  excess  Si02(c)  was  used  for  the  mixing  zone  model.  Formation  of 
solid  III-V  compound  was  postulated  to  be  kinetically  hindered  and  thus 
assumed  not  to  exist  in  the  mixing  zone.  As  a  result,  it  was  possible 
for  this  region  to  be  supersaturated.  The  model  also  allowed  the  addi¬ 
tion  of  various  species  (i.e.  HCl,  H2O,  VCl^,  VH^)  in  order  to  study 
their  effects  on  silicon  activity.  The  gas  reactant  input  for  the  mix¬ 
ing  zone  was  identical  to  that  calculated  from  the  equilibrium  source 
2one(s).  Consistent  with  the  source  and  mixing  zone  models,  the  vapor 
phase  of  the  deposition  zone  was  assumed  to  be  an  ideal  mixture.  Due 
to  the  large  volumetric  flowrate  of  gases  and  the  relatively  small 
deposition  rates  in  these  CVD  processes,  the  depletion  of  group  III, 
group  V  and  silicon  species  in  the  vapor  phase  as  a  result  of  film 
deposition  or  wall  interaction  was  neglected.  The  equilibrium  mixing 
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zone  gas  mixture  served  as  the  input  to  the  deposition  zone.  Essen¬ 
tially,  the  above  assumption  fixed  the  moles  of  each  element  in  the 
vapor  and  the  new  equilibrium  composition  was  calculated  as  a  result 
of  temperature  change  only.  This  model  provided  an  upper  bound  for 
the  computed  value  of  the  Si  activity  since  the  lower  temperatures 
found  in  the  deposition  zone  shifted  the  wall  interaction  towards 
Si02(c)  formation  and  including  the  III-V  compound  deposition  with 
Si  incorporation  would  remove  Si  from  the  gas  phase.  This  model, 
therefore,  assumes  that  the  epi-film  grows  from  a  supersaturated  vapor 
mixture  of  pseudo-steady  state  composition.  Furthermore,  this  pro¬ 
cedure  avoided  having  knowledge  of  the  solid  silicon  activity  co¬ 
efficient.  In  order  to  implement  this  model,  the  III-V  solid  phase 
was  not  included  in  the  deposition  zone,  thus  allowing  calculation  of 
the  degree  of  supersaturation  in  this  zone. 

The  effect  of  not  accounting  for  depletion  of  the  group  III  and  V 
atoms  from  the  vapor  phase  can  be  tested  by  the  following  simple  anal¬ 
ysis.  The  molar  growth  rate  of  an  epitaxial  layer  is 

V  •  3-5 

Where 

g^jj  =  molar  growth  rate  (moles/time) 
g^  =  linear  growth  rate  (length/time) 

p  =  compound  molar  density  (moles/1 enqth^) 

2 

A  *  substrate  area  (length  ) 

A  typical  set  of  operating  parameters  for  a  hydride  CVO  process 
would  specify  a  total  volumetric  flowrate  of  500  SCCM  through  each 
source  zone  having  an  inlet  composition  of  IS  HCl  to  the  group  III 
source  zone  and  IS  VH2  composition  to  the  group  V  source  zone.  Assuming 


that  all  of  the  HCl  react  to  form  III-Cl  results  in  3.7  u-moles/s 
of  group  III  atoms  transported.  The  molar  flowrate  of  group  V  atoms 
would  also  be  3.7  y-moles/s.  Choosing  as  typical  deposition  parameters 
a  2.54  cm  diameter  circular  substrate,  a  linear  growth  rate  of  1  ym/min 
and  given  the  molar  density  of  GaAs  as  0.0367  moles/cm^  [72],  the  re¬ 
sulting  molar  growth  rate  is  0.31  y-moles/s.  Thus,  in  the  worst  case, 
less  than  10%  of  the  III  and  V  atoms  are  depleted.  The  smaller  the 
growth  rate  and  substrate  surface  area  or  the  larger  the  volumetric 
flowrate,  the  better  the  approximation  becomes.  If  reaction  depletion 
was  indeed  important,  highly  non-uniform  film  thickness  would  occur. 
However,  this  is  not  experimentally  observed.  Similar  analyses  applied 
to  the  GaAs  chloride  system  and  the  analogous  InP  systems  yield  equiv¬ 
alent  results. 

The  activity  of  silicon  in  the  epi-layer  was  further  studied  in 
the  presence  and  absence  of  the  Si02  reactor  wall.  Since  the  deposition 
zone  is  typically  operated  at  a  lower  temperature  than  the  source  and 
mixing  zones,  inclusion  of  the  reactor  wall  would  be  expected  to  de¬ 
crease  the  silicon  activity  as  some  of  the  silicon  in  the  vapor  phase 
is  redeposited  on  the  reactor  wall  in  the  deposition  zone,  therefore, 
providing  an  additional  method  of  bounding  the  maximum  value  of  the 
silicon  activity  in  the  epitaxial  layer.  Justification  for  neglecting 
the  reactor  wall  lies  in  the  heterogeneous  nature  of  the  gas-wall  re¬ 
action.  Due  to  the  lower  temperature  of  the  deposition  zone  ('•873  K) , 
it  is  expected  that  this  heterogeneous  reaction  does  not  equilibrate 
as  rapidly  as  it  should  in  the  source  and  mixing  zones  (■-973  K).  This 
expectation  arises  from  the  fact  that  adsorption  reaction  rates  de¬ 
crease  strongly,  and  to  a  lesser  extent  molecular  diffusivities,  with 
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decreasing  temperature.  Additionally,  the  mean  residence  time  is 
typically  much  smaller  in  the  deposition  zone.  Thus,  the  reactor  wall 
in  the  deposition  zone  should  not  interact  with  the  vapor  phase  as 
strongly  as  it  does  in  the  source  and  mixing  zones. 

In  order  to  carry  out  parametric  analyses  of  the  two  processes, 
"base  cases"  were  chosen  for  each  system  around  which  each  parameter 
could  be  varied.  The  base  cases  were  chosen  from  commonly  used  oper¬ 
ating  parameters  reported  in  the  literature  shown  in  Tables  3-1  and 
3-2,  thus  providing  results  which  may  be  compared  to  experimental 
results.  The  chloride  system  base  case  parameters  were  as  follows. 

Source  Zone  Temperature  =  Mixing  Zone  Temperature  :  973  K 

Deposition  Zone  Temperature  :  873  K 
Inlet  VClj  Composition  :  1% 

Carrier  Gas  : 

Pressure  :  100  kPa 

For  the  hydride  system  the  following  base  case  was  chosen. 

Source  Zone  Temperature  =  Mixing  Zone  Temperature  :  973  K 

Deposition  Zone  Temperature  :  873  K 
Inlet  HCl  Concentration  *  Inlet  VH^  Concentration  :  1% 

Carrier  Gas  : 

Pressure  :  100  kPa 

Typically,  the  source  zone  of  the  hydride  system  is  operated  at 
a  higher  temperature  than  that  of  the  chloride  system  in  order  to  aug¬ 
ment  the  decomposition  of  VH^.  Due  to  the  strong  influence  of  temper¬ 
ature  on  the  species  present,  the  same  temperatures  were  used  in  both 
systems  in  order  to  provide  direct  comparison  between  the  chloride  and 
hydride  CVC  systems. 


Survey  of  Typical  Operating  Parameters  for  the  GaAs  and  InP  Chloride  Systems 
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Solid  State  Defect  Chemistry 

The  model  developed  for  the  deposition  zone  specified  that  reactant 
depletion  was  not  a  significant  phenomenon.  During  steady  state  growth, 
the  compound  film  is  exposed  to  a  vapor  phase  that  is  invariant  with 
respect  to  composition.  It  follows  directly  that  the  vapor  phase  fu- 
gacity  of  silicon  must  also  be  constant  and  therefore,  at  equilibrium, 
the  activity  of  solid  silicon  in  the  epitaxial  layer  must  be  constant. 
Thus,  prediction  of  the  amount  of  silicon  deposited  in  an  epitaxial 
layer  may  be  accomplished  by  determining  the  solid  state  silicon  con¬ 
centration  as  functions  of  the  fixed  vapor  composition  and  temperature. 

A  model  which  relates  the  point  defect  structure  of  the  solid  to  the 
solid  state  silicon  concentration  imist  therefore  be  developed  and  cou¬ 
pled  to  the  calculated  vapor  phase  composition. 

Hurle  [73-76]  has  proposed  models  for  native  defects  in  GaAs  and 
for  Te,  Sn  and  Ge  doped  GaAs.  These  models  are  extended  here  to  in¬ 
clude  the  formation  of  antisite  defects  and  to  account  for  the  doping 
of  GaAs  by  Si.  The  native  point  defect  model  allows  for  the  formation 
of  Frenkel  and  Schottky  disorder  on  both  the  group  III  and  group  V  sub¬ 
lattices.  Silicon  is  allowed  to  reside  either  on  group  III  or  group  V 
lattice  sites  or  as  an  interstitial.  Furthermore,  Si  substituted  on  a 
group  III  site  is  allowed  to  form  a  complex  with  a  group  III  vacancy 
or  a  Si  atom  on  a  group  V  site. 

A  set  of  independent  formation  reactions  for  neutral  species  is 
shown  below. 

1/2  As2(g)  +  V.  =  As.  3-6 

1/2  Asjtg)  ♦  3-7 

ASa3  *  '(  ■  ''A,  t  V^3 


3-8 


I 


» 
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Gaga  +  V.  .  Ga.  +  Vg^ 

3-9 

^®Ga  *  ''as  *  ‘^^As  ''Ga 

3-10 

0  •  ''Ga  "  ''As 

3-11 

sng)  +  Vg^  =  Sig^ 

3-12 

Si(g)  *si^3 

3-13 

Si(g)  +  V.  «  Si^ 

3-14 

The  notation  used  is  consistent  with  that  of  Hurle  [76].  As  an 
example  of  this  notation,  equation  3-6  combines  arsenic  dimer  in  the 
vapor  with  a  vacancy  on  an  interstitial  site  to  yield  an  arsenic  atom 
on  an  interstitial  site.  Thus,  in  the  solid  state,  the  subscripts  de¬ 
note  whether  an  atom  (or  vacancy)  is  residing  on  the  Ga  or  As  sublattice 
or  in  an  interstitial  location.  Equation  3-6  couples  the  point  defect 
structure  to  the  vapor  phase  via  the  formation  of  As  interstitials. 
Frenkel  defects  are  accounted  for  through  equations  3-8  and  3-9. 

Schottky  defect  formation  is  given  by  equation  3-11  and  the  formation 
of  antisite  defects  is  represented  by  equations  3-7  and  3-10.  Silicon 
incorporation  is  shown  in  equations  3-12,  3-13  and  3-14. 

Consistent  with  Hurle's  model  [76],  the  interstitial  species  are 
assumed  to  remain  electrically  neutral  while  the  other  defects  may 
ionize.  The  following  set  of  formation  reactions  may  be  written  to  re¬ 
present  the  ionized  species. 


''^Ga  *  ''^Ga 

+  e" 

3-15 

''as  *  ''is  " 

e' 

3-16 

'^^As  “  ^^As 

+  h"^ 

3-17 

''Ga  ‘  ''Ga^^ 

h" 

3-18 

^’Ga  *  ^’Ga 

+  e“ 

3-19 

^Vs  *  ^’As 

+  h^ 

3-20 

’.V-.v* 

Cv.'v'', 


The  ionization  of  native  point  defects  is  illustrated  by  equations 
3-15  through  3-18  while  the  amphoteric  behavior  of  silicon  is  represented 
by  equations  3-19  and  3-20-  Equations  3-21  and  3-22  represent  silicon 
complex  formation  and  equation  3-23  accounts  for  the  generation  and  re¬ 
combination  of  electrons  and  holes. 

All  that  is  needed  to  complete  the  model  is  to  couple  the  arsenic 
partial  pressure  to  that  of  gallium  through  the  sublimation  reaction 

GaAs(c)  =  Ga(q)  +  1/2  AS2  3-24 

and  to  write  the  electroneutrality  condition 

[h*l  »  Csi‘,]  *  [vj^]  *  [As*.]  .  [e-]  »  [Si^.]  -  [v^]  * 

3-25 

where  [  ]  denotes  concentration. 

Expressions  for  the  equilibrium  constants  for  equations  3-6  through 
3-24  may  be  written  in  the  usual  manner.  As  examples,  the  equilibrium 
constants  for  equations  3-6  and  3-7  are 

^  ■  ''As,C''=1^  "isf  3-26  ^ 

^7  ■  'isf  3-27 

where  Ya  represents  the  activity  coefficient  for  species  j  in  the  solid 
solution. 

The  model  represented  by  equations  3-6  through  3-25  may  be  simpli¬ 
fied  by  considering  the  electroneutrality  condition.  Hurle  [76]  has 
shown  that,  at  the  temperature  of  growth,  the  Frenkel  defect  on  the  ar¬ 
senic  sublattice  dominates  the  electron  concentration.  Thus,  the 


^  .’nr 
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3-28 


electroneutrality  condition  becomes 
[e-]  -  [v;^] 

and  the  electron  concentration  at  the  growth  temperature  is  given  by  the 
relation  1/2 


[e-] 


'^4'^12  1  p-1' 

K2Yas%-J  ''^2 


-1/4 
2 


3-29 


Assuming  that  the  defect  structure  of  the  epitaxial  layer  is 
"frozen-in"  at  the  growth  temperature,  the  electron  concentration  given 
by  equation  3-29  is  used  to  determine  the  concentration  of  each  defect 
present.  At  room  temperature,  the  compensation  ratio,  which  is  defined 
as  the  ratio  of  donors  to  acceptors  (Mq/Na) .  is  dominated  by  the  ionized 
silicon  impurities.  Using  the  equilibrium  constant  expressions  for 
equations  3-6  through  3-24  along  with  equation  3-29  yields 


1 


3-30 


Ga'Ga- 


where 


a  = 


'^s’^u'^lS^e'  ^Sic 


'^6'^ll'^12'^19'^23^SiA5 


b  =  '^ll'^16'^18'^2lV^SiSa 

The  dependence  of  the  compensation  ratio  upon  AS2  partial  pressure 
for  VPE  grown  GaAs  has  been  found  to  be  very  small  (Nq/Na*3).  Thus, 
consistent  with  Hurle's  results  for  Sn  and  Ge  doped  material,  equation 
3-30  indicates  that  the  dominant  acceptor  in  VPE  GaAs  is  the  Sig^Vg^ 
complex. 

An  expression  for  the  total  silicon  present  in  the  epitaxial  layer 
is  given  by 


*  [5ie.v^]  *  2Csi5.si^] 

[S'3TotaI  •  l(»''tef  *  *  »  *  "I’tej  ♦ 

*'*  9*SiI*Si 

where:  d  •  (*^|^^l'l3l<i6^'^24VsiJ;  Ife'3/f'^6^^'^n''l3''l9'*'vt  ''^Sl,  ' 

ua  As  As 

f  =  (Kg^%4K]5%,-+  Yg-^)/(K|^^K^^K^2*^^gY5,-  Yv+  ) 
ua  1  As 

g  =  (Kg^^K^3Kj^^K2QK22Ys^-+  Yg-^)/(K5^%4Y5d  54  y^+  ) 

ua  Ga  As  As 

ag^  =  activity  of  silicon  in  the  solid  phase. 


Currently,  the  thermodynamic  data  necessary  for  the  evaluation  of 
the  equilibrium  constants  and  activity  coefficients  is  not  available. 
Therefore,  a  quantitative  application  of  equation  3-32  is  not  possible. 
It  can  be  seen,  however,  that  the  Incorporation  of  Si  into  an  epitaxial 
layer  increases  with  increasing  a^.  and  depends,  in  a  complicated  way, 
upon  As2  partial  pressure. 

This  model  was  developed  for  GaAs  in  order  to  keep  the  notation 
manageable.  An  analogous  model  can,  of  course,  be  constructed  for  InP, 
which  would  yield  identical  dependencies  upon  a^.  and  the  group  V  dimer 
partial  pressure. 


CHAPTER  FOUR 

THERMOCHEMICAL  PROPERTIES 
INTRODUCTION 

Suinnarized  in  this  chapter  are  the  themiochemical  properties 
used  for  the  complex  chemical  equilibrium  analysis.  The  proper  se¬ 
lection  of  a  consistent  data  set  is  of  extreme  importance  as  a  small 
error  in  a  property  value  can  greatly  influence  the  eventual  calcu¬ 
lated  equilibrium  composition.  That  this  sensitivity  can  be  important 
is  nicely  illustrated  in  the  Si-Cl-H  subsystem  as  discussedf  later. 
Essentially,  what  is  required  for  these  calculations  is  a  means  of 
specifying  the  partial  molar  Gibbs  energy  of  each  species  believed  to 
be  present  as  a  function  of  temperature,  pressure  and  composition. 
Approximately  150  species  were  initially  examined  while  only  those 
species  that  would  be  present  at  a  mole  fraction  >  10'^^  were  included 
in  the  analysis.  The  scheme  of  representing  the  data  was  to  fix  the 
zero  enthalpy  scale  at  298  K  and  1  atm  with  the  pure  components  (stan¬ 
dard  states)  Ga(c),  As(c),  In(c),  white  P{c),  HgCg),  Si(c),  Cl2(g)  and 
enthalpy  of  formation  of  the  remaining  components  at  the 
standard  conditions  from  the  reference  components  was  then  determined. 

In  addition,  the  absolute  entropy  at  the  standard  conditions  for  each 
species  was  selected  which  allowed  a  calculation  of  the  standard  Gibbs 
energy  change  for  all  possible  reactions  at  298  K  and  1  atm.  Finally, 
knowledge  of  the  constant  pressure  heat  capacity  and  assuming  ideal  gas 
behavior  allows  the  Gibbs  energy  to  be  determined  at  any  temperature  and 
pressure.  The  gas  phase  was  assumed  to  be  a  solution  of  ideal  gases  due 


to  the  low  pressure  and  high  temperatures  investigated.  For  the  con¬ 
densed  solutions,  the  pressure  dependence  of  the  thermochemical  pro¬ 
perties  was  neglected.  However,  the  moderate  negative  deviations  from 
ideal  behavior  in  composition  dependence  for  the  liquid  source  mixtures 
were  accounted  for  and  represent  one  of  the  significant  refinements 
contained  in  these  calculations. 

Three  pieces  of  information  were  required  for  each  species  with 
the  standard  enthalpy  of  formation  being  the  most  critical,  particularly 
at  the  lower  temperatures.  It  is  for  this  quantity  that  the  most  un¬ 
certainty  exists  in  the  reported  value.  On  the  other  hand,  the  standard 
entropy  is  generally  quite  accurately  known,  either  from  low  temperature 
heat  capacity  measurements  or  spectroscopic  studies.  The  high  tempera¬ 
ture  heat  capacities  were  sometimes  estimated,  but  there  exists  a  par¬ 
tial  cancellation  of  its  effect  when  calculating  Gibbs  energy  changes. 
Presented  below  is  a  discussion  of  the  properties  selected.  It  is  noted 
that  in  many  instances  the  thermochemical  data  presented  in  the  JANAF 
tables  [77-79]  were  used  and  are  discussed  in  these  tables,  therefore 
precluding  a  discussion  here. 

Pseudo-Steady  State  Constraint  for  Liquid  Source  Boat 

It  has  been  observed  that  during  VPE  of  GaAs  and  InP  using  a  liquid 
source  boat  of  pure  group  III  metal  in  the  chloride  process,  an  initial 
transient  period  exists  in  which  the  composition  of  the  gas  stream 
leaving  the  source  region  is  a  function  of  time.  Initially,  the  pure 
metal  boat  is  dissolving  group  V  atoms  thus  producing  an  excess  of  group 
Hi  chloride.  As  the  metal  becomes  saturated  with  the  group  V  element, 

3  thin  crust  of  the  compound  is  formed  at  the  top  surface  since  the 
<iensity  of  the  compound  is  less  than  that  of  the  saturated  liquid.  It 


is  observed  that  the  crust  thickness  reaches  a  steady  state  value  and 
therefore,  from  a  simple  mass  balance,  the  vapor  phase  will  contain 
all  of  the  group  V  atoms  in  the  input  stream  plus  the  amount  of  group 
V  atoms  generated  from  the  saturated  liquid  (due  to  reaction  of  chlorine 
with  the  group  III  atom).  The  exact  amount  of  group  V  element  pro¬ 
duced  from  the  source  boat  is  therefore  a  function  of  the  temperature 
(i.e.  the  equilibrium  group  V  mole  fraction  in  the  liquid  is  a  function 
of  temperature  and  the  amount  of  group  III  atoms  leaving  depends  on  the 
form  in  which  they  leave  (e.g.  IIICl,  IIICl^,  III,  etc))  and  the  flow 
rate  (i.e.  mass  transfer  efficiency).  The  mechanism  by  which  the  group 
III  and  V  atoms  reach  the  gas/solid  interface  is  not  known  but  is  not 
required  for  the  thermodynamic  model  presented  here  as  mass  transfer 
barriers  (e.g.  the  crust)  are  assumed  not  to  be  present.  All  that  is 
required  is  to  assume  a  new  species  exists  having  a  stoichiometry 
equivalent  to  the  saturated  liquid  composition. 

The  thermodynamic  properties  of  the  hypothetical  liquid  species, 
Ai_xBx(l),  can  be  estimated  in  the  following  manner.  Letting  A  repre¬ 
sent  the  group  III  atom  and  B  the  group  V  atom,  consider  the  reaction 
sequence 


(1-x)  A(c)  »  (1-x)  A(l) 

4-1 

X  B(c)  *  X  B(l) 

4-2 

(1-x)  A(l)  +  x8(l)  » 

4-3 

3l1  occurring  at  the  source  temperature,  T.  Since  A(c)  and  B(c)  are  in 
the  pure  state,  the  Gibbs  energy  changes  for  reactions  4-1  and  4-2  are 
the  Gibbs  energies  of  formation  for  the  liquid  species  and  can  be  cal¬ 
culated  from  the  thermodynamic  sequence:  solid  element  A  or  B  is  taken 

from  T  to  its  melting  temperature  T^  or  T®,  the  solid  element  is  melted, 

mm 


the  liquid  element  is  taken  from  the  melting  temperature  to  the  ori¬ 
ginal  temperature  of  interest.  Approximating  the  heat  capacity 
difference  between  the  pure  liquid  and  pure  solid,  ACp,  as  a  constant, 
the  Gibbs  energy  change  for  reactions  4-1  and  4-2  are 

=  (1-x){AhJ;  (1  -  ^  +  acJ  (T  -  tJ;  -  T  In  C^]) }  4-4 

m  m 

and 

0  -  ^  P  -  t8  -  T  in  [^1)1  ,.5 

m  ‘m 

The  Gibbs  energy  change  of  the  third  reaction  is  simply  the  Gibbs  energy 
of  mixing.  Assuming  that  a  simple  solution  model  describes  the  liquid 
behavior  results  in 

AG^  *  (a  +  bT)  X  (1-x)  +  RT[xlnx  +  (1-x)  In(l-x)]  4-6 

where  a  and  b  are  adjustable  parameters  determined  in  conjunction  with 
solid-liquid  equilibrium  data. 

The  sum  of  reactions  4-1  to  4-3  is  the  desired  formation  reaction 
(l-x)A(c)  +  X  8(c)  =  A^_^e^(l)  4-7 

while  the  corresponding  Gibbs  energy  of  formation  of  A^  x®x^^^ 
sum  of  AG^  to  AG^.  Given  the  source  temperature,  T,  the  procedure  is 
to  first  calculate  the  liquidus  composition,  x,  from  the  implicit 
equation 

-Jb  hJS  -  R  lnC4x(l-x)]  +  bc{  -  x2  -  (l.x)^]  »  - 

m 

aCj  -  x2  -  (l-x)2]  4-8 

AB  AB 

where  and  T”  are  the  enthalpy  of  fusion  and  melting  temperature  of 
the  solid  compound  AB  and  R  is  the  gas  constant.  Once  the  equilibrium 
group  V  composition  is  determined  the  standard  Gibbs  energy  of  formation 


Property 


GaAs  Reference 


InP  Reference 


1511 

72 

1332 

71,72,74 

(kcal/mole) 

25.14 

72 

14.4 

71,72,73 

a  (cal /mole) 

4666 

see  text 

5055 

75 

b  (cal /mole  -  K) 

-8.74 

see  text 

-5.0 

75 

Table  4-2 

Thermochemical  Properties  of  the  Elements  Ga,  In,  As  and  P 
Required  for  Calculating  d6^[Ai  ^B^(l),  T] 


Property 

Gd 

Ref. 

In 

Ref. 

As 

Ref. 

P 

Ref. 

302.9 

84 

429.76 

84 

1090 

17 

313.3 

29 

4H^(kcal/mole) 

1.336 

84 

0.78 

84 

5.123 

17 

0.157 

29 

4Cp(cal/mole  - 

K)  0.27 

84 

-0.2 

84 

1.0 

estimated 

0.472 

84 
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of  can  be  calculated  from  equations  4-4  to  4-6  given  T,  x 

and  the  required  thermochemical  properties.  Tables  4-1  and  4-2  summa¬ 
rize  the  thermochemical  properties  used  in  these  calculations.  The 
adjustable  parameters  a  and  b  for  GaAs  were  determined  by  reduction  of 
the  liquidus  measurements  of  Clariou  et  al.  [80],  Hall  [81],  Kuszynski 
and  Riabeev  [82]  and  Osamura  et  al.  [83]  using  a  maximum  likelihood 


algorithm. 


The  Ga-As-Cl-H  System 


The  enthalpy  of  vaporization  of  As(c)  has  been  investigated  quite 
extensively  with  a  reported  range  [84-99]  at  298  K  of  34.4  to  38.54 
kcal/mole  corresponding  to  the  standard  enthalpy  of  formation  for 
As^(g).  The  literature  has  been  summarized  by  Hultgren  et  al.  [84] 
to  1973,  while  a  more  recent  measurement  of  Rusin  et  al .  [85]  by  a 
static  method  produced  a  value  of  H^(As^,  g,  298  K)  ®  38.14  kcal/mole. 
In  addition,  Rau  [100]  has  measured  the  total  vapor  pressure  over  solid 
and  liquid  arsenic  from  850  to  1400  K  with  a  Bourdon  gauge.  Analysis 
of  the  low  temperature  results  indicated  H^(As^,  g,  298  K)  =  37.34  + 

0.2  kcal/mole.  The  value  selected  was  37.5  kcal/mole  on  the  basis 
that  the  static  methods  are  believed  to  be  more  reliable. 

The  dissociation  enthalpy,  H^(As^  =  2As2,  298  K),  has  received 
considerable  investigation  with  early  mass  spectrometric  studies  pro¬ 
ducing  values  in  the  range  of  61.5  to  73.5  kcal/mole  [101-105].  These 
measurements  are  suspected  of  overestimating  the  As^  partial  pressure 
a  result  of  the  use  of  a  low  condensation  coefficient  for  As^.  More 
recent  determinations  using  M0AS2,  MO2ASJ,  GaAs,  InAs  and  InAs  +  InSb 
sources  with  improved  Knudsen-cell  mass-spectrometer  designs  [106-108] 
and  reduction  of  PVT  measurements  [100]  have  indicated  much  lower  values 
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for  the  dissociation  enthalpy  (54.21  -  1.5,  54.26  i  1.1,  54.2  +1.4  and 
54.8  +1.0  kcal/mole,  respectively).  The  value  selected  here  is  the 
average  of  these  four  measurements,  54.4  ^1.5  kcal/mole.  Using  the 
selected  values  of  the  standard  enthalpies,  AH°{GaAs,  c,  298  K)  =  -19.52 
kcal/mole,  AH°(As^,  g,  298  K)  ~  37.5  kcal/mole  and  the  above  dissociation 
enthalpy  yields  AH°(GaAs(c)  =  Ga(c)  +  jAs2(g).  298  K)  =  42.5  kcal/mole. 

This  can  be  compared  to  the  experimentally  determined  values  of  44.85 
[106],  45.06  [107]  and  45.4  [109]  kcal/mole. 

The  value  adopted  for  AH®(GaAs,  c,  298  K)  is  -19.52  kcal/mole  as 
determined  by  Martosudirdjo  and  Pratt  [110]  with  a  precipitation  calori¬ 
metric  technique.  This  value  can  be  compared  to  the  emf  work  of  Abbasov 
[111]  and  Sirota  and  Yushevich  [112]  in  which  values  of  -19.4  and  -20.96 
kcal/mole  were  reported,  respectively.  These  latter  two  values  are  expected 
to  include  uncertainties  due  to  the  assumed  valency  of  Ga  in  the  gal¬ 
vanic  cell  and  the  inability  to  accurately  determine  the  full  tempera¬ 
ture  dependence  over  the  relatively  small  temperature  range  of  measure¬ 
ment.  In  addition,  a  considerable  number  of  dissociation  pressure 
studies  [102-107,  113-118]  and  flow  equil ibration  investigations  with 
reactive  gases  [109,  119]  have  been  performed  which  contain  information 
about  solid  GaAs.  However,  knowledge  of  the  thermodynamic  properties  of 
other  species  is  required  (i.e.  As2(g),  As^(q),  GaCl(g),  etc.)  to  speci¬ 
fy  the  properties  of  GaAs{c).  Thus,  the  uncertainty  in  these  proper¬ 
ties  must  be  considered  in  addition  to  those  associated  with  the  measure¬ 
ments.  This  work  was,  however,  used  as  a  test  of  internal  consistency 
in  the  total  data  set.  The  standard  entropy  of  GaA*^ ,  S°(GaAs,c,298  K) 
was  taken  from  the  low  temperature  heat  capacity  measurements  of  Pies- 
bergen  [120]  while  the  high  temperature  heat  capacity  was  taken  from  the 


measurements  of  Lichter  and  Sommelet  [121]  and  are  in  good  agreement 
with  the  work  of  Cox  and  Pool  [122]  and  the  estimates  of  Marina  et  al . 
[123]  and  Maslov  and  Maslov  [124]. 

Very  little  experimental  information  is  available  for  the  arsenic 
chlorides.  The  reported  range  for  the  enthalpy  of  formation  of  AsCl^Cg) 
is  -52  to  -72  kcal/mole  [97,  99,  125-130].  The  value  adopted  was 
4H^(AsCl2,  g,  298  K)  =  -62.7  kcal/mole,  taken  from  the  enthalpy  of 
formation  of  the  liquid  and  the  enthalpy  of  vaporization.  The  enthalpy 
of  formation  of  the  mono  and  dichlorides  was  taken  from  the  estimates 
of  Shaulov  and  Mosin  [127]  as  was  the  standard  entropy  and  heat  capacity. 
The  enthalpy  of  fonmation  of  arsine  was  taken  as  aH^(AsH2,  g,  298  K)  = 
16.0  +  1.5  kcal/mole  based  on  the  work  of  Gunn  [131]  and  reported  tabu¬ 
lations  [97-99].  Finally,  the  thermochemical  properties  of  the  remain¬ 
ing  arsenic  hydrides  were  estimated  by  comparison  with  the  hydrides  of 
N,  P  and  Sb  [77-79]. 

The  thermodynamic  information  available  for  the  chlorides  of  gallium 
is  somewhat  scarce  and  inconsistent.  The  enthalpy  of  formation  for 
GaCl2(g)  was  determined  from  AH^(GaCl^,  c,  298  K)  »  -125.0  kcal/mole 
[132]  and  the  enthalpy  of  sublimation  taken  from  the  vapor  pressure 
measurements  of  Kuniya  et  al .  [133],  fiH^CGaCl^,  g,  298  K)  =  -107.3  +  3 
kcal/mole.  The  enthalpy  of  formation  of  gallium  monochloride  has  the 
reported  value  AH^(GaCl ,  g,  298  K)  =  -19.5  kcal/mole  and  is  taken  from 
the  dissociation  energy  of  Barrow  [134],  However,  a  value  of  AH^(GaCl , 
g,  298  K)  =  -12.0  kcal/mole  is  obtained  using  AH^(GaAs(c)  +  HCl(g)  * 

Gael (g)  +  j  As^(g)  +  ^  H2(g) ,  950  K)  =  37.52  +  8  kcal/mole  determined 
by  Battat  et  al .  [119]  and  the  thermochemical  data  selected  for  the  other 
species.  These  results  are  in  sharp  contrast  to  the  vapor  pressure 


57 


measurements  of  Kuniya  et  al .  [133]  who  report  a  second  law  AH“(GaCl3(g) 
=  GaCl(g)  +  Cl2(g),  1083  K)  =  45.912  kcal/mole.  The  value  selected  was 
AH^(GaCl ,  g,  298  K)  =  -17.0  5  kcal/mole  based  on  the  first  two  re¬ 

ports,  considering  the  value  of  Barrow  [134]  slightly  more  due  to  the 
uncertainties  found  in  the  enthalpy  of  formation  for  the  other  species 
in  the  reaction  studied  by  Battat  et  al .  [119].  The  enthalpy  of  forma¬ 
tion  for  gallium  dichloride  was  taken  from  the  measurements  of  Battat 
et  al.  [119]  using  the  thermochemical  data  selected  here  and  correcting 
the  second  law  entropy  to  that  calculated  by  Shaulov  and  Mosin  [135]. 

The  enthalpy  of  dimerization  has  been  investigated  by  several  authors 
[133,  136-139]  with  the  reported  enthalpy  and  entropy  of  dimerization 
in  the  range,  AH“(2GaCl3(g)  =  Ga2Clg(g),  298  K)  =  -22.6  to  -20.9  kcal/ 
mole  and  S‘’(2GaCl3(g)  =  GagClgCg),  298  K)  =  -31.66  to  -36.0  cal/mole-K. 
Accepting  the  enthalpy  and  entropy  of  dimerization  as  -21.0  kcal/mole 
and  -33.0  cal/mole-K,  respectively,  and  combining  these  results  with 
the  selected  thermochemical  data  for  GaCl3  produces  AH^(Ga2Clg,  g,  298 
K)  =  -235.6  +  10  kcal/mole  and  S°(Ga2Clg,  g,  298  K)  =  127  +  6.0  cal/ 
mole-K.  The  standard  entropy  and  heat  capacity  for  GaCl ,  GaCl2  and 
6aCl3  were  taken  from  Shaulov  and  Mosin  [135]  while  the  heat  capacity 
of  Ga2Clg  was  approximated  by  the  value  for  Al2Clg  [77].  In  addition, 
other  species  are  expected  to  exist  (i.e.  Ga2Cl^,  Ga2Cl2)  [lAO,  141], 
but  no  thermochemical  data  is  available. 

The  In-P-Cl-H  System 

The  standard  entropy  at  298  K  and  the  constant  pressure  heat  capac¬ 
ity  of  solid  and  vapor  In  were  taken  from  Hultgren  et  al .  [84],  As  sum- 
marized  by  Hultgren  et  al .  [84],  the  standard  enthalpy  of  vaporization 
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to  the  vapor  pressure  measurements  produces  the  range  of  49.8  to  58.1 
kcal/mole  for  AH^(In,  g,  298  K).  More  recent  mass  spectrometric  re¬ 
sults  of  Panish  and  Arthur  [132]  and  Farrow  [142]  suggest  the  values 
of  56.6  and  58.03  kcal/mole,  respectively,  with  the  average  of  these 
two  values  adopted  here.  In  a  similar  fashion,  the  thermochemical 
properties  of  phosphorous  selected  by  Hultgren  et  al.  [84]  or  the 
JANAF  tables  [77]  are  in  agreement  with  the  more  recent  results  [132] 
and  were  adopted  for  this  study.  However,  there  exists  a  small  dif¬ 
ference  in  the  reported  aH‘*(298  K)  of  the  reaction;  P^(g)  =  2P2(g)- 
Foxon  et  al .  [143]  report  a  value  of  aH“(298  K)  =  57.9  +  1  kcal/mole 
while  Panish  and  Arthur  [132]  reported  aH®{298  K)  ■  53.8  +  0.8  kcal/ 
mole  from  third  law  calculations  of  their  mass  spectrometric  results. 
The  third  law  reduction  of  the  mass  spectrometric  results  of  Farrow 
[142]  produces  a  value  of  ZiH®(298  K)  «  58.04  +  0.3  kcal/mole,  while 
the  JANAF  tables  [77]  suggest  AH®(298  K)  »  54.6  +  1.1  kcal/mole.  An 
average  value  was  selected  of  AH®(298  K)  *  56.1  +  2.0  kcal/mole. 

A  rather  wide  range  in  the  reported  values  for  the  standard  en¬ 
thalpy  of  formation  of  solid  InP  exists  (-13.52  to  -22.3  kcal/mole). 

As  shown  in  Table  4-3,  the  value  selected  was  the  average  of  the  two 
solution  calorimetric  determinations  as  thu  is  a  direct  determination 
of  the  property.  The  results  from  the  vapor  pressure  measurements  are 
subject  to  uncertainties  in  the  properties  of  the  vapor  phase  species 
and  also  the  heat  capacity  of  solid  InP  (e.g.  Panish  and  Arthur  [132] 
used  Cp  for  AlSb  which  produces  a  decrease  in  AH^(InP,  c,  298  K)  of 
0.5  kcal/mole,  when  compared  to  Cp(InP,  c,  T)  of  Pankratz  [144]).  The 
standard  entropy  of  InP(c)  was  taken  from  the  low  temperature  heat 
capacity  measurements  of  Piesbergen  [120]  while  the  heat  capacity 
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Table  4-3 

The  Reported  Standard  Enthalpy  of  Fonriation  of 


InP(c).  AH°(InP.  c.  298  K) 


298^^*^^^^^^  Method  Reference 

kcal  /tnol  e _ 


-18.83  +  O.tJ^ 
-18.58  +  0.7^“^ 

flow  equilibration 

132 

mass  spectrometry 

132 

-22.3  +  1.5-,  X 
-19.33  +  0.1 
-17.83  +  1 .4^'’^ 

mass  spectrometry 

104 

vapor  pressure 

109,  145 

calculated 

109 

(e) 

-20.33'®' 

calculated 

146 

mass  spectrometry 

142,143 

-13.52  +  0.26 

solution  calorimetry 

no 

-21.0  +  2 

bomb  calorimetry 

147 

-21.5  +  1.5 

bomb  calorimetry 

144,  148 

-14.5  +  0.44 

solution  calorimetry 

referenced  in 

no 

^  (a)  InP(c)  »  In(c)  +7  p2{g).  H-gg  =  360  kcal/mole 

*  298^^2'^^^  taken  as  34,34  kcal/mole 

S;  (b)  InP(c)  =  In(c)  +^p4(a).  Higg  =  22.1  kcal/mole 

298^^4^^^^  taken  as  14.08  kcal/mole 

B  (c)  InP(c)  =  In(c)  +2  P2(g)»  H^gg  =  36,5  kcal/mole 

298^**2^^^^  taken  as  34.34  kcal/mole 

(d)  InP(c)  =  In(c)  +  \  P2(g),  H|gg  =  35.0  kcal/mole 
298^'^2^^^^  taken  as  34.34  kcal/mole 

r  (e)  InP(c)  =  In(c)  +-|p2(g).  H“gg  =  37.50  +  0.1  kcal/mole 

i;  -  298^*^2^^^^  taken  as  34.34  kcal/mole 
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adopted  was  that  measured  by  Pankratz  [144].  The  result  is  in  good 
agreement  with  the  298  K  value  of  Piesbergen  [120]  and  in  fair  agree¬ 
ment  with  the  high  temperature  measurements  of  Cox  and  Pool  [122] 
and  the  suggested  value  of  Maslov  and  Maslov  [124]. 

Barrow  [134]  reports  the  dissociation  energy  of  InCl  to  be  102.5 
kcal/mole  and  combining  this  with  the  value  of  the  enthalpy  of  forma¬ 
tion  of  In(g)  and  Cl(g)  yields  AH^(InCl,  g,  298  K)  =  -16.21  kcal/mole. 
However,  the  atomic  fluorescence  value  for  the  dissociation  energy 
also  reported  by  Barrow  [134]  (0^  *  104.6  kcal/mole)  yields  AH^(InCl, 
g,  298  K)  *  -18,31  kcal/mole.  Klemm  and  Brautigan  [149]  reported 
AH^(InCl,  c,  273  K)  =  -44.6  kcal/mole  and  when  combined  with  the 
enthalpy  of  sublimation,  4H^(InCl,  c,  298  K)  =  27.8  kcal/mole  [150] 
gives  4H^(InCl ,  g,  298  K)  =*  -16.8  kcal/mole  and  is  the  value  adopted 
here.  The  standard  entropy  of  InCl  was  taken  from  the  calculations 
of  Malakova  and  Pashinkin  [151]  while  the  heat  capacity  is  that 
recommended  by  Kelly  [152].  The  standard  enthalpy  of  formation  and 
entropy  of  InCl 2  was  taken  from  the  estimate  of  Glassner  [153]  and  the 
heat  capacity  is  the  same  as  that  listed  for  GaCl2.  The  standard  en¬ 
thalpy  of  formation  and  entropy  of  GaCU  was  taken  to  be  the  values 
suggested  by  Mull  in  and  Hurle  [97]  and  the  constant  pressure  he4t  ca¬ 
pacity  estimated  by  Shaw  [98].  The  thermochemical  properties  of  the 
dimer,  In2Clg,  were  taken  from  the  values  suggested  by  Schafer  and 
Binnewies  [139]. 

The  standard  enthalpy  of  formation  of  phosphine  was  taken  from 
the  decomposition  studies  of  Gunn  and  Green  [154]  and  the  remaining 
properties  from  the  JANAF  tables  [77,  79].  The  JAMAF  tables  were  also 
used  for  the  other  phosphorous  hydrides,  chlorides  and  oxide  vapor 
phase  species. 
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The  thennocheniical  properties  of  Si  have  been  reviewed  by  Hultgren 
et  al .  [84]  and  the  JANAF  tables  [77].  In  particular,  there  exists  a 
rather  large  range  in  the  reported  third  law  values  of  the  standard  en¬ 
thalpy  of  vaporization,  AH^(Si,  g.  298  K)  =  86.75  to  109.06  kcal/mole. 
The  value  selected  was  in  between  the  Knudsen  studies  of  Oavis  et  al. 
[155]  and  Grieveson  and  Alcock  [156]. 

The  Si-Cl-H  system  has  received  considerable  attention  due  to  its 
importance  in  the  semiconductor  industry.  An  excellent  review  of  the 
literature  for  this  system  with  equilibrium  calculations  presented  is 
given  by  Hunt  and  Sirtl  [157]  and  Sirtl  et  al .  [158].  The  posture  taken 
here  is  to  assume  that  SiCl^  has  the  most  reliable  thermodynamic  data 
with  these  values  being  fixed  by  the  JANAF  tables  [77].  The  reaction 
Si(g)  +  SiCl^(g)  *  2SiCl2(g)  4-9 

has  been  investigated  extensively  [159-164].  Employing  the  thermo¬ 
chemical  data  for  the  three  species  in  reaction  4-1  from  the  JANAF 
tables  [77],  third  law  values  of  AH^(SiCl2.  9.  298  K)  were  calculated 
from  the  experimental  data.  The  effusion-mass  spectrometric  determin¬ 
ation  of  Farber  and  Srivastava  [164]  over  the  temperature  range  of 
1593  K  to  1792  K  produced  the  value,  AHp(SiCl2.  9.  298  K)  =  -40'.39  + 

0.3  kcal/mole  and  showed  no  temperature  dependence.  This  result  is 
in  good  agreement  with  flow  equilibration  data  of  Schafer  et  al .  [159] 
(1273  K  to  1473  K) ,  Teichmann  and  Wolf  [160]  (1223  K  to  1575  K)  and  the 
static  measurements  of  Schafer  and  Nicki  [163]  (1398  K  to  1573  K) ,  with 
third  law  values  of  -40.62,  -40.54  and  -40.44  kcal/mole,  respectively. 
The  flow  equilibrium  values  of  Antipin  and  Sergeev  [161]  (1273  K  to 
1673  K)  and  the  static  values  of  Ishino  et  al .  [162]  (1443  K  to  1573  K) 
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were  more  negative  and  exhibited  a  marked  temperature  dependence.  On 
the  basis  of  these  calculations,  the  value  AH|(SiCl2t  9t  298  K)  * 

-40.4  kcal/mole  was  selected.  The  values  for  the  standard  enthalpies 
of  formation  of  the  less  stable  chlorine  SiCl  and  SiCl^  were  computed 
from  the  high  temperature  flow  equilibrium  studies  of  Farber  and  Srivas- 
tava  [164].  In  the  third  law  analysis,  the  data  previously  discussed 
was  used  in  conjunction  with  the  heat  capacity  for  SiCl  and  SiClj. 
suggested  by  the  JANAF  tables  [77]  and  produced  the  value  of  47.4  +  0.6 
and  -93.3  +  0.5  kcal/mole,  respectively.  These  results  are  in  agreement 
with  the  analysis  of  Rusin  et  al.  [165-167]  on  the  total  pressure  meas¬ 
urements  of  Schafer  and  Nicki  [163].  No  additional  thermodynamic 
studies  of  Si2Clg  are  known  to  exist  and  thus  the  properties  suggested 
by  Hunt  and  Sirtl  [157]  were  adopted. 

It  was  pointed  out  by  Hunt  and  Sirtl  [157]  that  the  mole  ratio  of 
SiCl^  to  SiHCl^  is  very  sensitive  to  the  assumed  value  of  the  standard 
enthalpy  of  formation  of  SiHCl^.  Indeed,  this  ratio  is  seen  to  vary 
by  nearly  four  orders  of  magnitude  at  1000  K  when  bounded  by  the  ex¬ 
perimental  determinations  of  iH^CSiHCl^,  g,  298  K).  Since  the  work  of 
Hunt  and  Sirtl  [157]  was  published,  two  additional  experimental  in¬ 
vestigations  of  the  thermodynamic  properties  of  SiHCl^  have  been  per¬ 
formed.  Farber  and  Srivastava  [168],  from  effusion-mass  spectrometric 
measurements,  determined  the  reaction  enthalpy  for 

SiCl^(g)  +  H2(g)  =«  SiHCl3(g)  +  HCl(g)  4-10 

over  the  temperature  range  1155  K  to  1500  K.  Employing  the  thermo¬ 
dynamic  data  listed  in  Table  4-4  and  these  results,  a  relative  temper¬ 
ature  insensitive  third  law  value  for  ^H^CSiHCl^,  g,  298  K)  =  -119.30 
1  1.0  kcal/mole  is  obtained.  Using  both  static  and  dynamic  methods. 


Wolf  and  Teichmann  [169]  investigated  reaction  4-10  and  the  reaction 

4SiHCl3(g)  «  3SiCl4(g)  +  Si(c)  +  2H2(g)  4-11 

Third  law  values  for  dH^(SiHCl3,  g,  298  K)  were  calculated  from  the 
original  results  of  these  authors.  The  values  obtained  for  reaction 
4-11  and  for  the  three  data  sets  with  reaction  4-10  were  -119.47  + 

0.4,  -119.83  +  0.9,  -119.58  +  0.2  and  -119.50  ^  0.6  kcal/mole  and  the 
results  are  seen  to  be  in  good  agreement  with  the  measurements  of 
Farber  and  Srivastava  [168].  Since  these  values  were  nearly  3  kcal/mole 
more  negative  than  those  developed  by  Hunt  and  Sirtl  [157],  values  of 
^^(SiHCl^,  g,  298  K)  were  calculated  for  various  experimental  SiCl^/ 
SiHCl^  ratios  in  a  fashion  similar  to  Hunt  and  Sirtl,  The  experimental 
data  consisted  of  a  variety  of  feed  mixtures  (e.g.  SiCl^/H2,  H2/HCI , 
SiHCl2/H2)  which  were  contacted  with  Si(c)  at  different  temperatures 
during  a  deposition/etching  process.  The  results  of  these  calculations 
for  14  data  sets  suggested  /iH*(SiHCl3,  g,  298  K)  =  -118.16  +  1.70  kcal/ 
mole.  Based  on  these  results  and  the  new  experimental  determinations, 
the  value  adopted  was  AH^(SiHCl3,  g,  298  K)  *  -119.5  +  1.5  kcal/mole. 

The  standard  enthalpy  of  formation  of  the  di-  and  mono-chloro- 
silanes  was  taken  from  the  recent  measurements  of  Farber  and  Srivastava 
[168].  In  order  to  obtain  a  consistent  data  set,  third  law  values  of 
these  quantities  were  calculated  from  the  original  experimental  data 
while  using  the  data  base  adopted  here.  The  adopted  values  were 
^^(SiH2Cl2t  g,  298  K)  =  -75.5  +  2  kcal/mole  and  aH“(SiH3Cl,  g,  298  K) 

*  -32.7  +2.5  kcal/mole. 

No  additional  experimental  information  on  the  thermochemistry  of 
SiH^  and  SiH  exists  and  thus  the  JANAF  tables  recommendation  was  adopted. 
The  standard  enthalpy  of  formation  of  disilane  was  taken  from  the 
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calculations  of  Potzinger  et  al .  [170],  AH^(Si2Hg,  g,  298  K)  *  17.1  + 

3  kcal/mole  and  is  compared  with  the  calculations  of  O'Neal  and  Ring 
[171],  AH^(Si2Hg,  g,  298  K)  »  19.1  kcal/mole  and  the  value  of  16.0  kcal/ 
mole  obtained  from  the  estimated  enthalpy  of  formation  for  SiH^  (35 
kcal/mole)  and  Si-Si  (-54  kcal/mole).  The  standard  entropy  and  heat 
capacity  of  SigHg  were  obtained  by  comparison  with  C2Hg- 

A  summary  of  the  selected  thermochemical  properties  (with  two 
additional  references  [172,  173])  is  presented  in  Table  4-4.  In  or¬ 
der  that  the  stability  of  the  various  species  might  easily  be  compared, 
the  standard  molar  Gibbs  energies  of  formation  are  plotted  in  Figures 
4-1  through  4-5. 


Selected  Therciochemlcal  Values 
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Table  4-4  -  continued 
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CHAPTER  FIVE 

EXPERIMENTAL  METHOD  FOR  THE  INVESTIGATION 
OF  GROUP  V  TRIHYDRIOE  DECOMPOSITION 

Introduction 

The  decomposition  of  the  tri hydride  of  the  group  V  elements  N,  P 
and  As  has,  in  the  past,  been  studied  using  manometric  methods  in 
closed  systems  [50-55,  57-59]  and  also  by  infrared  spectrometry  in 
open  systems  [56,  60,  61,  64-66].  A  major  disadvantage  which  is  en¬ 
countered  with  measurements  based  on  manometric  techniques  is  that  the 
system  total  pressure  depends  on  all  of  the  species  present.  It  is 
therefore  difficult  to  remove  the  influence  of  other  reactions  in  the 
system  from  the  observed  data.  Usually  this  difficulty  is  addressed 
by  assuming  that  a  single  reaction  step  is  rate  limiting  and  that  the 
remaining  reaction  products  are  at  equilibrium.  This  technique  has 
been  employed  for  the  decomposition  of  PH^  and  AsH^  by  employing  the 
overall  reaction  [57] 

4VH3  —  >  V4  +  6H2  5-1 

This  reaction  is  not  applicable  to  the  NH^  system,  however,  since  no 
known  tetramers  of  N  exist.  For  NH3  decomposition,  the  following  over¬ 
all  reaction  has  been  applied  [51] 

2VH3  — >  Vg  +  3H2  5-2 

The  existence  of  one  or  more  slow  reactions  in  the  sequence  between  the 
disappearance  of  VH3  and  the  formation  of  V2  or  V^  can  cause  the  ini¬ 
tiating  step  of  the  reaction  to  appear  to  be  slov^er  than  it  actually  is. 
This  results  in  an  over  estimation  of  the  activation  energy  associated 


V^•^ 


£jiiJ 


74 

with  the  reaction.  A  major  advantage  inherent  in  manometric  methods  is 
that  pressure  can  be  measured  to  a  very  high  degree  of  accuracy. 

Spectrometric  investigations  provide  a  means  for  measuring,  direct¬ 
ly,  the  rate  of  appearance  or  disappearance  of  individual  chemical 
species  in  the  system.  Frequently,  more  than  one  signal  can  be  moni¬ 
tored  during  the  course  of  the  experiment  and,  therefore,  the  opportunity 
for  determining  the  entire  kinetic  sequence  is  greatly  enhanced.  The 
sensitivity  of  spectrometer  instruments  varies  depending  on  the  type  of 
spectrometer  and  supporting  equipment  employed,  but  it  is  not  unusual  to 
find  mass  spectrometers  which  have  detection  limits  below  1  ppm. 

A  major  disadvantage  connected  with  quantitative  composition  meas¬ 
urements  based  on  spectrometric  instruments  is  that  of  calibration.  Most 
spectrometers  provide  an  output  signal  which  is  proportional  to  the 
amount  of  the  species  which  is  present  at  the  detector.  The  value  of  the 
proportionality  constant  is  rarely  known  and  generally  depends  on  the 
specific  chemical  species  and  the  energy  (i.e.  temperature)  of  that 
species.  Moreover,  if  a  sample  must  be  removed  from  the  system  for 
analysis,  a  method  of  sampling  must  be  chosen,  such  that  the  sample 
composition  accurately  represents  the  system  from  which  it  was  removed. 
Also,  the  sampling  technique  must  either  not  perturb  the  system  signi¬ 
ficantly,  or  it  must  perturb  the  system  in  a  way  which  is  known  and 
can  be  corrected  for  during  the  data  reduction. 

Even  though  there  are  many  variables  regarding  the  application  of 
spectrometric  techniques  for  the  measurement  of  composition,  these 
techniques  are  highly  desirable  because  of  the  ability  to  follow  sig¬ 
nals  representative  of  individual  chemical  species.  Thus,  the  system 
employed  here  for  the  investigation  of  NH^,  PH^  and  AsH^  decomposition 
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is  based  on  a  quadrupole  mass  spectrometer  coupled  to  a  constant 
volume  reaction  tube  through  a  sampling  orifice.  This  technique 
provided  real  time  monitoring  of  the  reaction  gas  phase  composition 
while  perturbing  the  reacting  system  in  a  known  manner  which  was 
easily  corrected  for  during  data  analysis. 

Experimental  Apparatus  and  Method 

A  schematic  representation  of  the  equipment  used  during  the  in¬ 
vestigation  of  NH^,  PH2  and  AsHj  thermal  decomposition  is  shown  in 
Figure  5-1.  Due  to  the  extremely  toxic  and  flammable  nature  of  the 
gases  involved,  hooded  enclosures  were  constructed  around  the  storage 
area  for  the  gas  cylinders  and  the  reaction  chamber.  These  enclosures, 
along  with  the  vacuum  pump  exhaust  from  the  gas  sampling  system,  were 
vented  through  the  laboratory  exhaust  hood.  The  exhaust  gas  from  the 
reaction  tube  was  first  passed  through  a  Draeger  class  B3-P  filter 
before  being  vented  through  the  laboratory  exhaust  hood.  This  was  done 
to  remove  any  residual  VH^  and  its  toxic  reaction  products  from  the 
vented  gases. 

High  pressure  gas  cylinders  containing  4.3S  NH^,  10.07%  PH^  and 
10.03%  AsHj  in  Hg  were  connected  to  a  common  stainless  steel  gas  line 
and  solenoid  operated  valve  for  inlet  to  the  reaction  tube.  Hydrogen 
was  provided  as  an  additional  inlet  to  the  reaction  tube  through  a 
separate  gas  line  and  valve.  This  arrangement  allowed  the  reactor  to 
be  charged  with  any  of  the  available  gases  and  allowed  the  gas  in  the 
reaction  chamber  to  be  diluted  with  H2  if  desired.  Purge  gas  consisting 
of  or  He  was  available  through  the  AsH^^,  PH^  or  H2  purge  assemblies 
and  gas  lines.  A  capacitance  manometer  was  used  to  monitor  the  pressure 
in  the  reaction  tube  over  a  range  of  100  Pa  to  10^  Pa  with  a  precision 
of  0.1%. 
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The  reaction  tube  consisted  of  a  54  cm  long  by  6  cm  O.D.  quartz 
tube  placed  inside  of  a  three  zone  Marshall  furnace  (model  1169-5)  and 
is  shown  in  greater  detail  in  Figure  5-2.  A  temperature  profile 
which  was  constant  to  within  2  K  across  the  length  of  the  reaction 
tube  was  obtained  by  placing  the  tube  inside  of  two  Dynatherm  liquid 
sodium  furnace  liners  and  by  controlling  each  of  the  electrically 
heated  furnace  zones  with  individual  Linberg  Model  59344  heater  con¬ 
trollers.  The  flattest  temperature  profiles  were  obtained  when  furnace 
zones  2  and  3  were  operated  with  identical  setpoints  and  zone  1  was 
operated  at  a  setpoint  20  K  below  that  of  the  other  two  zones. 

Thermocouples  were  located  at  positions  12  cm,  25  cm,  50  cm  and  62 
cm  into  the  furnace.  These  positions  were  chosen  based  on  previous 
measurements  which  demonstrated  that  the  temperature  profile  from  16  cm 
through  54  cm  was  flat  to  within  0.5  K.  At  positions  less  than  16  cm 
into  the  furnace,  the  temperature  drops  off  due  to  heat  losses  near 
this  end  of  the  furnace.  Positions  from  54  cm  to  the  end  of  the  re¬ 
action  tube  typically  showed  the  highest  temperatures  in  the  furnace, 
but  this  was  compensated  for  by  lowering  the  setpoint  in  zone  1.  Fire¬ 
brick  insulation  was  placed  at  each  end  of  the  furnace  in  order  to 
minimize  the  heat  losses.  The  system  was  limited  to  an  operating  range 
of  673  K  to  1273  K  due  to  functional  and  safety  constraints  imposed  by 
the  furnace  liners. 

The  inlet  tube  to  the  reactor  consisted  of  a  40  cm  long,  6  mm  o.d. 
quartz  tube  which  extended  outside  of  the  furnace  and  was  mated  to  a 
stainless  steel  tube  through  a  stainless  steel  fitting  using  viton  o- 
rings  as  seals.  Sampling  of  the  gas  in  the  reaction  tube  was  accom¬ 
plished  by  continuously  drawing  a  sample  through  a  small  orifice 


(nominally  0.1  mm  in  diameter)  at  the  end  of  the  reaction  tube. 

Further  discussion  of  the  gas  sampling  system  and  mass  spectrometer 
is  provided  in  Appendix  B. 

The  use  of  an  orifice  for  obtaining  continuous  gas  samples  from 
the  reaction  tube  represents  a  significant  perturbation  on  the  re¬ 
acting  system  and,  therefore,  must  be  considered  in  the  data  analysis- 
As  is  described  in  Appendix  B,  the  mass  spectrometer  provides  an  output 
signal  which  is  proportional  to  the  partial  pressure  of  the  chemical 
species  present  in  the  reaction  chamber.  The  partial  pressure  of  each 
species  changes  due  to  participation  in  chemical  reactions  and  the 
continual  bleed  on  the  reactor  caused  by  the  gas  sampling  system.  A 
species  balance  on  the  reaction  volume  yields  the  following  equation 


r 

^i  j- 


5-3 


where:  =  pressure  of  species  i 

r^  =  chemical  reaction  rate  for  species  i 
n-  =  rate  of  molar  loss  of  species  i  through  the  sampling 
orifice 


V  =  system  volume 

Coulson  et  al .  [174]  have  analyzed  the  flow  of  a  compressible  fluid 
through  an  orifice  and  have  shown  that  for  an  isentropic  process,  crit¬ 
ical  flow  is  obtained  when 


Pd  o  k/(k-l) 

r  ■ 
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where  P^  is  the  downstream  pressure,  P^  is  the  upstream  pressure  and  k 
is  the  ratio  of  the  heat  capacity  at  constant  pressure  to  the  heat 


capacity  at  constant  volume.  The  heat  capacity  ratio  ranges  from  1.13 
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(PHj)  to  1.66  (He)  for  the  gases  used  in  these  experiments.  Therefore, 
keeping  the  ratio  of  the  reaction  tube  pressure  to  the  first  vacuum 
stage  pressure  below  0.49  will  assure  that  all  of  the  gases  flowing 
through  the  sample  orifice  are  at  critical  flow.  During  the  experi¬ 
ments,  the  partial  pressure  of  any  measurable  gas  in  the  reaction 
chamber  was  greater  than  500  Pa.  As  is  described  in  Appendix  B,  the 
operating  pressure  of  the  first  vacuum  stage  was  between  10  and  50  Pa. 
Thus,  the  requirement  for  critical  flow  was  always  met. 

Applying  the  results  of  Coulson  et  al.  [174]  to  equation  5-3 
yields  the  following  relationship 


The  variable  A^,  which  represents  the  product  of  the  orifice  dia¬ 
meter  and  discharge  coefficient,  is  unknown  since  a  sufficiently  accu¬ 
rate  value  for  the  orifice  diameter  is  not  known  and  the  discharge  co¬ 
efficient  is  a  function  of  the  orifice  Reynolds  number.  A  relationship 
between  A^,  gas  molecular  weight  and  system  temperature  was  therefore 
determined  experimentally  using  H2,  He,  N2  and  CO2. 

It  has  long  been  recognized  that  catalytic  surfaces  become  less 
active  the  longer  they  are  in  use,  but  usually  exhibit  a  relatively 
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constant  period  of  activity  following  an  initial  period  of  deactivation. 
The  quartz  reactor  surface  was  pretreated  prior  to  the  experiments  by 
pressurizing  the  tube  to  130000  Pa  with  hydrogen  at  a  temperature  of 
1073  K  for  48  hours.  The  NH2  and  AsH^  decomposition  experiments  were 
repeated  several  times  (days  apart)  in  order  to  look  for  changes  in  the 
catalytic  activity.  No  changes  in  activity  were  observed.  Between  ex¬ 
periments,  the  reaction  tube  was  either  maintained  under  vacuum  or  a 
helium  pressure  of  approximately  10^  Pa. 

The  procedure  employed  for  the  collection  of  rate  data  was  the  same 
for  each  group  V  trihydride.  Each  trihydride  was  investigated  over  a 
full  range  of  temperatures  before  the  next  was  admitted  to  the  system. 
Ammonia  was  studied  first  followed  by  phosphine  and  finally  arsine. 

The  desire  to  obtain  a  temperature  profile  which  was  as  flat  as 
possible  along  the  length  of  the  reaction  tube  required  the  suppression 
of  all  of  the  heat  losses  in  the  system.  Thus,  increasing  the  temper¬ 
ature  of  the  reactor  could  be  accomplished  quickly  (-ZOO  K/hr)  while 
decreasing  the  reactor  temperature  was  a  slow  process.  Therefore,  most 
of  the  data  taken  for  each  tri hydride  was  acquired  in  increasing  order 
of  temperature.  This  procedure  was  not  strictly  adhered  to  however, 
since  eventually,  the  next  higher  temperature  investigated  resulted  in 
reaction  rates  too  fast  to  be  followed  with  the  current  mass  spectro- 
meter  configuration.  It  was  therefore  not  unusual  to  allow  the  reactor 
to  cool  so  that  rate  data  at  intermediate  temperatures  could  be  gathered. 
This  occurred  most  frequently  for  AsH^  and  occasionally  for  PH^. 

The  laboratory  did  not  have  the  capability  for  automated  control 
0^  the  mass  spectrometer.  Each  scan  of  mass  to  charge  ratio,  there¬ 
fore,  vias  initiated  manually.  The  fastest  rate  at  which  data  could  be 
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reliably  scanned  and  averaged  by  the  mass  spectrometer  and  chart  re¬ 
corder  was  one  scan  every  10  seconds.  This  made  it  difficult  to  follow 
reactions  with  half  lives  less  than  30  seconds  since  relatively  few 
data  points  could  be  collected  before  the  species  signal  was  comparable 
to  the  background  signal  or  instrument  noise. 

The  decomposition  of  NH2  was  monitored  by  following  the  NH^  peak 
at  m/e  =  16.  This  peak  was  followed  rather  than  since  the  desorp¬ 
tion  of  HgO  off  the  walls  of  the  second  vacuum  stage  in  the  mass  spec¬ 
trometer  caused  a  large  OH*  background  peak  to  be  present  at  m/e  =  17. 

The  decomposition  of  PH^  and  AsH^  was  followed  by  observing  the 
entire  fragment  ion  pattern  V*,  VH*.  VH^  and  VH^  at  m/e  =  31,  32,  33 
and  34  for  PH^  and  m/e  =  75,  76,  77  and  78  for  AsHj-  No  significant 
background  was  observed  at  these  mass  to  charge  ratios.  The  levels 
of  confidence  for  the  PH^  and  AsH^  results  are,  therefore,  much  higher 
than  that  for  NH^. 

The  procedure  employed  for  these  experiments  was  to  first  evacuate 
the  reaction  tube  and  then  bring  the  reactor  to  the  desired  temperature. 
The  reactor  was  then  charged  with  the  desired  tri hydride  by  opening  the 
appropriate  solenoid  operated  valves  and  monitoring  the  system  pressure 
by  means  of  the  capacitance  manometer.  The  system  was  charged  to  a 
pressure  of  1.3x10^  +  10^  Pa  for  the  experiments  and  9.2x10^  +  10^ 
Pa  for  the  PH^  and  AsH^  experiments.  The  amount  of  time  required  to 
charge  the  reactor  was  between  3  and  5  seconds.  The  first  mass  spec¬ 
trometer  scan  was  begun  15  seconds  after  the  valves  were  closed.  This 
delay  was  primarily  due  to  the  restrictions  imposed  by  manual  operation 
of  the  system  (time  required  to  close  valves,  start  the  chart  recorder 
and  initiate  the  mass  spectrometer  scan),  but  also  provided  sufficient 
time  for  thermal  equilibrium  to  be  established. 
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The  amount  of  time  required  for  the  radial  temperature  profile 
in  the  reaction  tube  to  be  flat  within  0.1  K  may  be  estimated  from 
the  following  analysis.  Neglecting  natural  convection  and  heat  losses 
from  the  ends  of  the  reaction  tube,  the  radial  temperature  profile 
as  a  function  of  time  is  found  from  the  solution  of 


1_  tfT  ,  dn  ^  i  tfT 
a  dt  ~  dr?  7  dr 


5-7 


where  a  is  the  thermal  diffusivity  of  the  gas  in  the  reaction  tube. 
Carslaw  and  Jaeger  [175]  have  solved  equation  5-7  subject  to  the 
boundary  condition  of  zero  initial  temperature,  constant  wall  temper¬ 
ature  and  radial  symmetry.  The  result  is 

5-8 
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In  this  equation,  T  is  the  reactor  wall  temperature,  R  is  the  reactor 
w 

tube  radius  and  the  eigenvalues,  a^,  are  the  roots  of  the  equation 


J„(a„R)  .  0 
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Since  the  reaction  tube  contained  at  least  90%  H2  at  all  times,  the 

thermal  diffusivity  of  H2  was  used  to  evaluate  equation  5-8.  A  further 

assumption  inherent  in  equation  5-8  is  that  a  is  invariant  with  respect 

2 

to  temperature.  This  is  not  true  for  H2  as  a  goes  as  approximately  T 
[77].  However,  a  worse  case  calculation  can  be  performed  by  evaluating 
o  at  the  initial  temperature  (300  K)  of  the  gas.  The  average  and  cen¬ 
terline  temperatures  in  the  reaction  tube  are  shown  in  Figure  5-3. 

Five  eigenvalues  were  used  in  the  evaluation  of  equation  5-8  to  achieve 
this  result.  The  reactor  centerline  temperature  was  found  to  be  within 
0.01%  of  the  wall  temperature  (e.g.  0.1  K  at  T^^  =  1000  K)  for  times 
greater  than  7  seconds.  Evaluating  a  at  higher  temperatures  decreased 
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Figure  5-3 

Reaction  Tube  Temperature  During  Heating 
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this  time  sharply  (e.g.  evaluating  a  at  1000  K  required  only  1.9  s  to 
achieve  the  same  results).  The  influence  exerted  on  the  thermal  dif- 
fusivity  by  the  presence  of  VH^  in  the  system  is  on  the  order  of  10% 
to  15%.  This  is  considered  insignificant  relative  to  the  choice  of 
an  appropriate  temperature  for  evaluating  a.  Thus,  based  on  the  re¬ 
sults  depicted  in  Figure  5-3,  it  is  concluded  that  the  15  s  delay, 
between  charging  the  reactor  and  the  initiation  of  data  collection, 
was  sufficient  to  allow  the  gas  in  the  tube  to  reach  thermal  equilibrium. 

Mass  spectrometer  scans  were  made  every  10,  30,  60  or  120  seconds, 
depending  or.  the  reaction  rate,  and  data  was  taken  either  for  1200  se¬ 
conds  or  until  the  signal  was  too  small  to  be  reliably  measured.  Back¬ 
ground  signal  measurements  were  made  before  and  after  the  decomposition 
data  were  acquired  at  each  temperature  in  order  to  determine  whether  or 
not  the  VHj  background  signals  were  increasing  with  system  exposure  to 
the  VH^  species.  A  slight  increase  in  PH2  and  A.sH^  background  was  no¬ 
ticed  over  the  course  of  the  experiments,  but  the  change  at  any  one 
data  gathering  session  was  insignificant.  Changes  in  the  NH^  background 
could  not  be  observed  due  to  the  large  H2O  background  signals  present. 


rHAPTPO 

RESULTS  AND  DISCUSSION 
Chemical  Eoui librium  Investigation 


Introduction 

The  product  of  the  equilibrium  calculations  was  the  composition 
of  the  vapor  phase  in  the  presence  of  excess  condensed  phases.  The 
composition  was  investigated  as  a  function  of  temperature,  pressure 
and  inlet  gas  composition.  The  usual  procedure  was  to  vary  one  of 
the  operating  conditions  while  holding  the  remaining  ones  at  their 
base  values.  The  graphical  representation  of  these  results  illustrates 
the  equilibrium  vapor  mole  fraction  of  each  of  the  species  versus  the 
parameter  varied.  Since  the  primary  objective  of  this  study  was  to 
examine  unintentional  Si  incorporation  levels,  mole  fractions  are 
shown  typically  down  to  a  level  of  10  (0.1  ppb).  The  Si  species 

were  always  found  to  be  below  10  ^  mole  fraction.  Therefore,  only 
the  lower  five  orders  of  magnitude  were  shown  in  many  cases  unless  the 
upper  range  was  necessary  to  understand  the  results.  A  full  parametric 
analysis  was  performed  and  over  160  plots  were  generated.  In  many 
cases,  the  results  were  similar  to  analogous  system,  thus,  this  chap¬ 
ter  includes  only  those  graphs  necessary  for  understanding  the  prin¬ 
cipal  phenomena  predicted.  In  interpreting  these  plots,  it  should  be 
realized  that  an  excess  specie  serves  to  hold  the  activity  of  that 
specie  at  a  constant  value.  For  example,  with  solid  Si02  present,  the 
activity  of  Si02  is  fixed  at  unity  and  therefore  the  product  of  the  Si 


i 
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and  O2  partial  pressures  is  also  fixed.  Thus«  changes  in  operating 
parameters  that  alter  the  Og  fugacity  will  alter  the  Si  activity  by 
the  same  degree  in  an  inverse  fashion. 

The  GaAs  Chloride  System 

The  effect  of  temperature  on  the  species  present  in  the  GaAs 
chloride  system  source  zone  (100  kPa  pressure,  inlet  composition: 

IS  AsClj  in  Hg)  is  shown  in  Figures  6-1  through  6-4.  Figures  6-1  and 
6-2  apply  to  the  system  which  used  a  liquid  group  III  source,  Ga^As-j  j^ 
and  Figures  6-3  and  6-4  represent  the  results  for  the  system  which  em¬ 
ployed  GaAs(c)  as  the  group  III  source  material.  At  low  temperatures, 
GaClg  and  GaClj  became  relatively  important  gallium  vapor  species  along 
with  GaCl  in  the  solid  source  system.  In  the  liquid  source  system, 

GaCI  is  the  dominant  gallium  species  over  the  entire  temperature  range 
examined  (873  <T  <1173  K).  In  both  systems,  As^  was  the  dominant 
arsenic  species  at  low  temperature  while  AS2  became  important  at  high 
temperatures-  In  contrast  to  previous  studies  [20,  21],  the  trimer, 

As^,  mole  fraction  was  not  negligible.  In  general,  comparison  of  the 
silicon  activity  for  the  two  source  zones  revealed  that  the  silicon 
activity  associated  with  the  6aAs(c)  source  material  was  much  lower 
than  that  which  resulted  when  a  liquid  source  material  was  employed. 

The  predominant  silicon  species  in  the  vapor  phase  of  the  system  which 
used  a  solid  source  were  the  higher  silicon  chlorides  in  contrast  to 
the  hydrogen  rich  silicon  species  found  in  the  system  which  used  a  li¬ 
quid  source.  An  additional  interesting  feature  is  that  the  total  mole 
fraction  of  silicon  compounds  in  the  vapor  for  the  system  which  employed 
a  solid  source  was  greater  than  that  for  the  system  which  employed  a 
liquid  source.  At  first  glance,  this  fact  seems  contradictory  to  the 
lower  observed  silicon  activity. 
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Figure  6-1 

Effect  of  Temperature  in  the  GaAs  Chloride  System  Source 
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The  following  reaction  equations  may  be  written  to  describe  the 


formation  of  silicon  chlorides,  chi orosi lanes  and  silane  resulting 
from  reactions  with  the  quartz  reactor  wall. 

Si02(c)  +  nHCl  +  (4-n)H2  *  2H2O  +  SiH^_^Clj^  n=0,l,2,3.4  6-1 

kSi02(c)  +  krnHCl  +  k(2-m/2)H2  »  2kH20  +  51,^01^1^^ 

k*l,2  m»0. 1,2.3  6-2 

Reactions  6-1  and  5-2  represent  a  set  of  independent  formation  reactions 
which  describe  the  interplay  between  the  dominant  vapor  phase  silicon 
species  present  in  the  system.  Assuming  ideal  gas  behavior,  the  equili¬ 
brium  constants  for  these  reactions  are  as  follows: 

n,n'  ^H-0  ^SiH.  „C1  '  ^HCl 
c  4-n  n  i 


and 


6-4 


2  ^HCl 

where  y^=  vapor  phase  mole  fraction  of  species  i 

P.j.=  system  total  pressure  ratio  (total  pressure/reference 

state  pressure) 

The  activity  of  silicon  residing  in  a  condensed  phase  which  is  at 
equilibrium  with  the  vapor  phase  may  be  calculated  from  any  reaction 
using  a  vapor  Si  species  reactant  and  solid  Si  product.  For  example, 
consider  the  following  reactions  and  subsequent  equilibrium  expressions 
for  the  activity  of  Si(c) 


Si(c)  +  2H2  =  SiH^ 


Si(c)  +  4HC1  »  SiCl^  +  2H2 


*Si  - ^ 

Kc  P 


6-5 


‘Si' 


-  y. 


5  T 

y2 

SiCU 


^6  ^HCl  ^T 
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Here  K.  is  the  equilibrium  constant  for  reaction  (i«5,6).  Other  equi¬ 
valent  relations  may  be  written  in  order  to  calculate  the  condensed 
phase  silicon  activity  but  the  models  suggested  in  equations  6-5  and 
6-6  serve  as  convenient  points  of  focus  since  either  SiH^  or  SiCl^  is 
usually  a  significant  silicon  vapor  species  in  the  systems  studied. 

In  particular,  for  those  systems  using  Hg  as  the  carrier  gas,  the  mole 
fraction  of  H2  is  nearly  constant  with  a  value  close  to  unity.  There¬ 
fore,  a^^.  will  track  the  SiH^  mole  fraction  and  is  inversely  propor¬ 
tional  to  the  system  pressure.  Both  models  of  course  yield  identical 
values  for  the  silicon  activity  when  applied  to  the  same  situation.  The 
activity  of  Si  presented  in  these  plots  can  be  viewed  as  the  value  found 
in  a  solid  phase  in  equilibrium  with  a  vapor  having  the  composition 
shown.  In  order  to  translate  this  into  a  solubility,  the  nature  of  the 
solid  phase  must  be  considered  (i.e.  the  activity  coefficient  must  be 
known).  Due  to  the  low  degree  of  doping  encountered  in  these  systems, 
(e.g.  typically  <  10^^  cm"^,  which  is  <100  ppb),  the  activity  coeffi¬ 
cient  can  be  represented  by  Henry's  constant,  which  is  invariant  with 
respect  to  composition.  Therefore,  an  increase  in  Si  activity  corre¬ 
sponds  to  an  increase  in  solubility.  Thus,  the  models  provided  by 
equations  6-5  and  6-6  may  be  used  to  predict  the  direction  of  change 
in  the  silicon  concentration  based  on  the  knowledge  of  the  vapor  phase 
equilibrium  composition. 

The  lower  silicon  activity  associated  with  the  solid  GaAs  source 
can  therefore  be  viewed  as  due  to  a  suppressed  SiH^  concentration  when 
compared  to  the  liquid  source  (equation  6-5).  In  the  source  zone, 
which  employed  solid  GaAs,  the  presence  of  primarily  higher  chlorides 
and  chlorosilanes  at  the  lower  source  zone  temperatures  was  a  result 


of  less  gallium  being  present  in  the  vapor  phase  than  was  present  when 
a  liquid  source  was  employed.  Since  6a  was  in  excess  in  both  the  li¬ 
quid  and  solid  source  systems,  the  activity  of  Ga  was  constrained  with 
the  liquid  source  exhibiting  a  higher  Ga  activity.  Thus,  sufficient 
HCl  was  available  due  to  the  decomposition  of  AsClj  and  the  absence  of 
enough  Ga  to  consume  it  to  enhance  reactions  6-1  and  6-2  for  large  k 
and  ro  values.  Figure  6-5,  which  shows  the  chloride  system  pre-source 
zone  (IS  AsClj  in  H2.  no  group  III  source  material  present),  further 
supports  this  analysis.  The  absence  of  group  III  chlorides  caused  the 
total  amount  of  silicon  in  the  vapor  to  increase  above  the  level  ob¬ 
served  in  the  solid  source  system  while  the  condensed  phase  silicon 
activity  decreased  even  further. 

Table  6-1  lists  the  enthalpy  of  formation  at  298  K  and  the  Gibbs 
energy  of  reaction  at  973  K  for  some  of  the  vapor  species  described  by 
reactions  6-1  and  6-2.  The  large  negative  enthalpies  of  formation  are 
indicative  of  strong  interatomic  bonds  and  therefore  stable  species. 

Since  equilibrium  represents  the  state  having  the  lowest  value  for  the 
Gibbs  energy  of  the  system,  species  with  a  lower  Gibbs  energy  of  re¬ 
action  are  favored.  Therefore,  providing  sufficient  chlorine  to  react 
with  the  silicon  species  results  in  a  higher  total  silicon  concentration 
in  the  vapor  phase  but,  due  to  the  stability  of  these  species,  a  lower 
activity  of  solid  silicon  in  the  condensed  phase.  The  relative  stability 
of  silicon  halides,  when  compared  to  silicon  hydrides  was  also  recog¬ 
nized  by  Rai-Choudhury  [12]. 

The  outlet  equilibrium  compositions  of  the  source  zones  using 
solid  or  liquid  source  materials  at  973  K  were  used  as  input  to  the 
mixing  zone  and  the  effect  of  mixing  zone  temperature  was  investigated. 


Mole  Fraction 


*Reference  state;  Si(c),  Clo(v),  H-(v)  at  298  K  and  100  kPa 
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The  system  using  a  liquid  source  material  displayed  behavior  which 
was  nearly  identical  to  the  source  zone  behavior.  This  supports 
Weiner's  model  [30]  which  suggested  using  the  outlet  composition  of 
the  source  zone  as  the  inlet  composition  to  the  deposition  zone.  This 
model  is  only  applicable,  of  course,  if  the  source  and  mixing  zones 
are  operated  at  the  same  temperature.  Justification  for  isothermal 
operation  of  the  source  and  mixing  zones  comes  from  noting  that  the 
condensed  phase  silicon  activity  increases  with  temperature.  There¬ 
fore,  it  is  desirable  to  operate  the  mixing  zone  at  the  source  zone 
temperature  in  order  to  minimize  the  silicon  activity  and  at  the  same 
time  prevent  deposition  of  GaAs  in  this  zone. 

The  behavior  of  the  mixing  zone,  when  fed  from  a  source  zone  using 
solid  GaAs  as  the  source  material,  differed  from  that  of  the  source 
zone  alone  in  that  the  mole  fractions  of  all  of  the  chlorinated  vapor 
phase  silicon  species  increased  with  temperature,  as  did  the  condensed 
phase  silicon  activity.  Although  the  silicon  activity  in  the  mixing 
zone  increased  more  slowly  with  temperature  than  it  did  in  the  source 
zone,  the  results  of  this  analysis  suggest  operating  the  mixing  zone 
at  a  temperature  equal  to  or  less  than  that  of  the  source  zone,  in  or¬ 
der  to  attain  the  lowest  possible  silicon  activity. 

The  effect  of  temperature  on  the  deposition  zone,  which  was  fed 
from  source  and  mixing  zones  using  temperatures  of  973  K,  pressures  of 
100  kPa,  a  liquid  source  material  and  1%  AsCl^  in  H2  as  inlet  into  the 
reactor,  has  shown  GaCl  and  As^  to  be  the  dominant  group  III  and  V 
species.  A  measure  of  the  supersaturation  of  the  vapor  was  defined 


based  on  the  reaction 

GaAs(c)  =  Ga(g)  +  As(g) 


6-7 
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Using  the  equilibrium  relationship, 
fined  as 


P  .  ^'as 

^at 


the  saturation  ratio,  Is  de- 


6-8 


where  =  exp(^G;:^yRT)  6-9 

This  saturation  ratio  v/as  observed  to  decrease  from  4x10^  at  773  K  to 
700  at  973  K,  which  shows  that  the  vapor  phase  was  highly  supersaturated. 
When  S102(c)  was  not  included  In  the  deposition  zone,  the  condensed 
phase  silicon  activity  changed  only  slightly  from  9x10“^  at  773  K  to 
6x10“^  at  973  K.  Including  Si02(c)  in  the  deposition  zone  calculation 
resulted  in  the  condensed  phase  silicon  activity  becoming  a  strong 
function  of  temperature  (due  to  reversal  of  reactions  6-1  and  6-2)  with 
the  activity. value  at  973  K  remaining  unchanged  and  the  773  K  value 
fall ing  to  SxlO'^. 

When  solid  GaAs  was  used  as  the  group  III  source  material,  the 
following  results  were  obtained  for  the  deposition  zone.  The  saturation 
ratio  fell  from  a  value  of  200  at  773  K  to  the  expected  value  of  1  at 
973  K,  thus  revealing  the  system  to  be  much  less  supersaturated  than 
the  liquid  source  material  counterpart.  This  lower  degree  of  super¬ 
saturation  was  due  to  much  less  GaCl  being  present  in  the  vapor.  The 
dominant  group  V  specie  was  As^,  but  GaCl,  GaCl2  and  GaCl^  were  all  im¬ 
portant  contributors  to  the  group  III  vapor  species.  The  condensed 
phase  silicon  activity  was  found  to  increase  with  temperature  from 
9x10"®  at  773  K  to  7x10"^  at  973  K  for  the  case  where  SiOgCc)  was  not 
included  in  the  deposition  zone.  When  Si02(c)  was  included,  the  sili¬ 
con  activity  at  773  K  fell  to  2x10"^. 
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The  effect  of  pressure  was  investigated  over  the  range  of  1  to 
1000  kPa  (temperature  973  K,  inlet  composition:  IS  AsCl^  in  H2)  for 
both  the  solid  and  liquid  group  III  source  materials.  Over  the  en¬ 
tire  range  studied  GaCl  was  the  dominant  group  III  vapor  specie  in 
the  system  using  a  liquid  source  material,  v/hile  for  the  system  using 
solid  GaAs  as  the  source  material.  GaCl  became  the  dominant  group  III 
specie  at  pressures  below  100  kPa,  but  competed  with  GaClg  and  GaCl^ 
at  the  higher  pressures.  The  dominant  group  V  vapor  specie  was  As^ 
at  pressures  above  10  kPa  with  AS2  becoming  important  below  this  pres¬ 
sure,  in  agreement  with  Gentner  et  al.  [21],  in  both  the  liouid  and  solid 
source  systems. 

Figures  6-6  and  6-7  show  the  lower  five  orders  of  magnitude  in 
mole  fraction  and  the  condensed  phase  silicon  activities  in  the  liquid 
and  solid  material  source  zones.  The  silicon  activity  in  the  system, 
which  used  a  liquid  source  material,  reached  a  maximum  at  a  pressure 
of  4  kPa  and  then  decreased  with  increasing  pressure.  This  behavior 
has  not  previously  been  reported  in  the  literature  presumably  due  to 
the  constrained  nature  of  previous  equilibrium  calculations.  Referring 
to  equation  6-5  reveals  that  the  maximum  in  the  silicon  activity  was 
due  to  the  SiH^  mole  fraction  rising  faster  than  Py  in  the  1  to  10  kPa 
range.  Applying  equation  6-3  to  the  specie  SiH^  (n=0)  and  referring 
to  Figure  6-6  shows  that  the  SiH^  mole  fraction  dependence  on  pressure 
deviates  from  linear  behavior  due  to  the  H2O  mole  fraction  changes  in 
this  range  =1).  The  change  in  H2O  mole  fraction  was  due  to 
changes  in  the  total  amount  of  Si02(c)  which  reacted  with  the  vapor. 
Reaction  6-1  was  important  in  this  system  and,  as  the  pressure  in¬ 
creased,  caused  more  Si02(c)  to  react  which  generated  more  H2O. 
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The  dominant  silicon  vapor  species  present  at  the  low  end  of  the 
pressure  range  was  SiO  which  was  formed  via  the  reaction 

2Si02(c)  +  Hg  =  2SiO  +  H2O  6-10 

with  the  corresponding  equilibrium  relationship 
'<3  '  ^SiO  ’’t  /  \ 

Thus,  the  observed  minimum  in  the  H2O  mole  fraction  was  due  to  the 
interaction  between  the  decreasing  mole  fraction  of  SiO  with  the  in¬ 
creasing  system  pressure  from  reaction  6-7  along  with  H2O  generation 
from  reaction  6-1.  Reaction  6-2  was  not  important  in  this  situation. 

The  source  zone  which  used  solid  6aAs  as  the  group  III  source 

/ 

material  showed  a  strictly  decreasing  condensed  phase  silicon  activity 
with  increasing  system  pressure  and  is  best  described  via  reactions.  6-1 
and  6-2  in  conjunction  with  the  silicon  activity  model  provided  by  equa¬ 
tion  6-6.  The  decrease  in  silicon  activity  was  due  to  the  3  order  of 
magnitude  increase  in  SiCl^  mole  fraction  being  offset  by  an  order  of 
magnitude  increase  in  HCl  mole  fraction  )  and  the  Py^  depen¬ 

dence  of  the  silicon  activity.  Once  again,  the  activity  of  silicon  in 
the  system  using  solid  GaAs  as  a  source  material  was  much  less  than 
the  activity  resulting  from  using  liquid  Ga^As.j_^  as  the  group  III 
source  material. 

The  mixing  and  deposition  zones  were  studied  using  only  the  liquid 
group  III  source  material  in  the  source  zone.  The  mixing  zone  results 
were  again  essentially  the  same  as  those  of  the  source  zone  (operated 
at  the  same  temperature)  and  therefore  do  not  require  further  discussion 
The  deposition  zone,  shown  in  Figures  6-8  and  6-9,  exhibited  a 
saturation  ratio  of  approximately  2x10^  at  a  pressure  of  1  kPa  rising 
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to  approx i ir.a tel y  2x10  at  1000  IcPa.  Thus,  the  deposition  zone  was 
supersaturated  over  the  entire  pressure  range  investigated. 

The  maximum  value  of  the  condensed  phase  silicon  activity  ob¬ 
served  for  the  source  and  mixing  zones  was  not  prevalent  in  the  de¬ 
position  zone  since  Si02(c}  was  not  included  in  the  deposition  zone 
model.  Therefore,  the  silicon  activity  decreased  with  increasing 
pressure  in  accord  with  equation  6-S. 

The  mole  fraction  of  AsCl^  present  in  the  feed  stream  was  varied 
from  O.IS  to  lOS  in  order  to  determine  its  effect  on  the  condensed 
phase  silicon  activity.  For  the  source  zone  utilizing  a  liquid  group 
III  source  material  most  of  the  chlorine  atoms  provided  by  the  decom¬ 
position  of  AsCl^  were  used  to  generate  GaCl.  Therefore,  the  condensed 
phase  silicon  activity  was  not  appreciably  affected  until  large  concen¬ 
trations  of  AsClj  were  reached.  The  silicon  activity  was  found  to  de¬ 
crease  from  7x10“^  at  0.1*  AsCl^  to  6x10“^  at  15  AsCl^  and  finally  to 
7x10"^  at  105  AsClj.  These  results  agree  qualitatively  with  previous 
calculations  [11,  12]  and  observations  [8,  9,  10,  13]. 

A  very  marked  effect  on  the  silicon  activity  was  observed  as  the 

ASCI2  inlet  concentration  was  varied  in  the  system  using  solid  GaAs  as 

the  group  III  source  material.  Since  less  GaCl  was  generated  in  this 

system  when  compared  to  the  system  using  a  liquid  group  III  source, 

more  HCl  was  available  from  the  decomposition  of  AsCl^  to  react  with  the 

silicon  vapor  species.  Thus,  the  condensed  phase  silicon  activity  was 
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found  to  decrease  uniformly  from  9x10”  at  0.15  AsCl^  to  1.5x10"  at 
105  AsCl^  inlet  concentration. 

The  reason  less  GaCl  was  generated  in  the  system  using  solid  GaAs 
as  the  group  III  source  is  explained  by  the  following  reactions. 
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i  Ga/s^^^(l)  +  HCl(g)  =  GaCl{g)  +  ^  H2(g)  +  ^  As^Cg) 

Reaction  6-12  refers  to  the  system  using  the  solid  source  and,  at  973  K, 
has  a  Gibbs  energy  change  of  3.2  kcal/mole  while  reaction  6-13,  repre¬ 
senting  the  system  with  a  liquid  source,  undergoes  a  Gibbs  energy  change 
of  -12.5  kcal/mole.  The  negative  Gibbs  energy  change  of  reaction  6-13 
causes  the  products  of  the  reaction  to  be  favored. 

The  influence  of  temperature  and  pressure  on  the  vapor  phase  III/V 
ratio  for  solid  and  liquid  group  III  sources  is  shown  in  Figure  6-10. 
When  a  liquid  source  is  used,  the  III/V  ratio  is  limited  to  a  maximum 
of  3.  This  is  due  to  most  of  the  Ga  being  transported  as  GaCl  and  the 
3  to  1  ratio  of  Cl  atoms  to  As  atoms  in  AsCl^.  Increasing  temperatures 
cause  the  III/V  ratio  to  decrease  when  a  liquid  source  is  used  due  to 
the  increased  amount  of  As  present  in  the  liquid.  The  III/V  ratio  in¬ 
creases  with  increasing  temperature  for  solid  sources  due  to  a  lowering 
of  the  Gibbs  energy  (due  to  entropic  effects)  of  reaction  6-12.  In¬ 
creasing  pressure  causes  a  decrease  in  the  III/V  ratio  for  both  the 
solid  and  liquid  source  systems  due  to  reversal  of  reactions  6-12  and 
6-13. 

The  source  zone  results  were  carried  through  the  mixing  zone  and 
deposition  zone  for  the  system  using  a  liquid  group  III  source.  The 
mixing  zone  yielded  the  same  results  as  did  the  source  zone  (both  zones 
operating  at  a  temperature  of  973  K)  and  the  deposition  zone  showed  the 
same  trends  as  were  observed  in  the  source  and  mixing  zones  except  that 
the  values  of  silicon  activity  were  slightly  higher  due  to  the  lower 
deposition  zone  temperature  (873  K). 
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Figure  6-10 

Temperature  and  Pressure  Effects  on  the  Vapor  III/V  Ratio 
in  the  GaAs  Chloride  System  Source  Zones  (Solid  and  Liquid  Sources) 


A  very  effective  way  to  reduce  the  silicon  activity  in  systems 
using  a  liquid  group  III  source  is  to  add  HCl  (or  AsCl^)  downstream 
of  the  source  zone.  This  allows  the  chlorine  atoms  to  react  with  the 
silicon  species  instead  of  generating  additional  GaCl. 

Figure  6-11  shows  the  effect  of  adding  small  quantities  of  HCl 
to  the  mixing  zone  on  the  condensed  phase  silicon  activity  in  the 
deposition  zone  (basis:  1  mole  of  vapor  in  the  mixing  zone).  The 
initial  HCl  mole  fraction  in  the  deposition  zone  prior  to  the  addi¬ 
tion  of  any  HCl  was  4.5x10'^.  In  accord  with  reactions  6-1,  6-2,  6-5 
and  6-6,  the  silicon  compounds  shifted  from  being  hydrogen  rich  to 
chlorine  rich  and  the  silicon  activity  decreased  markedly.  The  addi¬ 
tion  of  AsCl^  had  the  same  effect  except  that  the  activity  decrease 
is  slightly  more  pronounced  since  there  are  three  chlorine  atoms  per 
molecule  of  AsCl^  compared  to  one  for  HCl. 

Another  method  of  decreasing  the  condensed  phase  silicon  activity 
is  to  add  H^O  to  the  system.  This  causes  a  decrease  in  the  total  amount 
of  silicon  in  the  vapor  by  shifting  reactions  6-1  and  6-2  in  favor  of 
Si02(c).  This  effect  is  demonstrated  in  Figure  6-12  for  small  addi¬ 
tions  of  H2O  to  the  mixing  zone  (basis:  1  mole  of  vapor  in  the  mixing 
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zone)  where  the  mole  fraction  of  H2O  prior  to  the  additions  was  5.5x10 
This  effect  was  predicted  by  Rai-Choudhury  [12]  and  observed  by  Palm 
et  al .  [17]. 

Replacing  the  hydrogen  carrier  gas  with  an  inert  gas  is  yet 
another  method  of  reducing  the  condensed  phase  silicon  activity.  This 
method  also  reduces  the  total  amount  of  silicon  in  the  vapor  by  shift¬ 
ing  reactions  6-1  and  6-2  in  favor  of  Si02Cc)  as  shown  in  Figure  6-13. 
Reduction  of  silicon  incorporation  using  an  inert  to  replace  H,  was 


nitial  HC1  before  additions 


studied  by  Seki  et  al .  [18]  and  observed  by  Ozeki  et  al .  [19].  The 
curvature  in  the  silicon  activity  is  best  understood  by  referring  to 
reaction  6-5.  As  the  hydrogen  carrier  gas  is  replaced  by  an  inert, 
the  mole  fractions  of  H2  and  SiH^  decrease.  The  competing  nature  of 
these  two  mole  fractions  causes  a  maximum  in  the  condensed  phase  sili¬ 
con  activity  to  occur  at  approximately  90%  inerts  after  which  the  si¬ 
lane  mole  fraction  rapidly  goes  to  zero  and  the  silicon  activity  de¬ 
creases  to  a  very  small  value.  The  silicon  activity  will  never  reach 
zero  as  predicted  by  6-6  since  reactions  6-1,  6-2,  6-5  and  6-6  are  not 
valid  models  in  systems  devoid  of  hydrogen.  Instead,  small  concen¬ 
trations  of  Si(g),  SiO(g),  etc.  will  remain  in  the  vapor  to  provide  a 
nonzero  but  very  small  condensed  phase  silicon  activity. 

The  use  of  solid  GaAs  as  the  group  III  source  material  appears  to 
offer  an  advantage  over  the  liquid  group  III  source  in  that  lower  con¬ 
densed,  phase  silicon  activities  were  predicted  by  these  thermodynamic 
models.  It  must  be  emphasized,  however,  that  the  solid  GaAs  source 
was  assumed  to  be  pure  (i.e.  devoid  of  Si  and  other  contaminants)  and, 
from  a  theniiodynamic  point  of  view,  that  the  purity  of  an  epitaxial 
layer  can  be  no  better  than  that  of  the  source  material  unless  methods 
are  employed  to  improve  the  purity  (e.g.  additions  of  HCl ,  H2O,  etc.) 
during  the  CVO  process.’ 

The  GaAs  Hydride  System 

The  effect  of  temperature  on  the  silicon  species  present  in  the 
group  III  and  group  V  source  zones  of  the  GaAs  hydride  system  is  shown 
in  Figures  6-14  and  6-15.  The  condensed  phase  silicon  activity  for  the 
group  V  source  zone  was  essentially  the  same  as  that  which  would  be 


observed  in  the  group  V  pre-source  zone  (AsH^  absent)  since  arsine  did 
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Effect  of  Temperature  on  the  Group  III  Source  Zone  of  the 

GaAs  Hydride  System 


not  compete  with  silicon  for  any  atoms  in  the  vapor  other  than  hydro¬ 
gen.  The  group  III  source  zone  showed  a  much  greater  silicon  activity 
than  the  group  III  pre-source  zone  (6a(l)  absent),  shown  in  Figure  6-16, 
since  the  liquid  gallium  source  reacted  strongly  with  HCl  to  form  GaCl . 
This  forced  the  silicon  species  to  be  rich  in  hydrogen  and  therefore 
the  condensed  phase  silicon  activity  was  larger.  As  can  be  seen  from 
Figures  6-14  and  6-15,  the  group  V  source  zone  was  primarily  responsible 
for  the  silicon  activity  at  low  temperatures  while  at  high  temperatures, 
the  group  III  source  zone  contribution  to  the  silicon  activity  also 
became  important.  As  was  observed  in  the  chloride  system  for  the  li¬ 
quid  Ga^As^  ^  source,  the  dominant  group  III  vapor  species  was  GaCl 
with  As^  being  the  dominant  group  V  vapor  species  at  temperatures  below 
1073  K  and  As2  dominant  above  this  temperature. 

The  results  of  these  group  III  and  V  source  zone  calculations  at 
973  K  were  then  combined  and  the  effect  of  temperature  in  the  mixing 
zone  was  investigated.  Figure  6-17  shows  this  effect  on  the  silicon 
activity  and  species  in  the  lower  five  orders  of  magnitude  in  mole 
fraction.  The  silicon  activity  and  silane  mole  fraction  were  found 
to  be  lower  than  the  values  for  the  group  V  source  zone  alone  at  low 
temperatures  due  to  the  dilution  effect  of  adding  the  two  source  zone 
streams  together  (equal  molar  flowrates  were  assumed  in  each  source 
zone).  Since  the  silane  mole  fraction,  and  therefore  the  silicon 
activity,  in  each  source  zone  was  approximately  the  same  at  1173  K, 
the  resulting  silicon  activity  in  the  mixing  zone  was  the  same  as  that 
at  the  outlet  of  either  source  zone. 

Using  the  results  from  the  mixing  zone  at  973  K,  the  deposition 
zone  v/as  studied  in  the  absence  of  Si02(c).  The  saturation  ratio,  as 
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Effect  of  Temperature  on  the  Mixing  Zone  of  the  GaAs  Hydride  System 


defined  by  equation  6-8,  was  found  to  decrease  from  2x10  at  773  K 

to  100  at  973  K,  indicating  that  the  vapor  was  supersaturated  over 

the  entire  temperature  range.  The  condensed  phase  silicon  activity 
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varied  only  slightly,  increasing  from  2x10  at  773  K  to  5x10  at 
973  K,  which  is  consistent  with  the  effects  observed  in  the  chloride 
system  using  a  liquid  6a^As^_^  group  III  source.  The  lesser  degree 
of  supersaturation  found  for  the  hydride  system,  when  compared  to  the 
chloride  system,  was  due  to  the  lesser  amounts  of  group  III  and  V 
species  present  in  the  vapor.  The  dilution  effect  of  adding  the  two 
source  zone  streams  together  causes  the  group  V  vapor  species  to  be 
one-half  the  mole  fraction  observed  in  the^chloride  system.  The  use 
of  HCl  coupled  with  this  dilution  effect  reduced  the  amount  of  group 
III  species  transport  to  one-sixth  of  that  in  the  chloride  system. 

The  group  III  and  V  source  zones  were  investigated  as  a  function 
of  pressure  at  a  temperature  of  973  K.  The  silicon  activity  in  the 
group  V  source  zone  remained  constant  to  10  kPa  then  fell  from  a  value 
of  3x10'^  to  4x10'^  at  1000  kPa.  The  group  III  source  zone  exhibited 
a  maximum  in  the  condensed  phase  silicon  activity  at  a  pressure  of  4 
kPa,  as  did  the  chloride  system  source  and  mixing  zones  using  a  liquid 
group  III  source.  Upon  combining  the  two  hydride  system  source  zones 
and  performing  the  mixing  zone  equilibrium  calculation,  the  silicon 
activity  in  the  mixing  zone  became  a  decreasing  function  of  pressure. 

The  effect  of  pressure  on  the  deposition  zone  of  the  hydride 
system  very  closely  matched  that  of  the  chloride  system.  This  result 
was  expected  since  the  source  zones  of  the  two  systems  are  the  only 
differences  between  the  two  and  once  downstream  of  the  source  zones, 
the  equilibrium  chemistry  of  the  hydride  and  chloride  systems  are  the 
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The  concentration  of  AsH^  in  the  feed  gas  stream  of  the  group 
V  source  zone  was  found  to  have  no  effect  on  the  condensed  phase  sili¬ 
con  activity.  This  was  due  to  the  silicon  species  being  rich  in  hydro¬ 
gen  and,  in  the  presence  of  a  large  fraction  of  hydrogen  carrier  gas 
(>  90«),  the  hydrogen  atoms  released  from  the  decomposition  of  AsH^  did 
not  contribute  significantly  to  the  overall  system  hydrogen  content. 
Therefore,  the  silane  mole  fraction  was  not  significantly  affected.  In 
contrast,  Pogge  and  Kemlage  [37]  found  that  increased  AsH^  concentra¬ 
tions  decreased  the  free  carrier  concentrations  in  epitaxial  GaAs. 

They  cited  kinetic  effects,  however,  not  thermodynamic  limitations  as 
the  reason  for  their  observations. 

Increasing  the  concentration  of  HCl  in  the'  feed  stream  to  the 
group  III  source  considerably  decreased  the  condensed  phase  silicon 
activity  by  forming  chlorine  rich,  as  opposed  to  hydrogen  rich,  silicon 
species.  Thus,  one  method  of  decreasing  the  silicon  activity  while 
maintaining  a  constant  vapor  III/V  ratio  is  to  increase  both  the  HCl 
and  AsH^  mole  fractions  in  the  feed  streams  to  each  source  zone  to¬ 
gether.  This  effect  is  demonstrated  for  the  deposition  zone  in 
Figure  6-18.  Additionally,  if  various  III/V  ratios  are  desired,  it 
is  advisable  to  operate  the  system  with  a  large  HCl  concentration  in 
order  to  realize  a  low  silicon  activity  and  vary  the  AsH^  concentra¬ 
tion  until  the  appropriate  III/V  ratio  is  attained. 

As  was  discussed  for  the  chloride  system,  a  very  effective  method 
for  reducing  the  silicon  activity  is  to  make  small  additions  of  HCl  or 
H2O  to  the  mixing  zone.  This  preserves  the  system  III/V  ratio  from 
the  hydride  system  source  zones  and  shifts  reactions  6-1  and  6-2  in 
favor  of  Si02(c).  The  results  for  the  hydride  system  were  essentially 
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the  same  as  those  obtained  for  the  chloride  system  and  will  not  be 
discussed  further. 

Replacement  of  the  hydrogen  carrier  gas  by  an  inert  gas  was  less 
effective  in  the  hydride  system  than  it  was  in  the  chloride  system 
since  hydrogen  was  provided  by  the  decomposition  of  HCl  and  AsH^. 

Figure  6-19  shows  this  effect  for  the  hydride  system  deposition  zone. 
Although  a  very  sharp  bend  is  observed  in  the  silicon  activity  even 
when  all  of  the  hydrogen  was  replaced  by  inerts,  the  condensed  phase 
silicon  activity  was  approximately  4x10“^.  Therefore,  replacing  the 
hydrogen  carrier  gas  by  an  inert  gas  was  not  an  acceptable  method  to 
achieve  low  silicon  activities.  The  addition  of  small  amounts  of  HCl 
or  H2O  to  the  mixing  zone  appears  to  be  the  most  promising  method  of 
attaining  very  low  silicon  activities  in  the  hydride  system. 

The  InP  Chloride  System 

The  InP  chloride  system  was  investigated  under  the  same  conditions 
as  the  GaAs  chloride  system  and  was  found  to  behave  similarly  in  many 
respects.  The  source  zone  using  liquid  as  the  group  III  source 

material  was  quite  similar  to  the  GaAs  system  in  that  InCl  and  were 
the  dominant  group  III  and  V  vapor  species.  One  difference  observed 
was  that  P2  was  not  as  significant  in  the  InP  system  as  AS2  was  in  the 
GaAs  system.  This  may  be  due,  however,  to  the  equilibrium  calculation 
being  constrained  as  a  result  of  the  lack  of  a  thermodynamic  data  set 
for  the  species  Pj.  The  condensed  phase  silicon  activities  for  the 
GaAs  and  InP  systems  were  found  to  be  essentially  the  same.  These  same 
comments  also  apply  to  the  mixing  zone  which  was  fed  from  this  source 
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Consistent  with  the  calculational  technique  employed  for  the 
GaAs  chloride  system,  the  vapor  phase  inlet  composition  to  the  depo¬ 
sition  zone  was  determined  by  the  equilibrium  composition  of  the  mixing 
zone.  The  mixing  zone  inlet  composition  was  determined  by  the  equilibrium 
source  zone  composition.  The  deposition  zone  of  the  InP  chloride  system 
was  studied  for  a  source  zone  using  a  liquid  group  II!  source 
with  a  vapor  inlet  composition  of  1«  PCl^  in  H2.  The  source  and  mixing 
zone  temperatures  were  both  973  K  and  the  pressures  in  the  source, 
mixing  and  deposition  zones  were  100  kPa. 

The  deposition  zone,  shown  in  Figures  6-20  and  6-21  for  the  sit¬ 
uation  where  Si02(c)  was  not  included  in*the  calculation,  demonstrates 
several  differences  between  the  GaAs  and  InP  chloride  systems.  First, 
the  saturation  ratio  for  the  InP  systan,  defined  in  analogy  with  equa¬ 
tions  6-7  and  6-8,  was  much  smaller  than  that  of  the  GaAs  system.  This 
is  due  to  the  Gibbs  energy  change  for  reaction  6-7  at  873K  being  95.2 
kcal/mole  for  GaAs  while  the  Gibbs  energy  change  for  the  analogous  InP 
system  is  only  69.5  kcal/mole.  The  value  of  the  equilibrium  constant 
for  the  InP  system  is,  therefore,  much  larger  than  that  for  the  GaAs 
system.  Thus,  larger  vapor  phase  In  and  P  partial  pressures  must  be 
present  in  order  to  attain  equal  degrees  of  supersaturation  in  the  two 
systems.  Since  the  saturation  ratio  definition  is  independent  of  the 
process  to  which  it  is  applied,  it  is  expected  that  the  InP  hydride  sys¬ 
tem  will  also  show  a  lesser  degree  of  supersaturation  than  the  GaAs  hydride 
system.  These  results  are  consistent  with  those  of  Shaw  [23]  who  also 
found,  through  thermodynamic  calculations,  that  InP  systems  are  less  super¬ 
saturated  than  GaAs  systems. 
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The  condensed  phase  silicon  activity  for  the  InP  system  was 
slightly  lower  than  that  of  the  GaAs  system,  due  to  the  lower  SiH^ 
mole  fraction,  but  increased  slightly  with  temperature.  Recall  that 
the  silicon  activity  in  the  GaAs  system  displayed  a  slight  decrease 
in  temperature  when  Si02(c)  was  not  included  in  the  calculation. 

The  InP  source  zone  using  solid  InP  as  the  group  III  source 
material  is  shown  in  Figures  6-22  and  6-23.  The  behavior  of  this 
system  was  markedly  different  than  that  of  the  analogous  GaAs  system 
in  that  InCl  was  clearly  the  dominant  group  III  vapor  specie  due  to  the 
favorable  Gibbs  energy  change  (-8  kcal/mole  at  973  K)  of  the  following 
reaction  when  compared  to  reaction  6-12  (3.2  kcal/mole  at  973  K). 

InP(c)  +  HCl (g)  =  InCl (g)  +  ^  P^Cg)  +  \  H2(g)  6-14 

The  formation  of  a  larger  amount  of  group  III  monochloride  consumed  much 
of  the  HCl  initially  present  in  the  system  and  caused  the  vapor  phase 
silicon  species  to  be  higher  in  hydrogen  content  than  those  species  in 
the  analogous  GaAs  system.  As  a  result,  the  condensed  phase  silicon 
activity  in  the  InP  chloride  system,  using  solid  InP  in  the  source  zone, 
was  found  to  be  much  larger  than  that  of  the  GaAs  system.  These  same 
comments  also  apply  to  the  mixing  and  deposition  zones  of  the  InP  sys¬ 
tem  which  follow  the  source  zone  using  solid  InP  as  the  group  III  source 
material . 

The  effect  of  pressure  on  the  InP  chloride  system  followed  closely 
that  observed  for  the  GaAs  system,  relative  to  the  previous  discussed 
differences.  Thus,  for  the  system  employing  liquid  In^^P^.^  as  the  group 
III  source,  the  condensed  phase  silicon  activity  in  each  zone  was 
slightly  lower  than  that  of  the  GaAs  system  and  the  InP  saturation  ratio 
in  the  deposition  zone  was  much  less  than  that  of  the  GaAs  system.  As 
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in  the  GaAs  system,  maxima  were  observed  In  the  condensed  phase 
silicon  activities  at  4  kPa  in  the  source  and  mixing  zones.  The 
chloride  system,  using  solid  InP  as  the  group  III  source,  also  dis¬ 
played  a  pressure  dependence  which  was  similar  to  that  of  the  anal¬ 
ogous  GaAs  system.  The  major  differences  being  that  InCl  was  clearly 
the  dominant  group  III  species  and  the  condensed  phase  silicon  acti¬ 
vity  was  much  higher  in  the  InP  system,  due  to  reasons  previously 
discussed  for  GaAs. 

The  influence  of  PCl^  inlet  concentration  on  the  InP  system  was 
similar  to  that  of  AsCl^  in  the  GaAs  system  with  the  following  differ¬ 
ences.  When  liquid  In^P^  was  used  as  the  group  III  source  material, 
the  condensed  phase  silicon  activity  was  always  20*  to  SOS  less  than 
that  of  the  corresponding  GaAs  systan.  This  was  due  to  InCl^  being 
present  in  a  smaller  amount  in  the  InP  system  than  GaClj  was  in  the 
GaAs  system.  Thus,  more  chlorine  was  available  to  react  with  the  sili¬ 
con  vapor  species.  When  solid  InP  was  used  as  the  group  III  source 
material,  the  silicon  activity  in  the  InP  system  was  greater  than 
that  of  the  corresponding  GaAs  system  due  to  the  large  amount  of  InCl 
formed. 

The  replacement  of  the  hydrogen  carrier  gas  by  an  inert  gas  was 
also  investigated  in  the  InP  system  with  the  results  following  the 
same  trends  as  did  the  analogous  GaAs  system.  The  differences  be¬ 
tween  the  systems  were  consistent  with  the  previous  discussions,  i.e. 
the  condensed  phase  silicon  activity  for  the  InP  system  using  solid 
InP  for  a  source  was  larger  than  the  GaAs  system  and  the  liquid  source 
InP  system  showed  a  slightly  reduced  silicon  activity  in  the  deposition 
zone  relative  to  the  GaAs  system. 
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Additions  of  HC1 ,  H2O  and  PClj  to  the  mixing  zone  of  the  InP 
system  were  also  studied.  The  trends  observed  were  the  same  as  those 
in  the  GaAs  system. 

The  InP  Hydride  System 

The  results  of  a  parametric  analysis  of  the  InP  hydride  system 
were  similar  to  those  of  the  GaAs  hydride  system  discussed  previously. 
The  condensed  phase  silicon  activity  in  the  InP  system  was  found  to 
be  consistently  less  than  that  in  the  GaAs  system  under  all  analogous 
conditions.  At  very  low  inlet  HCl  concentrations  (~0.1%)  the  silicon 
activities  were  nearly  the  same  while  at  large  inlet  concentrations 
('10%),  the  InP  system  exhibited  silicon  activities  which  were  half 
the  value  of  those  in  the  GaAs  system.  This  was  primarily  due  to  the 
more  negative  Gibbs  energy  of  formation  of  GaCl  (-41.5  kcal/mole  at 
973  K) ,  which  caused  a  greater  production  of  GaCl  via  reaction  6-13, 
than  InCl  in  the  analogous  reaction  where  In  is  the  group  III  species. 
The  differences  in  the  Gibbs  energy  of  formation  at  973  K  for  Ga(l)  and 
In(l )  were  not  significant.  Thus,  more  HCl  was  available  in  the  InP 
System  to  form  chlorine  rich  silicon  vapor  species  via  reaction  6-1, 
which  in  turn  lowered  the  condensed  phase  silicon  activity  relative 
to  the  GaAs  system. 

Figures  6-24  and  6-25  show  the  effect  of  temperature  on  the  de¬ 
position  zone  of  the  InP  system.  As  observed  in  the  chloride  systems, 
the  saturation  ratio  for  InP  was  much  smaller  than  that  for  GaAs.  The 
hydride  system  showed  that,  for  the  source  zone  conditions  shown  in 
the  figures,  the  InP  system  did  not  become  supersaturated  until  the 
temperature  was  below  860  K.  Alternatively,  supersaturation  at  873  K 
can  be  achieved  by  increasing  the  system  pressure  from  100  kPa  to 
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200  kPa  or  by  increasing  the  PH^  and  HCl  inlet  compositions.  These 
results  are  in  contrast  to  those  observed  in  the  InP  chloride  system 
and  both  GaAs  systems  which  displayed  supersaturated  conditions  in 
the  deposition  zone  under  all  situations  studied. 

Thermal  Oecomoosition  of  PH^  and  AsH^ 

Results  of  the  Thermal  Decomposition  Measurements 

Before  the  thermal  decomposition  rate  data  could  be  taken,  it  was 
necessary  to  characterize  the  parameter  A^,  which  represents  the  sample 
orifice  area  and  orifice  discharge  coefficient,  in  equation  5-5.  As 
this  parameter  was  expected  to  be  dependent  upon  molecular  weight  and 
temperature,  the  gases  H2,  He,  N2  end  CO2  were  used  to  determine, 
empirically,  a  correlation  for  A^  over  the  temperature  range  600  K  < 

T  <  1300  K.  For  a  nonreacting  system  equation  5-5  can  be  integrated  to 


yield 


w 
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where  P^  is  the  initial  pressure  of  species  i.  The  heat  capacity  ratio, 
k,  which  is  needed  to  evaluate  was  calculated  from  the  correlations 

shown  in  Table  6-2  and  from  using  the  relationship  Cy  =  Cp  -  R  for  an 
ideal  gas.  From  equation  6-15,  it  is  apparent  that  one  merely  needs  to 
follow  the  decay  of  pressure  in  the  reaction  tube  in  order  to  obtain 
a  value  for  A^.  This  was  done  for  the  gases  in  Table  6-2  over  the 
temperature  range  previously  mentioned  and  the  results  are  shown  in 
Figure  6-26.  The  lines  in  this  figure  represent  the  correlation 

6-16 


=  3.99  X  10'^(M^)°-°^  -  5.6  x  lO'^^T 

G  Vf 


where  has  units  of  This  correlation  predicted  values  of  A^  which 
were  lower  than  the  actual  values  for  high  molecular  weight  gases. 
However,  as  will  be  discussed  in  the  error  analysis,  this  shortcoming 
was  only  important  for  the  lowest  temperature  AsH^  data. 


The  decomposition  of  the  group  V  trihydrides  was  followed  in 
accordance  with  the  procedure  outlined  in  Chapter  Five.  The  compos¬ 
ition  versus  time  data  were  similar  for  all  of  the  group  V  tri hydrides 
studied,  thus  suggesting  a  common  mechanism  for  their  decomposition. 
This  contrasts  with  the  results  of  previous  investigators  [57-61,  65] 
who  concluded  that  the  decomposition  of  AsH^  and  PH^  were  strictly 
first  order  and  were  unaffected  by  the  presence  of  H2  in  the  system. 
The  results  presented  in  this  chapter  clearly  indicate  that  the  pre¬ 
sence  of  H2  strongly  influences  the  decomposition  rate  of  PH^  and 
AsH^. 

The  data  taken  was  insufficient  to  identify  all  of  the  reactions 
taking  place  in  the  decomposition  of  VH^  into  V2  and  V^.  However,  the 
following  segment  of  a  mechanism  suggested  by  Kelvin  [56]  for  NH^  de¬ 
composition  adequately  described  the  loss  of  from  the  system. 

VH3(g)  +  S  =  VH3(ad)  6-17 

VH3(ad)  VH(ad)  +  H2(g)  6-18 

kz 

VH(ad)  -^>  Products  (H2,  ^2' 


Equation  6-17  represents  a  Langmuir  adsorption  isotherm  [176]  which 
implies  that  the  concentration  of  VH^  adsorbed  on  the  surface  of  the 
vessel  is  at  equilibrium  with  the  vapor  phase.  The  fraction  of  sur¬ 
face  sites  covered  by  VH3  molecules  is  given  by  [176] 


KPVH3 

®  “  T+i^H3 
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where  K  is  the  equilibrium  constant 
6-13,  the  rate  of  VH3  decomposition 

'^l  '‘'^VH3 


for  reaction  6-17.  From  reaction 
is 


r  =  - 


T+kfO. 


+  k,P„  [VH(ad)] 
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where  [VH{ad)]  is  the  concentration  of  adsorbed  VH  molecules.  Assum¬ 
ing  that  CVH(ad)]  is  at  steady  state  equation  6-21  becomes 


r  =  - 


CiKP^ 

I+KPVH3 


k-jk2KPVH3PH2 


ipT^pkpH2) 


The  results  of  other  investigators  have  demonstrated  that  first  order 
kinetics  are  often  observed  for  NH^,  PH^  and  AsH^  decomposition  undef 
the  conditions  of  constant  pressure  [50-65].  Based  on  the  reaction 
rate  expression  in  equation  6-22,  apparent  first  order  behavior  can  be 
observed  only  if  KPy^g^*^^  •  T^ius, 

r  «  k^KPvH3(^PH2  +  6-23 

Combining  equation  6-23  with  equation  5-5  produces  the  following  ex¬ 
pression  for  the  total  loss  rate  of  VH^  from  the  reaction  tube 


‘^'’'^”3  _  p  r.  .,’'2  p  .  ,,-l  (RTxl/2n 

-3^  -  Ph2  ”  -7~  J 

0  w 


The  pressure  of  H2  in  the  system  was  controlled  by  the  bleed  im¬ 
posed  on  the  reactor  by  the  gas  sampling  system.  Therefore,  applying 
equation  6-15 

PH2  »  ’•Hi  e'“«2'  6-25 

where  aH2  represents  the  constants  in  equation  6-15  evaluated  for  H2  at 
the  system  temperature  and  Ph2  is  the  initial  pressure  of  H2  in  the 
reactor.  Strictly  speaking,  equation  6-25  should  contain  a  term  for  the 
production  of  H2  resulting  from  the  decomposition  of  VH^.  This  term  was 
neglected  since  most  of  the  hydrogen  in  the  system  was  present  as  a  re¬ 
sult  of  the  initial  gas  charge  (>90%  H^)  and  the  gas  sampling  system 
removed  H2  much  faster  than  it  was  produced  by  the  decomposition  reaction. 
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Substituting  equation  6-25  into  6-24  and  integrating  gives  the  result 


k,K 


l+?i445 


kTHz 


Equation  6-26  employs  four  reaction  constants  (K,  ,  k2,  k^), 

however,  only  two  parameters  (  k^K  and  k2/k2)  are  separable  for  fitting 

I 

the  rate  data  to  the  model.  For  all  t  >  0,  the  argument  of  the  logo- 
rithm  in  the  second  term  of  equation  6-26  is  always  less  than  unity. 
Therefore,  this  term  causes  the  rate  of  change  of  ln(P“/P)  to  be  less 
than  that  due  to  the  first  term  alone.  At  long  times,  the  second  term 
in  equation  6-26  becomes  approximately  constant  and  a  linear  relation¬ 
ship  between  InCP^/P)  and  t  is  obtained. 

The  observed  rate  data  are  plotted  in  Figures  6-27,  6-28  and  6-29 
for  NHj,  PHj  and  AsH^,  respectively.  These  data  were  fit  to  equation 
6-26  using  a  maximum  likelihood  generalized  least  squares  algorithm. 

The  solid  lines  represent  equation  6-26  with  the  parameters  listed  in 
Table  6-3.  For  an  irreversible  first  order  reaction,  a  plot  of  ln(P®/P) 
against  time  would  yield  a  straight  line.  This  type  of  relationship 
was  observed  for  AsH^  and  PH^  at  high  temperatures,  but  not  for  NH^ 
over  the  temperature  range  studied.  All  three  systems  studied  yielded 
straight  line  relationships  at  long  times  due  to  the  hydrogen  partial 
pressure  becoming  small  as  a  result  of  the  gas  sampling  system.  Gener¬ 
ally,  as  the  forward  reaction  rate  increased  (i.e.  k^ K  increased),  the 
amount  of  curvature  in  the  plots  decreased  due  to  the  dominance  of  the 
first  term  in  equation  6-26.  It  is  clear  from  Figures  6-28  and  6-29 
that  for  temperatures  in  excess  of  850  K  for  PH^  and  780  K  for  AsH^, 
these  decomposition  reactions  may  be  accurately  represented  by  an  ir¬ 
reversible  first  order  reaction.  Since  the  source  zones  of  the  GaAs 
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The  next  item  which  must  be  accounted  for  during  analysis  of  the 
rate  data  is  the  fact  that  the  sample  drawn  from  the  reaction  tube  may 
not  accurately  represent  the  average  concentration  in  the  tube.  Since 
these  decomposition  reactions  are  believed  to  be  heterogeneous  in  na¬ 
ture,  a  sample  obtained  from  the  vicinity  of  the  reactor  wall  may  be 
expected  to  show  a  greater  degree  of  conversion  than  one  in  the  center 
of  the  reactor.  Additionally,  if  a  mass  transfer  barrier  is  present, 
the  reaction  may  appear  to  be  first  order  due  to  the  rate  of  diffusion 
being  proportional  to  the  concentration  gradient.  This  latter  concern 
is  investigated  by  considering  the  radial  profile  in  an  infinitely  long 
cylinder  with  reaction  at  the  tube  wall.  The  analogous  heat  transfer 
problem  has  been  solved  by  Cars! aw  and  Jaeger  [175].  The  result,  trans 


formed  to  apply  to  a  mass  transport  process  is 

C, 


-  '^o^V^  e’°iV 


'(y  M  =  2AC  I  5 - 

Vl  (b^  +  A^)  J,(a^l 
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R) 


where 


A  -  «^xn/°i 
*>0  * 

and  the  eigenvalues,  a^,  are  the  roots  to  the  equation 

bJ^(b)  -  AJg(b)  *  0  6-28 

Diffusion  coefficients  for  NH^,  PH^  and  AsH-j  in  Hg,  D^,  may  be 
calculated  using  the  Chapman-Enskog  relation  [177]  and  viscosity  in¬ 
formation  [178].  The  diffusion  coefficients  for  these  gases  at  700  K 


2 

were  determined  to  be  9.5,  6.7  and  3.9  cm  /s,  respectively.  These 

3/2 

values  may  be  corrected  for  other  tanperatures  by  assuming  a  T  '  de¬ 
pendence  [177]. 

.  In  order  to  transform  the  reaction  rate  coefficients  obtained  in 
this  study  to  heterogeneous  rate  coefficients,  they  must  be  divided 
by  the  reaction  tube  surface  to  volume  ratio 

"rxn  ' 

Using  the  first  10  terms  of  equation  6-27,  the  radial  composition 
profiles  for  NH^,  PH^  and  AsH^  were  calculated  for  experimental  con¬ 
ditions  listed  in  Table  6-3.  Since  the  reaction  rate  coefficient  in¬ 
creased  much  faster  than  did  the  diffusion  coefficients,  the  greatest 
amount  of  curvature  in  the  radial  composition  profiles  was  found  to 
occur  at  the  highest  temperatures  studied.  Further,  the  curvature  in 
the  compositions  was  greatest  at  times  close  to  zero  due  to  the  com¬ 
position  dependence  of  the  reaction  rate.  Figure  6-30  shows  the  cal¬ 
culated  radial  composition  profiles  for  AsH^,  PH^  and  NH^  at  the 
highest  temperature  investigated  for  each  of  the  species,  for  a  time 
5s  into  each  experiment.  As  is  shown,  the  reaction  tube  centerline 
composition  was  approximately  1%  larger  than  the  composition  at  the 
wall.  This  small  disparity  between  the  centerline  and  wall  compositions 
indicates  that  mass  transport  was  not  a  limiting  factor  in  these  ex¬ 
periments.  Further,  the  composition  profiles  shown  in  Figure  6-30 
represent  the  greatest  deviations  from  uniform  radial  compositions  for 
all  of  the  experimental  conditions  encountered. 

The  difference  in  composition  between  a  sample  drawn  from  the  end 
of  the  reaction  tube  and  the  average  composition  in  the  reaction  tube 
is  investigated  by  determining  the  axial  composition  profile  in  the 
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tube.  This  model  assumes  that  only  diffusion  exists  in  the  tube 
(i.e.  flow  due  to  the  sample  orifice  is  neglected),  the  radial  com¬ 
position  profile  is  flat,  the  reactor  is  isothermal  and  reaction 
occurs  only  at  the  walls  (and  ends)  of  the  tube.  The  analogous  heat 
transfer  problem  has  been  solved  by  Carslaw  and  Jaeger  [175]  and  the 
result  is  presented  here  with  the  appropriate  mass  transport  proper¬ 
ties 

C(z.t)  .  C,j  cos(azn)  ^  5.30 

i\=l  [(h^  +  a^)l  +  h]  cos(la^) 

where:  h  = 

1=1/2  reactor  length 

The  eigenvalues  represented  by  a^^  are  the  roots  of  the  implicit 


equation 


a„  tan(la„)  =  h 
n  n 
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Figure  6-31  shows  the  axial  composition  profiles  for  AsH^,  PH^ 
and  NH^  at  816  K  in  the  reaction  tube  as  calculated  from  equation  6-30 
using  the  first  10  eigenvalues.  The  average  composition  in  the  tube 
corresponds  to  99%  conversion  of  the  VH^  initially  present  in  the  tube. 
The  time  required  to  reach  this  level  of  decomposition  was  determined 
to  be  190  s  for  AsH^,  1127  s  for  PH^  and  2.08x10^  s  for  NH^.  The  large 
gradient  in  AsH^  composition  is  due  to  the  high  rate  of  reaction  at 
this  temperature  and  the  small  diffusion  coefficient  for  VH^  relative 
to  PH^  and  NH^.  The  calculated  PH^  and  NH^  axial  composition  profiles 
at  909  K  and  1048  K  were  much  more  exaggerated  than  those  shown  in 
Figure  6-29,  but  were  less  dramatic  than  that  of  AsH^  at  816  K. 
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Due  to  the  shape  of  the  axial  composition  profile  at  high  re¬ 
action  rates,  the  effective  first  order  rate  constants  (k^K)  shown 
in  Table  6-3  are  greater  than  those  averaged  over  the  reaction  tube 
length.  The  largest  discrepancies  were  4.3«  for  AsH^  at  816  K,  3.3% 
for  PH^  at  909  K  and  2.8%  for  at  1048  K.  The  k^K  values  corrected  for 
axial  composition  profiles  are  presented  in  Table  6-4. 

The  uncertainties  associated  with  the  k^K  values  listed  in  Table 
6-4  are  due  to  the  uncertainty  associated  with  measuring  the  partial 
pressure  of  the  gas  with  the  mass  spectrometer,  the  presence  of  the 
inlet  tube  which  is  not  at  the  reaction  tube  temperature  and,  in  the 
case  of  AsHj,  the  underprediction  for  the  value  of  A^  by  equation 
6-16.'  Very  little  background  signal  was  present  in  the  m/e  ranges 
used  for  the  investigation  of  AsH^  and  PH^  decomposition.  Thus,  by 
making  repetitive  scans  of  the  fragment  ion  patterns  of  these  VH^ 
molecules  at  low  temperature,  it  was  determined  that  their  presence 
could  be  detected  to  within  an  uncertainty  of  65  Pa  in  the  reaction 
chamber.  Since  the  initial  pressure  of  these  gases  was  typically 
9200  Pa  and  the  final  pressure  after  decomposition  was  approximately 
130  Pa,  the  uncertainty  varied  from  a  minimum  0.7%  to  50%  for  each 
datum.  At  low  temperature,  where  the  reaction  rates  were  slowest, 
most  of  the  data  taken  was  at  relatively  high  VH^  partial  pressure 
and,  therefore,  the  uncertainty  of  these  data  is  the  smallest.  High 
temperature  and  fast  reaction  rates  require  that  many  data  points  at 
small  VH^  partial  pressure  be  taken.  This  causes  the  observed  in¬ 
crease  in  uncertainty  as  the  reaction  rate  increases. 

The  measurement  of  NH^  partial  pressure  was  less  precise  than 
that  of  PH^  and  AsH^  due  to  the  background  created  by  the  desorption 


Table  6-4 

Corrected  Rate  Coefficient  Results 


T(K) 

k^KCs'"’) 

750 

0.00284 

848 

0.00336 

898 

0.00405 

949 

0.0060 

992 

0.0114 

1023 

0.288 

1048 

0.0827 

712 

0.00167 

815 

0.00386 

846 

0.00929 

855 

0.0107 

909 

0.0424 

712 

0.00156 

726 

0.00311 

748 

0.00519 

762 

0.00798 

779 

0.00862 
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of  H2O  from  the  walls  of  the  second  stage  of  the  vacuum  system. 
Fluctuation  in  the  background  signal  caused  the  uncertainty  in  the 
measurement  of  partial  pressure  to  be  130  Pa.  This  uncertainty 
is  exacerbated  by  the  smaller  concentration  of  NH^  in  the  gas  cylin¬ 
der  (4. 3")  compared  to  that  of  PH^  and  AsK^  (10%).  The  maximum 
measured  partial  pressure  of  NH^.in  the  reaction  tube  was  5500  Pa 
which  yields  an  uncertainty  of  2.4%.  The  measurement  of  130  Pa  of 
NH^  carries  with  it  an  uncertainty  of  100%.  Therefore,  the  uncertainties 
associated  with  determining  the  compositions  are  much  higher  than 
those  for  PH^  and  AsH^. 

The  tube  provided  for  inlet  of  the  gases  to  the  reaction  chamber 
contained  a  volume  of  gas  which  was  0.4%  of  the  volume  of  the  reaction 
chamber.  The  surface  area  of  the  tube  was  approximately  5%  of  that  in 
the  reaction  chamber.  An  uncertainty  of  2%  was  assigned  to  each  rate 
constant  as  a  result  of  the  presence  of  this  tube. 

The  failure  of  equation  6-16  to  accurately  predict  A^  for  high 
molecular  weight  gases  causes  additional  uncertainty  in  the  values 
for  AsHj.  Based  on  the  results  for  CO2  shown  in  Figure  6-26,  it 
appears  that  the  predicted  A^  value  for  AsH^  should  be  at  most  10% 
low.  The  loss  of  gas  due  to  the  sampling  orifice  represents  1/3  of 
the  first  term  in  equation  6-21  for  AsH^  decomposition  at  712  K  and 
even  less  at  the  higher  temperatures.  An  uncertainty  of  3%  was 
added  to  the  uncertainty  at  712  K  and  1%  was  added  to  the  uncertain¬ 
ties  at  726  K,  748  K,  762  K,  and  779  K. 


Determination  of  Activation  Enerc 


The  temperature  dependence  of  the  constant  k^K  is  determined  by 
assuming  an  Arrhenius  type  dependence  for  k^ 
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and  relating  the  equilibrium  constant  K  to  the  enthalpy  and  entropy 
of  adsorption  through  the  Gibbs  energy 

Thus,  the  temperature  dependence  becomes 

k^K  *  (kQe^^ad/R)e'^^a  ^Had)/'^'^ 
where  E.  is  the  activation  energy  of  the  forward  reaction  6-18.  The 

a 

data  in  Table  6-4  were  fit  to  equation  6-34  and  the  result  is  shown 
in  Figure  6-32  and  Table  6-5.  Due  to  the  strong  influence  the  pre¬ 
sence  of  H2  had  on  the  decomposition  reactions  for  NH^  and  PH^,  the 
low  temperature  rate  constants  were  in  poor  agreement  with  equation 
6-33.  Therefore,  only  the  four  highest  temperature  reaction  rate  con¬ 
stants  were  included  in  the  fit  for  these  gases. 

The  apparent  activation  energies  found  by  other  investigators  for 
the  decomposition  of  NH^,  PH^  and  AsH^  on  quartz  are  sunnmarized  in 
Table  6-6.  The  apparent  activation  energy  found  for  AsH^  in  these 
experiments  is  in  reasonable  agreement  with  that  of  Frolov  et  al .  [65] 
who  used  an  open  system  and  infrared  spectrophotometer.  The  apparent 
activation  energy  found  for  PH^  is  approximately  10  kcal/mole  lower 
than  that  found  by  other  investigators.  The  results  of  van't  Hoff  and 
Kooj  [57]  and  Hinshelwood  and  Topley  [59]  were  based  on  raanometric 
methods  and,  therefore,  are  expected  to  over  estimate  the  activation 
energy.  This  is  due  to  reactions  in  the  decomposition  sequence  which 
may  not  have  achieved  equilibrium.  The  method  employed  by  Devyatykh 
et  al .  [60]  was  not  presented. 

The  apparent  activation  energy  found  for  the  decomposition  of 
NH^  was  significantly  greater  than  that  found  by  previous  investi¬ 
gators  [52,  53,  55].  All  of  these  investigators  employed  manometric 
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Table  6-5 


Pre-exDonential  Factors  and  Apparent  Activation  Energies 


Gas 

k^eaSad/R(3-l, 

Ea  +  AH^^(kcal/mole) 

Uncertainty  {%) 

NH3 

2.92x10^^ 

60.2 

19 

PH3 

2.42x10^ 

36.5 

9 

ASH3 

1 .62x10^ 

29.2 

7 

Table  6-6 

Apparent  Activation  Energies  for  VH.,  Decomposition  on  Quartz 

"  ^ 

Gas 

Ea  +  AH^^(kcal/mole) 

Reference 

ASH3 

32.6 

65 

PH3 

46.4 

57 

PH3 

46  ±  4 

59 

'’"a 

44.2 

60 

NH3 

43  ±  5 

52 

NH3 

38.2 

53 

NH3 

33.7  ±  2 

55 
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methods  for  their  experiments,  but,  since  N2  and  H2  are  the  only  ex¬ 
pected  products  of  the  decomposition  reaction,  it  is  doubtful  that 
their  values  are  systematically  high.  Only  4  data  points  were  used 
to  determine  the  apparent  activation  energy  for  NH^  decomposition  in 
this  study  and  the  highest  point  had  a  very  large  uncertainty  (26%). 
Thus,  the  uncertainty  associated  with  the  NH^  results  may  be  larger 
than  was  estimated. 

Assuming  that  the  rate  constants  IC2  and  have  Arrhenius  type 
temperature  dependencies,  it  is  expected  that  a  plot  of  In(k2/k2) 
versus  time  would  yield  a  straight  line.  Further,  the  slope  of  this 
line  should  be  the  difference  between  the  activation  energies  for 
the  two  reactions.  This  plot  is  shown  in  Figure  6-33.  The  AsH^ 
results  and  the  high  temperature  NH^  results  show  reasonable  agree¬ 
ment  with  the  suggested  model.  The  PH^  results  and  low  temperature 
NH^  results  are  anomalous.  Based  on  Figure  6-33,  the  activation 
energy  differences  were  59.6,  8.0  and  22  kcal/mole  for  NH^j,  PH^  and 
AsHj,  respectively. 

Implications  of  the  Measurements  on  Other  Reactions  in  the  Decom- . 

position  Chain 

The  relative  peak  intensities  for  the  mass  spectrometer  frag¬ 
ment  ion  pattern  of  NH^,  PH^  and  AsH^  are  shown  in  Table  6-7.  The 
PH^  and  AsH^  patterns  found  in  this  study  agree  reasonably  well  with 
those  in  the  reference  [179].  The  poor  agreement  for  the  NH^  results 
is  due  to  the  presence  of  the  H2O  background  which  dominates  the  sig¬ 
nal  at  m/e  *  17  due  to  the  OH^  fragment  ion.  Generally,  some  dis¬ 
crepancies  are  expected  in  the  relative  peak  intensities  reported  by 
different  institutions  due  to  the  differences  in  operating  parameters 
and  geometries  of  each  mass  spectrometer. 
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The  mass  spectrometer  signals  for  PH^,  P  ,  P2  and  P^  are  shown  in 
Figure  6-34  for  the  decomposition  of  PH^  at  855  K.  Since  the  fragment 
ion  pattern  of  PH^  was  known,  the  signal  for  P*  was  corrected  for  the 
presence  of  PH^.  Unfortunately,  the  fragment  ion  patterns  for  P^  and 
Pg  were  not  known  and  therefore,  it  was  not  possible  to  correct  the  P^ 
and  Pg  signals  for  the  presence  of  P^  and  the  P^  signals  for  the  pre¬ 
sence  of  P2.  If  these  patterns  were  known,  it  would  be  possible  to 
determine  the  rates  of  formation  of  P,  P2  and  P^  from  the  decomposi¬ 
tion  of  PH^.  Only  a  very  small  signal  corresponding  to  P^  was  ob¬ 
served. 

Figure  6-35  depicts  the  decomposition  of  ASH2  at  779  K  and  shows 
the  mass  spectrometer  signals  for  AsH^.  As^,  As^  and  As^.  The  signals 
in  the  As  system  are  smaller  than  those  in  the  P  system  due  to  the  de¬ 
creased  electron  multilplier  gain  at  large  m/e  values.  This  is  clearly 
demonstrated  for  AsJ  (m/e  =  300)  where  the  maximum  signal  observed 
(2mV)  was  barely  recognizable  over  the  noise  at  the  baseline.  As  in 
the  P  system,  the  As^  signal  was  corrected  for  the  presence  of  AsH^. 
but  the  influence  of  AS2  was  not  accounted  for  since  its  fragment  ion 
pattern  was  unknown.  No  signal  corresponding  to  As^  was  observed. 

In  both  the  As  and  P  systems,  the  V^,  and  signals  reach 
maxima  at  the  same  time.  This  may  be  due  to  the  and  Vg  signals 
being  fragments  of  since  is  thermodynamically  favored  at  the 
temperatures  and  pressures  employed  in  the  reaction  tube.  If  these 
maxima  did  not  occur  in  coincidence,  one  would  be  compelled  to  con¬ 
clude  that  additional  rate  influential  reactions  may  exist  in  the 
formation  of  and  from  the  decomposition  of  VH^.  Since  this  is 
not  the  case,  it  appears  reasonable  to  assume  that  the  decomposition  of 
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Figure  6-35 

AsH,  Decomposition  at  779K 


VH^Cad)  into  VH(ad)  is  the  rate  limiting  step  in  the  reaction  se¬ 
quence  and  that  V,  V2  and  may  be  assumed  to  be  at  equilibrium. 

Comparison  Between  Equilibrium  and  Kinetic  Models 
for  VHj  Decomposition 

The  applicability  of  an  equilibrium  model  to  describe  the  de¬ 
composition  of  VH^  in  a  hydride  system  reactor  can  best  be  investi¬ 
gated  by  modelling  the  group  V  source  zone.  Assuming  an  isothermal 
reactor,  cylinder  geometry,  plug  flow  and  reaction  only  at  the  re¬ 
actor  wall,  the  steady  state  axial  and  radial  VH^  composition  pro¬ 
files  are  found  by  solving  [177] 
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where  C  is  the  concentration  of  VH^,  v  is  the  velocity  of  the  plug  and 
D  is  the  diffusion  coefficient  for  VH^  in  Hg.  Assuming  that  C  may  be 
separated  into  r  and  2  components  (C  =  the  partial  differential 

equation  6-35  may  be  written  as  two  ordinary  differential  equations 


=  0 


6-36 


+ 


avr^  r 
~5~  ‘'r 


=  0 


6-37 


subject  to  the  boundary  conditions: 
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The  axial  composition  profile,  represented  by  equation  6-36,  has  the 


solution 


*  C^e^^  +  Cge'"-^ 


where 


-  vCl  ±  (1 


Applying  equation  6-39  requires  =  0  therefore, 


where  the  constant  is  evaluated  in  conjunction  with  the  radial 


solution. 


Equation  6-37  is  Bessel's  equation  which  has  the  solution 

where  J  and  are  the  Bessel  functions  of  the  first  and  second  kind, 

respectively.  Since  becomes  infinite  at  r=0,  the  boundary  condition 

represented  by  equation  6-41  requires  that  vanish.  Thus, 

Cr  “ 

where  the  eigenvalue,  a,  is  determined  from  the  solution  of  equation 
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Equation  6-47  is  rewritten  as 
Jl(aR)  '  ° 


which  has  an  infinite  number  of  roots.  Thus,  the  solution  to  equation 
6-35  is 


where  the  constant  C2  has  been  included  into  A^.  The  constants  A^  are 
evaluated  by  applying  the  principal  of  orthogonality  which  yields  the 
resul t 
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which  was  evaluated  by  numerical  integration.  Using  five  terms  of 
the  series  in  equation  6-49  allowed  the  solution  to  converge  to  with¬ 
in  at  least  three  significant  figures  for  all  of  the  cases  studied. 

Figures  6-36  and  6-37  show  the  residence  time  required  to  achieve 
99»  decomposition  of  ASH2  and  PH^  for  various  reactor  diameters.  This 
degree  of  decomposition  was  chosen  because  the  equilibrium  compositions 
of  V2  and  remain  essentially  unaffected  by  any  further  decomposition 
of  VH^.  Thus,  if  a  hydride  CVD  system  is  designed  which  yields  a  VH^ 
decomposition  of  at  least  99S,  an  analysis  of  the  gas  phase  based  on 
the  assumption  of  equilibrium  will  provide  accurate  results  for  the 
and  compositions.  As  is  shown,  higher  temperatures  require  shorter 
tube  lengths  in  order  to  reach  this  level  of  conversion.  Also,  AsH^ 
decomposes  much  faster  than  PH3  for  a  given  temperature  and  tube  dia¬ 
meter.  The  increase  in  residence  time  with  tube  diameter  is  due  to  the 
fact  that  the  volume  of  gas  per  unit  length  of  tube  increases  faster 
than  the  surface  area  as  the  tube  diameter  increases  (i.e.  for  a  cylin¬ 
der  the  surface  to  volume  ratio  goes  as  l/radius).  Thus,  the  reaction 
rate  per  volume  decreases  and  therefore  longer  times  are  required  to 
reach  the  same  degree  of  decomposition. 

The  importance  of  including  axial  diffusion  in  the  model  is  also 
demonstrated  by  Figures  6-36  and  6-37.  The  second  derivative  with 
respect  to  axial  position  in  equation  6-35  was  set  to  zero  and  the 
residence  times  were  calculated  at  873  K  as  functions  of  tube  diameter. 
As  shown,  the  model  neglecting  axial  diffusion  grossly  underpredicts 
the  residence  time  (  and  therefore  the  reactor  length)  necessary  to 
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achieve  the  desired  degree  of  decomposition.  These  results  were  cal¬ 
culated  at  a  velocity  of  1  cm/s  which  corresponds  to  0.1  <  Pe  <  2 
(Pe  =  2Rv/0).  At  these  Peclet  numbers  axial  diffusion  is  extremely 
important.  For  systems  employing  higher  velocities,  and  therefore 
larger  Peclet  numbers  (e.g.  Pe  >  10),  the  importance  of  including  axial 
diffusion  is  diminished  [180]. 

Hsieh  [181]  has  studied  the  Graetz  problem  for  heat  transfer  [182] 
(which  includes  a  parabolic  velocity  profile)  and  compared  the  results 
to  a  model  using  plug  flow.  He  found  that  the  plug  flow  model  gave  re¬ 
sults  which  were  conservative  (i.e.  predicted  longer  reactor  lengths 
to  achieve  the  same  average  radial  temperature)  when  compared  to  the 
solution  to  the  Graetz  problem.  Based  on  these  results,  it  is  believed 
that  Figures  6-36  and  6-37  are  also  conservative  in  nature.  Typical 
hydride  system  CVO  reactors  employ  a  tube  for  the  group  V  source  zone 
which  is  approximately  2  cm  in  diameter  and  use  a  gas  velocity  of  1 
cm/s.  Based  on  Figure  6-36  and  6-37,  tube  lengths  of  10  cm  and  30  cm 
are  required  to  reach  99%  VH^  decomposition  at  973  K  for  AsH^  and  PH^, 
respectively.  At  873  K,  these  tube  lengths  must  be  increased  by 
approximately  a  factor  of  three.  Many  hydride  system  reactors  are  be¬ 
tween  50  cm  and  100  cm  long.  Therefore,  it  is  expected  that  an  equi¬ 
librium  analysis  for  a  system  using  AsH^  will  yield  acceptable  results 
if  the  source  and  mixing  zone  temperatures  are  above  873  K.  Source  and 
mixing  zone  temperatures  greater  than  973  K  are  necessary  in  order  for 
an  equilibrium  analysis  to  apply  to  a  system  using  PH^. 

The  AsHj  radial  composition  profiles,  99%  decomposed,  are  shown 
for  temperatures  of  873  K,  973  K  and  1073  K  in  Figure  6-38  in  a  10  cm 
diameter  reactor.  The  large  diameter  reactor  was  used  in  this 


calculation  in  order  that  some  curvature  in  the  radial  composition 
profile  might  be  observed.  The  more  pronounced  curvature  at  high 
temperatures  is  due  to  the  reaction  rate  increasing  faster  than  the 
diffusion  coefficient  with  increasing  temperature.  Thus,  at  high 
temperatures  in  large  diameter  tubes,  mass  transport  limitations  are 
expected  to  become  increasingly  important. 


CHAPTER  SEVEN 

CONCLUSIONS  AND  RECOMMENDATIONS 
Conclusions 

Parametric  analyses  of  the  chemical  reaction  thermodynamics 
pertaining  to  the  chemical  vapor  deposition  of  GaAs  and  InP  in  the 
chloride  and  hydride  systems  have  revealed  some  dissimilarities  be¬ 
tween  these  systems.  Within  the  chloride  system  the  condensed  phase 
silicon  activity  resulting  from  the  use  of  solid  III-V  material  as 
the  group  III  source  was  shown  to  be  smaller  than  that  obtained  when 
a  liquid  group  III  source  was  used.  This  difference  was  more  pro¬ 
nounced  for  the  GaAs  system  than  it  was  for  the  InP  system.  The  sili¬ 
con  activity  for  the  InP  system  using  a  liquid  source  was  slightly 
less  than  that  of  the  analogous  GaAs  system,  while  the  GaAs  system 
using  the  solid  source  exhibited  a  much  smaller  silicon  activity  than 
did  the  InP  system  using  a  solid  source.  Small  additions  of  HCl , 

H^O  or  VCl^  to  the  mixing  zones  of  these  chloride  systems  were  very 
effective  in  reducing  the  silicon  activity.  Replacing  the  hydrogen 
carrier  gas  with  an  inert  gas  was  found  to  be  an  effective  method  of 
reducing  the  silicon  activity  only  if  greater  than  98%  of  the  hydro¬ 
gen  was  replaced.  The  degree  of  supersaturation  in  the  deposition 
zone  was  much  larger  for  the  GaAs  system  than  for  the  InP  system. 

Also,  the  use  of  a  liquid  group  III  source  material  led  to  a  greater 
degree  of  supersaturation  in  the  deposition  zone  than  did  the  use  of 
a  solid  III-V  compound  as  the  group  III  source. 


The  hydride  system  was  found  to  be  substantially  less  super¬ 
saturated  than  the  chloride  system  under  typical  operating  conditions 
and  the  condensed  phase  silicon  activity  was  less  by  approximately 
one  order  of  magnitude  for  the  hydride  system  when  compared  to  the 
chloride  system.  Replacement  of  the  hydrogen  carrier  gas  with  an 
inert  was  not  as  effective  a  method  of  reducing  the  silicon  activity 
in  the  hydride  system  as  it  was  in  the  chloride  system  due  to  the  H2 
liberated  from  the  decomposition  of  HCl  and  VH^.  The  condensed  phase 
silicon  activity  was  found  to  decrease  with  increasing  pressure  in  the 
deposition  zones  while  maxima  in  the  silicon  activity  were  found  in 
the  source  and  mixing  zones  at  approximately  4  kPa  pressure  for  sys¬ 
tems  using  liquid  source  materials.  Increasing  the  temperature  was 
found  to  strongly  increase  the  silicon  activity  in  the  source  and 
mixing  zones.  The  influence  of  inlet  composition  on  silicon  activity 
was  not  as  strong  as  that  of  temperature,  pressure,  replacement  of 
by  inerts  or  the  addition  of  species  downstream  of  the  source  zone. 
Reduced  silicon  activities  occurred  as  the  inlet  composition  of  chlo¬ 
rine  containing  species  was  increased. 

In  general,  it  was  found  that  shifting  the  vapor  phase  silicon 
species  from  those  which  are  hydrogen  rich  to  chlorine  rich  species 
markedly  reduced  the  condensed  phase  silicon  activity.  This  effect 
was  so  pronounced  that  even  in  systems  where  the  total  amount  of  si¬ 
licon  in  the  vapor  phase  as  chlorinated  species  was  much  greater  than 
that  of  hydrogen  rich  species  in  other  systems,  the  condensed  phase 
silicon  activity  was  still  sharply  reduced. 

Investigations  into  the  decomposition  kinetics  of  PH^  and 
AsH^  have  revealed  that  the  mechanisms  of  decomposition  for  these 
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molecules  appear  to  be  identical  in  at  least  the  first  steps  of  the 
reaction  sequence.  The  first  step  consists  of  Langmuir  adsorption 
onto  the  surface  of  the  vessel  followed  by  a  reversible  reaction 
between  the  adsorbed  VH^  molecule  and  its  products  VH  and  Hg.  The 
third  step  appears  to  be  an  irreversible  decomposition  of  the  VH 
molecule.  The  second  step  in  the  VH^  decomposition  reaction  se¬ 
quence  is  rate  limiting  and,  over  the  temperature  range  750  K  <  T 
<  1048  K,  the  reversibility  of  this  reaction  must  be  included  in  the 
decomposition  of  NH^.  It  was  found  that  for  temperatures  greater 
than  780  K  for  AsH^  and  850  K  for  PH^,  these  decompositions  could  be 
accurately  modeled  as  irreversible  first  order  reactions.  Below  these 
temperatures,  the  reverse  reaction  in  the  second  step  of  the  sequence 
must  be  included.  The  apparent  activation  energies  for  these  decom¬ 
position  reactions  were  determined  to  be  60.2,  36.5  and  29.2  kcal/mole 
for  NH^,  PH^  and  AsH^,  respectively.  Indirect  evidence  supports  the 
assumption  that  the  products  of  these  decomposition  reactions  (V,  V2 
and  V^)  may  be  considered  to  be  at  equilibrium. 

An  isothermal  plug  flow  model  for  the  group  V  hydride  source  zone 
demonstrated  that  the  application  of  an  analysis  based  on  the  assump¬ 
tion  of  equilibrium  is  justified  for  temperatures  greater  than  873  K 
in  the  AsH^  system  and  973  K  in  the  PH^  system.  This  conclusion  is 
based  on  a  gas  velocity  of  1  cm/s  in  a  2  cm  diameter  tube.  The  plug 
flow  model  also  showed  that  axial  diffusion  must  be  included  in  the 
model  if  accurate  results  are  to  be  obtained  when  Pe  <  10. 

Recommendations 

The  results  of  this  study  provide  a  basis  for  further  investi¬ 
gations  into  the  fundamental  chemistry  associated  with  the  hydride 
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and  chloride  processes.  It  is  recommended  that  the  decomposition 
rate  of  the  various  VH^  species  be  determined  at  other  surface  area 
to  volume  ratios  than  that  which  was  employed  in  these  experiments. 
This  information  is  necessary  to  quantify  the  heterogeneous  nature 
of  these  reactions.  Additionally,  the  fragment  ion  patterns  for  the 
V2  and  species  need  to  be  determined  in  order  to  further  charac¬ 
terize  the  decomposition  of  VH^. 

The  current  disposition  of  the  mass  spectrometer  and  reaction 
tube  is  such  that  many  investigations  pertinent  to  the  chemistry  of 
III-V  CVD  processes  may  be  carried  out  without  significant  modi¬ 
fications,  Fragment  ion  patterns  for  the  '^4  species  can  be 

obtained  by  placing  solid  group  V  material  into  the  reactor,  thus 
making  an  effusion  cell.  The  decomposition  kinetics  of  ASCI2.  PClj 
and  the  group  III  metalorganics  can  be  studied  using  the  same  methods 
which  were  employed  for  VH^.  A  bubbler  system  would  need  to  be  in¬ 
stalled  in  order  to  provide  a  source  for  the  VCl^  species. 

It  is  recommended  that  the  mass  spectrometer  system  be  upgraded 
to  a  true  molecular  beam  facility,  as  is  described  in  Appendix  B. 

The  resulting  system  would  be  expected  to  perform  with  sufficient 
sensitivity  (~1  ppm)  to  allow  the  detection  of  minor  species  in  the 
reaction  tube.  The  rate  of  formation  of  vapor  phase  silicon  species, 
arising  from  the  interaction  between  Hg,  HCl  and  the  reactor  wall, 
could  then  be  determined  and  compared  to  the  calculated  equilibrium 
results.  Conversion  of  the  current  constant  volume  reactor  to  an 
open  system  would  allow  the  study  of  actual  hydride  and  chloride  CVD 
systems.  Once  again,  the  mass  spectrometer  may  be  employed  to  deter¬ 
mine  the  vapor  phase  composition  and  these  results  could  then  be  com¬ 
pared  to  the  thermodynamic  results. 


Finally,  the  characterization  of  the  chemistry  surrounding 
the  epitaxial  growth  of  III-V  materials  requires  that  the  principal 
species  taking  part  in  the  deposition  reaction  be  identified.  Thus 
measurement  of  the  vapor  composition  in  the  immediate  vicinity  of 
the  substrate  (e.g.  through  a  hole  etched  in  the  substrate)  is  ex¬ 
pected  to  produce  local  information  directly  related  to  the  depos¬ 
ition  reaction.  Changes  in  operating  parameters  which  affect  the 
deposition  process  (e.g.  growth  rate,  epi-layer  quality)  can  then 
be  related  directly  to  the  local  chemistry  at  the  substrate  surface 
This  information  is  instrumental  for  the  understanding  of  III-V  CVD 


APPENDIX  A 

A  STOICHIOMETRIC  ALGORITHM  FOR  THE  CALCUUTION  OF  CHEMICAL  EQUILIBRIA 
The  stoichiometric  algorithms  employed  by  Cruise  [183]  and  Smith 
and  Missen  [184]  may  be  extended  to  include  a  solution  phase  in  addi¬ 
tion  to  vapor  and  pure  condensed  phases.  A  stoichiometric  algorithm  is 
one  in  which  a  set  of  formation  reaction  equations  are  used  to  provide 
a  driving  force  in  the  numerical  solution  of  the  equilibrium  problem. 
The  formation  reaction  for  any  species  i  in  a  system  containing  E  ele¬ 
ments  is  written  as 


A-1 


The  reactants  in  equation  A-1  are  called  the  basis  species  (thus, 
the  superscript  "b")  and  form  a  linearly  independent  set  of  species  from 
within  the  system  which  represent  all  of  the  elements  present  in  the 
system.  The  stoichiometric  coefficient  describes  the  number  of  mole¬ 
cules  of  basis  species  k  which  are  required  in  the  formation  of  species 
i . 


The  technique  used  in  this  stoichiometric  algorithm  to  calculate 
the  equilibrium  composition  of  the  system  is  as  follows.  First,  the 
equilibrium  constants  for  the  formation  reactions  in  equation  A-1  are 
calculated  from  the  Gibbs  energy  change  df'the  reactions.  Next,  the 
current  compositions  for  each  of  the  species  present  are  used  to  cal¬ 
culate  another  "equilibrium  constant".  The  values  of  the  two  constants 
are  then  used  in  finite  difference  equations  which  approximate  deriva¬ 
tives  of  the  equilibrium  constants  with  respect  to  the  reaction  extents. 
These  equations  predict  changes  in  the  reaction  extents  which  will 
yield  an  improved  approximation  to  the  equilibrium  composition  of  the 
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The  discrete  formulation  of  this  algorithm,  applied  to  a  system 
consisting  of  a  vapor  phase  containing  V  species,  a  solution  phase 
containing  S  species  and  C  pure  condensed  phases,  starts  with  the  de¬ 
finition  of  the  elemental  abundance  matrix  A,  a  vector  representing 
the  molar  amount  of  each  of  the  species  N,  and  a  vector  containing  the 
total  moles  of  each  element  present  in  the  system 

^  ^  r  ^ll'"  ^iE 
~  ^vsc,l  ■“^vsc,E  ■ 

-  ^"l  •••  "vsc^ 

=  (b,  ...  bg) 

where:  a..  =  number  of  atoms  of  element  j  in  species  i 

*  J 

n.  =  moles  of  specie  i  in  the  system 

b.  »  moles  of  element  j  in  the  system 

J 

vsc  »  V+S+C 

E  *  number  of  elements  in  the  system 
These  definitions  result  in  the  following  relationship  which  de¬ 
scribes  the  system  mass  balance. 

N  =  B  A-2 

In  order  to  construct  a  set  of  formation  reactions,  as  in  equation 
A-1 ,  a  set  of  basis  species  which  is  linearly  independent  from  each 
other  yet  represents  all  of  the  elements  present  in  the  system,  must  be 
found.  A  further  restriction  is  placed  upon  the  set  of  basis  species 
when  the  convergence  rate  of  the  algorithm  is  considered.  Obviously, 
the  molar  amounts  of  each  of  the  species  must  be  nonnegative  for  a 
physically  realistic  situation  to  exist.  If  a  basis  species  which  has 
a  very  small  initial  composition  is  chosen,  and  the  formation  of  an¬ 
other  species  requires  the  consumption  of  this  basis  species,  then  the 
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rate  of  convergence  will  be  slow  in  order  to  prevent  the  molar  amount  of 
the  basis  species  from  becoming  negative.  Thus,  a  rapidly  converging 
algorithm  will  employ  an  optimum  set  of  basis  species  which,  in  addi¬ 
tion  to  the  two  previously  mentioned  requirements,  are  present  in  the 
greatest  molar  amounts  available.  Also,  since  the  molar  amounts  of  each 
species  change  after  each  iteration,  this  optimum  set  of  basis  species 
may  need  to  be  chosen  at  each  Iteration. 

The  optimum  set  of  basis  species  is  chosen  by  first  sorting  vector 
^  into  descending  order  (n^>n2>...n)  and  ordering  the  rows  of  A  corre¬ 
spondingly.  The  first  row  of  A  is  then  transferred  into  the  first  row 
of  a  new  matrix  0.  The  second  row  of  A  is  then  transferred  into  the 


next  row  of  D  and  tested  for  linear  independence.  If  it  is  linearly 
independent,  the  row  is  kept  and  the  process  continues.  If  not,  the 
next  row  in  A  is  tried.  This  process  is  continued  until  D  is  filled. 
For  a  system  comprised  of  E  elements,  there  will  be  E  basis  species, 
therefore  matrix  £  will  always  be  square. 

An  efficient  method  of  testing  for  linear  independence  between  the 
rows  of  the  D  matrix  is  to  build  the  Gram-Schmidt  orthogonal i zed  matrix 
0^'  using  the  formulation  [185] 


i-1  I  d..  d,'. 

I!.  =  d..-  T  d^'. 

IJ  E  - 


k=l 


where:  i  =  2,3,...md 


j  »  1 ,2, . . .E 

md  *  current  row  in  matrix  0 


A  linear  dependence  between  row  md  of  £  and  the  remaining  rows  in 
0  exists  if  row  md  of  O'  contains  all  zeros. 


VT*  v-v 


The  matrix  containing  the  stoichiometric  coefficients  for  all  of 
the  formation  reactions  is  given  by 


'"n  •••  '^’E 


vsc,l  *■*  vsc,E 


=  A  D 


Equilibrium  constants  may  now  be  calculated  for  each  of  the 
formation  reaction  equations  using  the  Gibbs  energy  change  for  each 
reaction. 

"eqi  'RT  "k  -  “P  *-5 

where:  =  equilibrium  constant  for  reaction  i 

u?  =  standard  chemical  potential  of  species  i 

*  standard  chemical  potential  of  basis  species  k 
An  "equilibrium  constant"  for  each  formation  reaction  may  also  be 
calculated  from  the  current  estimate  to  the  equilibrium  composition. 


In  Q.  =  In  a-  -  v..  In  a.  A-6 

1  1  IK  K 

where:  a^.  =  y.  n-  Pp^/n-j-  (activity  of  species  i) 

=  activity  coefficient 

=  total  moles  of  species  in  the  same  phase  as  specie  i 

p  =  {P/Po  vapor  species 
^pi  “  ’■1  nonvapor  species 

Pjj  =  reference  state  pressure 
Defining  the  extent  of  the  reaction,  as 
£l  .  (aj  - 

where:  n^  =  composition  of  species  j  at  reaction  extent  Ej 

«J  * 

the  "equilibrium  constant"  as  calculated  from  the  current  composition 
may  be  written  as 
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In  Q,  .  In  [B,(n,  ♦  C,)]  lnCSi(n^  *  v,,,?,)] 


A-7 


where:  8^  *  Y^-  Pp^/nx 

As  an  approximation,  assume  that  6^.  and  Sj^  are  constant  for  small 
changes  in  Introducing  a  function  6.  to  account  for  nonvapor  phases, 
a  finite  difference  approximation  which  relates  changes  in  the  extent  of 
reaction  to  the  difference  between  actual  equilibrium  constant  and  the 
"equilibrium  constant"  calculated  from  the  current  composition  is  ob¬ 
tained  by  differentiation  of  equation  A-7 


In  K 


eqi 


A^i 


InQ.  ^dlnq.  6^4 


A-8 


n-  lc*l 


.  _  r  r  _  vapor  phase 

where:  5-,  -  {q  phase 

Thus  the  difference  between  the  equilibrium  constant  for  formation 
reaction  i  as  calculated  from  the  Gibbs  energy  and  as  calculated 

from  the  current  approximation  to  the  equilibrium  composition  (Q.)  re¬ 
sults  in  a  prediction  of  which  will  yield  a  better  approximation  to 
the  system  equilibrium  composition- 

In  order  to  prevent  the  occurrence  of  negative  molar  amounts,  the 
effect  of  on  the  new  composition  must  be  tested.  Thus,  for  the 
nonbasis  species 


nj  .  n,  as; 

and  for  the  basis  species 
vsc 

"k '  "k  -  ,1,  ''ik 

where;  AC^  =  <AE^ 
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In  order  to  assure  that  nonbasis  species  remain  positive,  A5^  is 
set  equal  to  -n.  if  AE^<0  and  lAE^|>  n^.  Nonnegativity  of  the  basis 


species  is  assured  by  applying  the  relation 
vsc 


A-11 


to  each  of  the  basis  species  in  the  system  and  choosing  the  smallest 
positive  value  of  <  which  results  while  tc  is  constrained  to  be  no 
greater  than  unity. 

In  order  to  assure  convergence  of  the  numerical  scheme,  it  is 

necessary  to  apply  a  convergence  forcer.  The  Gibbs  free  energy  of  the 

system  is  given  by 
vsc 


G  =  n.[p?  +  RT  In  a^] 
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The  procedure  used  for  determining  the  value  of  the  convergence 
forcer,  x,  is  as  follows  [183]: 

1)  Evaluate  the  derivative: 


'^1 

dXj 


vsc 


i.i  ■  jj, 
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where:  a^  *  6^  nl 


If 


''"i  '  "i  ■  "i 

■gx]x=l  -  * 

2)  If  >  0  then  evaluate: 

%  VSC 

f)j=0  =  ,1,  ''"(^>‘1  *  "(1 
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The  optimum  value  of  X  is  then  approximated  by: 


(dG/dX), 


dG 

dG 

[dX, 

1 

O 

11 

[dX^ 

X=1 

Equation  A-15  represents  a  single  step  of  a  Regula-Falsi  root  find¬ 
ing  algorithm  which  yields  a  sufficiently  accurate  estimate  to  the  op¬ 
timum  value  of  the  convergence  forcer.  Infrequently,  equation  A-15  will 
predict  a  negative  value  for  X  and  in  this  case  X  is  set  to  0.05  in  or¬ 
der  to  allow  the  1  iterative  solution  to  continue. 


The  convergence  forcer  is  applied  to  the  updating  of  the  nonbasis 
species  molar  amounts  as 

A-16 

nl  =  n.  + 

The  application  of  the  convergence  forcer  greatly  enhances  the  rate  of 
convergence  of  the  algorithm. 


APPENDIX  B 

GAS  SAMPLING  SYSTEM  AND  MASS  SPECTROMETER 
A  schematic  representation  of  the  mass  spectrometer  and  gas  sampl¬ 
ing  system  coupled  to  the  reaction  tube  used  for  the  decomposition  ex¬ 
periments  is  shown  in  Figure  B-1 .  The  overall  system  consisted  of  a 
three  stage  differentially  pumped  vacuum  system  which  allowed  a  quadru- 
pole  mass  spectrometer  (which  must  be  operated  at  pressure  <^10"^  torr) 
to  be  used  for  the  measurement  of  the  gas  phase  composition  in  the 
reaction  tube  (normally  operated  at  700-1000  torr).  The  objective  of 
the  design  was  to  provide  a  facility  which  would  draw  continuous  gas 
samples  from  the  reaction  tube  in  the  form  of  a  molecular  beam.  The 
formation  of  a  molecular  beam  precludes  the  opportunity  for  intermolecu- 
lar  collisions  or  wall  interactions  while  the  molecules  travel  from  the 
sampling  orifice  to  the  mass  spectrometer  ionizer.  This  condition 
causes  the  composition  of  the  beam  to  be  “frozen  in"  at  the  composition 
present  just  after  the  sampling  orifice  and,  therefore,  can  be  related 
to  the  composition  within  the  reactor.  Due  to  conductance  limitations 
in  the  tube  leading  to  the  skimmer  orifice,  the  pressure  was  too  high  in 
this  tube  to  achieve  the  long  mean  free  paths  required  for  the  formation 
of  a  molecular  beam.  This  Appendix,  therefore,  discusses  the  design 
philosophy  and  operating  parameters  of  the  system  and  suggests  design 
modifications  which  will  result  in  a  functional  modulated  molecular 
beam  facility.  It  should  be  noted  that  even  though  a  molecular  beam 
was  not  present  for  the  decomposition  experiments  in  this  work,  the 
large  signal  intensities  for  AsH^  and  PH^,  combined  with  the  extremely 
small  background  signals  at  these  mass  to  charge  ratios,  allowed  very 
accurate  determinations  of  the  vapor  phase  AsH^  and  PH^  compositions  to 
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be  made.  This  was  not  the  case  for  the  NH^  decomposition  experiments, 
since  a  large  background  due  to  the  fragment  ions  of  H^O  and  N^,  which 
were  desorbing  off  the  vacuum  system  walls,  were  present.  A  modulated 
molecular  beam  facility  would  be  expected  to  greatly  improve  the  results 
of  the  NHj  experiments. 

The  application  of  a  modulated  molecular  beam  mass  spectrometer  for 
the  measurement  of  gas  phase  compositions  in  reactors  operating  at  high 
pressure  (>ltorr)  has  received  considerable  interest  and  two  of  the  most 
complete  publications  regarding  these  systems  have  been  presented  by 
Dun  et  al .  [186]  and  Foxon  et  al .  [187].  Basically,  the  technique  em¬ 
ploys  a  very  small  sampling  orifice  (-0.1  mm  diameter)  across  which 

gases  are  expanded  from  the  high  pressure  reactor  into  the  first  stage 

.3 

of  a  vacuum  system  operating  at  a  pressure  of  approximately  10  torr. 
The  pressure  ratio  across  this  orifice  is  very  high  which  results  in  a 
supersonic  jet  of  gas  in  the  first  vacuum  system  stage  with  the  orifice 
being  at  critical  flow.  A  second  orifice,  known  as  the  skimmer  orifice, 
is  then  aligned  with  this  gas  jet  not  more  than  0.5  mm  downstream  of  the 
sampling  orifice.  The  skimner  orifice,  which  is  frequently  conically 
shaped  and  has  a  diameter  approximately  equal  to  that  of  the  sampling 
orifice,  allows  a  fraction  of  the  gas  jet  to  pass  into  the  second  stage. 
This  skimming  effect  causes  the  angular  velocity  distribution  of  the 
molecules  to  be  highly  peaked  in  the  direction  of  the  mass  spectrometer, 
and  therefore,  the  probability  that  these  molecules  will  collide  with 
the  vacuum  system  walls  and  participate  in  heterogeneous  reactions  is 
small.  Furthermore,  the  pressure  in  the  second  vacuum  stage  is  usually 
on  the  order  of  10“  torr.  The  long  molecular  mean  free  paths  (500  M 
for  N2)  at  this  pressure  preclude  the  chance  for  homogeneous  gas  phase 


reactions  to  occur. 


During  passage  through  the  second  vacuum  stage,  the  molecular  beam 
is  modulated  by  a  motor  driven  chopper  blade.  The  chopped  beam  then 
passes  through  a  collimating  orifice  (nominally  2  nn  in  diameter)  and 
finally  into  the  ultrahigh  vacuum  system  (operating  at  a  pressure  on 
the  order  of  10  torr)  where  it  encounters  the  ionizer  of  the  quadru- 
pole  mass  spectrometer. 

The  signal  from  the  mass  spectrometer  is  sent  to  a  lock-in  ampli¬ 
fier  which  is  locked  to  the  frequency  of  the  chopper.  The  output  of 
the  lock-in  amplifier  therefore  represents  a  frequency  and  phase  corre¬ 
lated  signal  which,  in  principal,  contains  information  only  about 
species  present  in  the  molecular  beam.  In  practice,  however,  molecular 
species  desorbing  from  the  vacuum  system  walls  which  arrive  in  coinci¬ 
dence  with  the  modulated  molecular  beam  are  also  detected.  The  objec¬ 
tive  of  using  a  phase  and  frequency  lock  detection  technique  is  to 
markedly  increase  the  signal  to  noise  ratio  (or  signal  to  background 
ratio)  over  that  which  is  possible  by  an  unmodulated  measurement.  This 
greatly  enhances  the  ability  of  the  instrument  to  make  quantitative 
measurements  down  to  very  low  levels  of  concentration  (1  ppm  or  below). 

The  mass  spectrometer  and  gas  sampling  system  shown  in  Figure  B-1 
was  mated  to  the  60  n»n  quartz  reactor  by  .fusing  a  38  mm  quartz  tube 
to  the  end  of  the  reaction  tube.  The  38  mm  tube  was  then  connected  to 
the  vacuum  system  using  a  compression  type  viton  o-ring  fitting. 

The  skimmer  orifice  was  fabricated  from  10  mm  quartz  tubing.  One 
end  was  drawn  to  a  point  and  then  scribed  and  broken.  The  desired  dia¬ 
meter  of  the  skimmer  orifice  was  obtained  by  sanding  the  drawn  end  of 
the  skimmer  tube  back  until  a  channel  was  clearly  visible  and  a  thin 
stream  of  alcohol  was  able  to  pass  through  it.  A  micrometer  and  lens 
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was  used  to  estimate  the  orifice  diameter  at  0.1  mm,  but  no  method 
of  accurately  determining  the  diameter  was  available.  The  skimmer 
orifice  was  mated  to  the  vacuum  systan  using  a  compression  type  viton 
o-ring  fitting.  Adjustment  of  the  distance  between  the  skimmer  orifice 
and  the  sampling  orifice  was  accomplished  by  varying  the  amount  of 
skimmer  tube  which  was  inserted  into  the  o-ring  fitting.  Alignment  of 
the  skimmer  and  sampling  orifices  was  made  possible  by  two  micrometers 
mounted  through  the  walls  of  the  first  vacuum  stage  resting  at  right 
angles  to  each  other  on  the  skimmer  tube  mount.  The  reaction  tube  was 
filled  with  helium  and  the  position  of  the  skimmer  orifice  was  rastered 
until  a  maximum  in  the  mass  spectrometer  signal  for  He^  (m/e  *  4)  was 
observed. 

A  chopper  motor  cooled  by  a  water  jacket  was  mounted  Inside  the 
second  vacuum  stage.  The  motor  (Globe  model  FC-75A1003-2)  was  designed 
to  operate  from  120  VAC  at  400  HZ.  An  audio  oscillator  driving  a  power 
amplifier  was  used  to  vary  the  speed  of  the  motor.  The  chopper  blade 
was  fabricated  from  29  mm  diameter  brass  by  drilling  a  7  mm  hole  in  the 
base  for  the  motor  shaft  and  cutting  a  slot  19  mm  wide  through  the  brass. 
This  formed  a  chopper  blade  with  a  50*  duty  cycle.  Set-screws  were  used 
to  hold  the  chopper  blade  onto  the  motor  shaft. 

The  UHV  system  was  isolated  from  the  first  two  vacuum  system  stages 
by  a  2  mm  collimating  orifice  and  a  gate  valve.  The  background  pressure 
in  the  UHV  system  when  isolated  from  the  second  vacuum  stage  was  7x10”^^ 


torr . 


The  quadrupole  mass  spectrometer  head  was  an  EAI  QUAD  250B  using  a 
cross  beam  electron  ionizer  and  dynode  style  electron  multiplier.  This 
mass  spectrometer  has  a  sensitivity  of  100  A/ton  for  N2  and  a  resolution 


of  fVAM>2M  up  to  m/e  *  SOO.  The  electronics  used  to  drive  the  mass 
spectrometer  head  were  those  of  a  modified  Finnigan  3000. 

The  operating  parameters  for  the  mass  spectrometer  are  shown  in 
Table  B-1 .  These  parameters  are  in  accord  with  those  suggested  in  the 
Finnigan  operator's  manual  except  that  the  ion  energy  was  increased  to 
7.1  V.  This  value  of  ion  energy  was  chosen  based  on  a  maximum  in  peak 
height  which  was  observed  for  all  m/e  values.  The  peak  height  was  also 
found  to  be  very  sensitive  to  the  emission  current  and  electron  multi¬ 
plier  voltage. 

The  first  vacuum  stage  in  the  system  was  pumped  by  a  60  1/s  roots 
blower  backed  by  a  12  1/s  rotary  vane  pump.  The  pressure  in  this  stage 
was  monitored  using  a  thermal  conductivity  guage. 

The  second  vacuum  stage  was  pumped  by  a  280  1/s  diffusion  pump 
backed  by  a  6  1/s  rotary  vane  pump.  A  liquid  nitrogen  cryogenic  trap 
was  placed  between  the  diffusion  pump  and  the  vacuum  system  in  order  to 
prevent  the  diffusion  pump  oil  from  backstreaming  into  the  vacuum  sys¬ 
tem.  Pressure  in  the  second  vacuum  stage  was  monitored  using  a  tubu¬ 
lated  ion  guage  connected  to  a  Granville  Phillips  series  260  controller. 
A  thermal  conductivity  guage  was  used  to  monitor  the  foreline  pressure 
of  the  diffusion  pump. 

The  UHV  system  was  pumped  by  a  400  1/s  triode  ion  pump.  A  triode 
style  pump  was  chosen  due  to  its  increased  ability  to  pump  inert  gases 
when  compared  to  diode  pumps.  Pressure  In  the  UHV  system  was  monitored 
using  a  nude  ion  guage  connected  to  a  Granville  Phillips  series  271  con¬ 
troller.  A  thermal  conductivity  guage  was  located  in  the  roughing  line 
for  monitoring  the  pressure  during  rough  pumping. 
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Table  B<1 


Mass  Spectrometer  Operating  Parameters 


Ion  Energy:  7.1  eV 
Electron  Energy:  70.0  eV 
Lens  Voltage:  10.3  V 
Extractor  Voltage:  4.3  V 
Collector  Voltage:  26.0  V 
Electron  Multiplier  High  Voltage:  -2.06  kV 
Emission  Current:  1.05  mA 

_7 

Preamplifier  Sensitivity:  10  A/V 
Filter:  300  AMU/s 
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Pressure  tapers  in  each  of  the  vacuum  stages  are  shown  in 
Figure  B-2  for  a  reaction  tube  pressurized  to  760  torr with  Ng  at  300 
K.  These  pressure  tapers  were  determined  by  noting  the  pressures  in¬ 
dicated  by  the  ion  guages  and  first  stage  thermal  conductivity  guage 
and  then  back  calculating  the  tapers  using  standard  conductivity  cal¬ 
culations  [188]. 

The  first  stage  is  seen  to  operate  at  a  pressure  of  0.05  torr  in 

the  region  near  the  sample  orifice  and  0.03  torrat  the  flange  leading 

to  the  roots  blower.  The  second  stage  pressure  at  the  outlet  of  the 

-3  6 

skinner  orifice  is  1.4x10  torrfalling  to  2x10"  torrat  the  end  of 
the  skimmer  tube.  The  high  pressure  region  following  the  skinner 
orifice  is  due  to  the  large  pressure  just  upstream  of  this  orifice 
and  the  poor  conductance  offered  by  the  10  tm  skimmer  tube.  It  is 
this  initial  high  pressure  region  which  makes  the  formation  of  a  mole¬ 
cular  beam  impossible  in  the  present  system.  From  the  kinetic  theory 
of  gases,  an  expression  for  the  mean  free  path  of  nitrogen  at  300  K 
is  [188] 

\  =»  5x10”Vp  B-1 

where:  X  =  mean  free  path  (M) 

P  =  pressure  (torr) 

Thus,  at  0.001  torr  the  mean  free  path  is  on  the  order  of  5  cm  which 
suggests  that  many  intermolecular  collisions  will  occur  in  the  upstream 
section  of  the  skimmer  tube.  These  collisions  will  cause  the  molecules 
to  be  scattered  and  the  molecular  beam  will  be  lost. 

Table  B-2  summarizes  the  pumping  speeds  and  conductances  for  each 
stage  in  the  vacuum  system.  A  severe  conductance  limitation  is  caused 
by  the  use  of  a  1  M  length  of  50  mm  tubing  to  connect  the  roots  blower 
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Figure  B-2 

CurreriL  Vacuum  Stage  Pressure  Tapers 


Vacuum 


to  the  first  stage  of  the  vacuum  system.  Shortening  this  piece  of 
tubing  and  increasing  its  diameter  would  greatly  Increase  the  pumping 
capacity  of  the  first  stage  and  consequently  lower  the  pressure  up¬ 
stream  of  the  skimmer  orifice. 

A  drastic  conductance  limitation  is  also  present  in  the  second 
vacuum  stage  due  to  the  small  diameter  of  the  skimmer  tube.  Increasing 
the  diameter  of  the  skimmer  tube  would  greatly  reduce  the  pressure  near 
the  skiimer  orifice.  This  in  turn  would  increase  the  molecular  mean 
free  path  and  allow  the  formation  of  a  molecular  beam. 

The  following  improvements  are  suggested  for  upgrading  the  first 
twd  vacuum  stages  in  order  to  attain  the  proper  operating  conditions 
for  a  molecular  beam: 

1)  Raise  the  mount  for  the  roots  blower  such  that  a  line 
less  than  0.5  M  long  can  be  used  to  connect  it  to  the 
first  vacuum  stage. 

2)  Increase  the  diameter  of  the  line  between  the  roots  blower 
and  the  first  vacuum  stage  from  50  mm  to  100  mm. 

3)  Increase  the  diameter  of  the  entire  first  vacuum  stage 
to  100  mm. 

4)  Increase  the  diameter  of  the  skimmer  tube  from  10  mm 
to  38  mm. 

A  summary  of  the  predicted  pumping  speeds  and  conductances  in  the 
vacuum  system  incorporating  the  previously  mentioned  improvements  is 
presented  in  Table  B-3.  Figure  B-3  shows  the  predicted  pressure  tapers 
for  the  upgraded  system  under  the  same  conditions  as  Figure  B-2.  Much 
lower  pressures  and  a  smaller  pressure  taper  are  predicted  to  exist  in 
the  skimmer  tube  due  to  the  lower  pressure  at  the  inlet  to  the  skimmer 
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Table  B-3 

Vacuum  System  Pumping  Speed  and  Conductance 
Summary  for  the  Improved  Design 


M 


‘.iV'.'V’,' 


Pumping  Speeds 
At  pump  flange 
At  vacuum  system  flange 
At  orifice  outlet 
Conductances 


1st  Stage 
60. 

33. 


2nd  Stage 
280. 
130. 


UHV  Stage 
400. 
310. 


From  pump  flange  to 
vacuum  system  flange 


From  vacuum  system  flange 
to  orifice 


Orifice 


24. 

17. 

14. 

74. 

240. 

1400. 

92. 

20. 

15. 

3.0003 

>0.005 

-0.4 

Figure  B-3 

Predicted  Vacuum  Stage  Pressure  Tapers  After  Conductance  Upgrade 
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orifice  and  the  increased  conductivity  in  the  skimmer  tube.  The  mean 
free  path  of  a  300  K  N2  molecule  exiting  the  skimmer  orifice  is  expected 
to  be  on  the  order  of  10  M.  Since  the  length  of  the  skimmer  orifice  is 
approximately  0.3  M,  it  is  expected  that  these  conditions  will  be  very 
favorable  for  the  formation  of  a  molecular  beam. 

Foxon  et  al .  [187]  have  developed  a  sophisticated  modulated  mole¬ 
cular  beam  mass  spectrometer  system  and  have  used  it  to  study  the 
effusion  products  of  GaAs,  GaP  and  InP.  They  found  that  a  cryogenic 
baffle  located  directly  behind  the  cross  beam  ionizer  of  the  mass 
spectrometer  greatly  reduced  the  background  signals  which  were  present 
due  to  species  desorbing  from  this  vacuum  system  wall.  The  observation 
of  chemical  species  produced  as  a  result  of  wall  interactions  upstream 
of  the  mass  spectrometer  was  reduced  by  varying  the  chopper  motor  fre¬ 
quency  and  applying  Fourier  transforms  to  the  data  analysis.  A  photo¬ 
diode  was  coupled  to  their  chopper  blade  in  order  to  provide  a  reference 
signal  for  the  lock-in  amplifier  since  the  chopper  motor  was  unable  to 
remain  synchronized  with  the  rapidly  changing  oscillator  frequency. 

Practically,  the  use  of  Fourier  transforms  during  data  analysis 
suggests  that  computer  controlled  data  acquisition  and  mass  spectrometer 
control  should  be  instituted.  Also,  pulse  counting  instead  of  analog 
signal  measurements  should  be  employed.  These  improvements  would  result 
in  the  data  being  immediately  available  to  the  computer  for  the  appli¬ 
cation  of  signal  averaging  and  data  reduction  techniques. 

Finally,  it  is  recommended  that  a  quadrupole  mass  spectrometer 
head  containing  a  Faraday  cup  and  channel  type  electron  multiplier  be 
used.  Electron  multipliers  which  employ  a  channel  structure  yield 
higher  electron  gains  and  the  outputs  are  more  amenable  for  pulse 
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counting  since  all  of  the  output  pulses  are  nearly  of  the  same  ampli¬ 
tude  [189]-  The  presence  of  a  Faraday  cup  allows  the  electron  multi¬ 
plier  gain  to  be  determined  on  a  day-to-day  basis.  This  is  a  very  im¬ 
portant  piece  of  information  if  a  library  of  fragment  ion  patterns  is 
to  be  collected  for  quantitative  investigations. 
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